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VTOL IN GROUND EFFECT FLOWS 
FOR CLOSELY SPACED JETS WITH TEMPERATURE 

W.G. Hill, R.C. Jenkins, S.G. Kalemaris and M.J. Siclari 

Grumman Aerospace Corporation 
Bethpage, New York 

1 - SUMMARY 

The primary purpose of this study was to obtain detailed pressure, 
temperature, «md velocity data for twin- fan configurations in- ground-effect and 
to develop fiow models to aid in predicting pressures and up wash forces on air- 
craft surfaces. For the basic experiments, 49.5mm-diameter jets were used, 
oriented normal to a simulated ground plane, with pressurized, heated air pro- 
viding a jet. The experimental data consisted of: (1) the effect of jet height 
and temperature on the ground, model, and upwash pressures and temper- 
atures, (2) the effect of simulated aircraft surfaces on the isolated flow 
field, (3) the jet -induced forces on a three-dimensional body with various 
strakes, (4) the effects of non-uniform coannular jets. 

For the uniform circular jets, temperature was varied from room 
temperature (24®C) to 232°C. Jet total pressure was varied between 9,300 
Pascals and 31,500 Pascals. For the coannular jets, intended to repr'^-'^t 
turbofan engines, fan temperature was maintained at room temperature while 
core temperature was varied from room temperature to 437®C . Fan and core 
total pressures were chosen to match the power settings most frequently used 
during the large-scale tests of Reference 1. 

In general, the test data correlate well with the induced forces pre- 
sented in Reference 1. No scale effects were found, and jet temperature and 
pressure did not affect the nondimensionalized induced lift. The induced lift 
in-ground-effect was found to be higher for the uniform circular jets than 
for the simulated fan jets. 

Modifications to an existing wall-jet transition model adequately predict 
the trends with height above ground of upwash temperatures and pressures. 
The addition of a recirculation model is necessary to predict upwash tempera- 
tures in the presence of an aircraft. 
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2 - INTRODUCTION 


Lift and control for V/STOL aircraft operating in- ground- effect pre- 
sent a critical condition in sizing the propulsion system . Ground proximity 
effects can restilt in largpe jet-induced lift losses or produce positive fountain 
lift, depending on the aeiropropulsion configuration. The complexity of the re- 
sulting flow field and the sensitivity to many design parameters give rise to a 
large body of experimental data which model various aspects of the flow field. 
The primary purpose of this study was to expand the existing data base for a 
twin- fan V/STOL aircraft by exploring the effects of temperatvire and non- 
uniformities (temperature and pressure) in the jet. These effects were then 
added to an existing computerized prediction method . 

Of almost equal importance , and possibly greater interest to the 
typical reader, this effox^ also investigated scale effects on jet-induced 
characteristics. To this end, the model used was a 1/24-scale replica of a TF 
34-powered, large-scale model tested recently at NASA- Ames (Reference 1). 
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3 - SYMBOLS 


A - Area; nozzle exit area 

b = Thickness of viscous layer 

D = Diameter: nozzle exit diameter 

f = Analytic functions respresenting viscous profiles 

F = Force; radial flux deflection function 

aF = Interference force (total force on body minus thrust) 

h = Nozzle height above ground 

H = Heat; distance to ground 

M = Mach number; momentum 

NPR = Nozzle pressure ratio relative to ambient 

NTR = Nozzle temperattire ratio relative to ambient 

P = Static pressure 

q = Dynamic pressure 

r = Jet radius; ground radius coordinate 

r^ = Ground impingement radius 

R = Radius from one jet impingement point 

5 = Jet spacing 

T = Temperature; static temperature 

Y = Velocity 

X = Distance perpendicular to Y and Z 

Y = Distance from midpoint of line joining jet centerline 

Z = Distance above ground; vertical jet or ground coordinate 

Z' = Distance downward from nozzle exit 

3 = Exponent of profile function 

t = Isentropic exponent; ratio of specific heats 

6 = Boundary layer thickness 

A = Difference between local and ambient conditions 
X, - Non-dimensional wall thickness 

n = Thrust efficiency factor; nondimensional thickness 
0 = Angular orientation around one jet impingement point 

p = Density 

Y = Angle in ground polar 
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a, A -- Ambient 

B = Body 

C = Core nozzle exit conditions; potential core 
F = Flux, fan nozzle exit conditions 

FD = Fully-developed 

g = Value at ground effect height 

H = Half value 

J = Nozzle value; jet exit conditions for open circixlar nozzle 

m, M = Maximum value; viscous layer 

N = Nozzle 

PC = Potential core 

rec = Recirculation 

S = Stagnation value; static 

T = Thermal layer; total; stagnation 

u = Upwash 

V = Velocity layer 

W = Wall jet 
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4 - EXPERIMENTAL BACKGROUND 


4. 1 Test Facility 

The experiments conducted under this contract were performed in the 
Grumman Environmental Test Facility (ETF) . This facility performs a wide 
rcuige of system and component testing for Grumman and has two compressor 
systems capable of supplying heated air: the ram air system for pressures to 
110 kPag and temperatures to 232^C , and the bleed air system for pressures to 
3.45 MPag and temperatures to 760°C. A general overview of the installation in 
this facility is given in Figure 1; a detailed photo of the test arrangfement is 
presented in Figure 2. 


4.2 Air Supply 

For most of the tests, the air was compressed by a Rootes-type blower 
(referred to in this facility as the "ram air compressor"). Compression heating 
of the flow was removed and the flow divided into two piping systems. One 
system was then heated by a controlled resistance electrical heating system with 
a maximum temperattire of 249°C . Both the hot and cold flows were piped to the 
test site by 20.3-cm diameter insulated piping systems. Supply pressure was 
controUed by a feedback control system at the compressor control panel and a 
second control valve at the user's outlet area (Figure 3) . For the experiments 
with core temperatures above 232®C , air was supplied by a piston-type com- 
pressor and gas and electric heaters in series (bleed-air systems) . 

4.3 Settling Chamber 

The settling chamber for the fan nozzle was a cylindrical steel tank of 
76.2-cm diameter and 218.4-cm length, adapted from underwater flow research. 
(See Figure 4). The first 91.4 cm of this chamber are occupied by a diffuser to 
spread the flow from the 20. 3-cm diameter inlet evenly over the chamber cross- 
section. This is followed by a tube bundle of 9.5-cm diameter stainless steel 
tubes and two fine mesh screens. The heated air supply for the core flow enters 
through the bottom of the chamber after a combined turn and transition from 
a single 20.3-cm pipe to a pair of separate 5.1-cm tubes. This supply is con- 


5 



©©©(D00Q©©©©©© 


LEGEND 


DISCONNECT SW. - COMPRESSOR 
SW. GEAR - COMPRESSOR SW. 

DISCONNECT SW, - HEATER 
DISCONNECT SW. - CONDENSER PUMP 
CONDENSER PUMP SW. 

DISCONNECT SW. - REFRIG. UNIT 
REFRIG. UNIT SW. 

DISCONNECT SW. - EXHAUST FAN (CELL #3) 
MOORE VALVE CONTROLLER 
FUSED SW. 

CONTROL POWER SW. 

EMERGENCY STOP SW. 

HEATER CONTROL POWER SW. 




Cell # 1 
Controli 


© 

PL #26 KEY SW. 


HEATER POWER 
ON-OFF BUTTONS 


T«st Call # 1 


Air Drop 
(Typ, 6 Places) 


0 

/ 


R81-1622-O01D 


I Aftarcoo<»'r 



Compfisor 

Room 


I Metar ) © 

I SCR I 






1 


Test Cell # 3 


CD© 





Main 

Power 

Panel 






Cell # 3 Controls 


Ground 

Plane 

Assembly 


Test Area # 2 


Figure 1 - Ram air system location plan, Plant 31. 
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tinued inside the chamber by flexible, insulated pipes to the nozzle entrances . 
For the open- jet experiments , this tubing system was entirely removed and the 
piping system was changed to bring heated air in through the diffuser inlet . 


The outer fan nozzle, which is also the hardware fcr the open-circular- 
jei experiments, is an ASME elliptical profile flow metering contour with an in- 
side diameter of 5,0 cm, (Figure 5). This nozzle is continued as a straight 
tube for 30.5 cm in a tradeoff of boundary layer growth and heat transfer in the 
nozzle flow against the need to separate the model from the front face of the 
chamber, thus preventing the chamber frmi interfering with the recirculation 
flow field. The inner core nozzle for the heated flow begins with a contraction 
from the 5.1-cm flexible supply tubes beginning 5.1 cm before the fan flow 
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Figure 3 - Ram air system sch«natic. 
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nozzle entrance. This contraction connects to a double wall tubing system with 
an outer diameter of 3.8 cm (the scale centerbody diameter) and an inner dia- 
meter of 3.0 cm. The area between the two walls of this contraction section and 
the constant area section running to the final core nozzle contraction were in- 
sulated with aliuninum oxide powder. The final contraction to the 20.7-mm core 
exit diameter and the core centerbody outer profile are accomplished in one 
solid machined segment. The core-flow tube is held in place by a machine screw 
in the simulated core pylon, and three thin finger-type supports near the 
nozzle entrance. 

4.5 Ground Plcne 

The simulated ground plane consists of a 2.5-cm thick aluminum plate 
121.9-cm square, with a circular insert of steel for the pressure tap and wall 
jet probe movmtings. (Figure 6a). This plate was mounted on a framework 
which rode on rails that were aligned parallel to the nozzles. The ground plane 
was moved by a screw that was driven by a variable speed motor /controller. 

The instrumentation locations are shown in Figure 6b, c, and d, (pressure 
taps, surface temperature thermocouples and wall jet rakes). Numbering of 
the instrumentation locations matches that of the full-scale ground plane 
(Ref. 1), with additional instrumentation designated by "R”. 

4.6 Aircraft Model 

The model (Figure 7) was a 1/24 scale reproduction of the Grumman 
Design 698-411B tested in the NASA-Navy-Grumman full-scale demonstrator 
program. This model was constructed of a partially-hollow high-temperature- 
plastic fuselage section housing the balance and pressure /temperature lines, a 
metal bottom plate for pressure and temperature instrximentation, metal strakes 
at the fuselage chines, and separately attachable plastic wings and empennage. 

4.7 Instrumentation 

4.7.1 Pressure. - Pressures were measured with Validyne model DP15 or DP103 
transducers, powered by single channel CD-15 or MCl-lO multichannel systems 
using CD 18 or CD 19 carrier demodulators. Transducers were calibrated versus 
an oil filled manometer, either singly (transducers used for probe measurer. cnts) 
or in groups with the manometer calibration pressure applied through a manifold 
(ground plane or aircraft surface). 
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relocated so that temperature (thermocouple) instrximentation could be added. 

4.7.? Temperature . - All temperature measurements were made with chromel- 
altunel (Type k) thermocouples connected by a multi-connection switch box to 
either an Omega model 2175A or a Fluke model 2166A linearized thermocouple 
readout amplifier. The ground plane surface temperatures and those on the 
flat segment of the aircraft bottom were measured with thermocouples spot 
welded to thin (0.4 mm) stainless steel areas. On the aircraft, the entire flat 
bottom was formed from one thin piece of stainless steel, and thicker pieces 
were spot-welded in areas for pressvure tap installation. (Figure 8). The 
30.5-cm diameter ground plane insert was 4. 2-cm thick, and the areas for 
temperature measurement were machined to the 0. 4-mm thickness for a distance 
of at least 1.3 cm from all thermocouples. (See Figure 6c). The thin skin 
allows gfreater spatial resolution in temperature measurements , which is impor- 
tant in small-scale testing. 

4.7.3 Probes . - Three separate probe assemblies were used for measurements 
in the free jet, wall jet, and upwash flows (Figure 9). In each probe assembly 
the temperature data were taken by a bare thermocouple probe protruding from 
a 1.5-mm diameter stainless steel sheath that was located close to one of the 
pressure probes. Bare thermocouples were chosen, rather than a more con- 
ventional shielded thermocouple probe, in order to obtain good spatial reso- 
lution and time response and to avoid flow disturbances that would be asso- 
ciated with a housed probe. Temperature measurements with bare thermocouples 
require correction to account for radiation losses. Temperature corrections 
were acquired for the thermocouple probes in each assembly by recording output 
in a flow field of known temperature over a wide temperature range. We used 
the open jet nozzle configuration for this work. Each probe assembly was 
placed on the centerline of one nozzle shortly downstream of the exit and a 
housed thermocouple probe (United Sensor t 3 rpe KT-12-C/A) was placed at a 
corresponding point in the flow from the other nozzle. 

The probe assembly used for free-jet data (Fig. 9a) consisted of 1.6-mm 
O.D. sharp -lipped pitot tube with a 0.5-mm diameter thermocouple displaced 
2.5 mm laterally from the tip of the pitot tube. Wall jet data were taken with a 
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1.3-mm O.D. pitot tube (Fig. 9b) that was flattened to 9. 5 x 1. 5-mm dimensions 
at the tip, with a 0.5-mm diameter thermocouple that was displaced laterally 
1.8 mm. The probe assemb’y Tor upwash measurements consisted of a 3.0-mm 
O.D. Kiel probe with a 1.6-mm O.D. static tube displaced 4.3 mm on one side, 
and 0.5-mm diameter thermocouple displaced 3.0 mm on the other side of the Kiel 
probe (Figure 9c). 

The free-jet probe assembly and the upwash probe assembly were 
supported by a motor-driven traverse that was mounted on the same rails as 
the ground plane carriage. Probes were traversed in the plane containing the 
two jet centerlines, along paths perpendicular to those centerlines. 

The support for the wall-jet probe assembly passed through a hole in 
the ground plane and attached to a manually driven screw tra^'-erse that was 
fastened to the back of the ground plane (Figure 10). The probe location re- 
lative to the wall was measured from the voltage output of a thin film linear re- 
sistance transducer that was attached to the screw drive. 

4. 8 Experimental Techniques 

All experiments were conducted under steady or pseudo-steady 
operating conditions. After initial startup time (approximately one hour) 
changes in pressure operating conditions took only a few minutes. Total 
temperature changes were much slower, primarily because of the heat capacity 
of the supply piping. 

Under most of the operating conditions for this contract, we utilized 
the techniques we developed earlier (Refs. 2, 3 and 4) of continuously recording 
test variables versus input parameters in an analog form on an X-Y recorder. 
This allowed a constant monitoring of test results. Operating conditions could 
be changed to allow for results observed (such as gradients or fluctuations). 

The ; .ost frequent examples of this are the recording of interference force, 
ground-plane pressures, and vehicle surface pressures versus height above 
ground, or probe temperature and pressures versus probe position. 

Probe and ground plane position were converted to voltage by a cable- 
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driver. I’otary potentioneter (Figure 11) for most cases, and a linear potent- 
iometer for the wali-jet measurements (Figure 10). Voltage was supplied by 
Hc'.vlotx Packard Model 6215A D C. power supplies and adjusted so that a 
"calibrs.ed” scale setting ‘jouid be used on the recorder (sllcwing changes in 
scale with ease) . Position was controlled by a variable speed DC mctor/cont- 
troUer and a screw drive. 

For the aspects of the experiments requiring mariy mefiisurements 
simultaneously (such as ground plane 6uid aircraft surface pressures) data were 
recorded simultaneously by an on-site minicomputer (Hewlett-Packard Model 
lOOO(X) through A/D converters. 

The aircrai force measurements were taken on separate runs fro.m the 
pressure and temperature measurements because of force interference from the 
pressure and te.nperature leads which could not be eliminated. The instru- 
mented lower surface and the instrumentat’on leads attached to it were removed 
for the force tests. A plain lower surface was fabricated to match the external 
contours of the instrumented lower surface and used to cover the cavity in 
the fuselage during the force tests. 
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5.1 Fan- Jet Nozzles 

The sketch of fan-jet impingement geometry in Fig. 12 illustrates the 
co-ordinate system used to describe probe and nozzle locations. All data were 
taken with the spacing between nozzles scaled to the full-scale model of 
Reference 1, with the jets normal to the ground and with the aircraft model 
present . 

5.1.1 Nozzle Flow . - Figures 13 and 14 show pressure profiles taken across the 
fan and nozzle exit, with baseline pressure contritions and a 24°C total tempera- 
ture in both flows. The total-pressure profile at the core exit is quite flat. The 
profile across the fan exit shows a well-defined maximum near the center of the 
annulus. The maximum value of the pitot pressure at the exits was used to 
define the nozzle pressure ratios (1.103 for the core and 1.195 for the fan at 
baseline conditons) . Most of the fan-jet impingement data were taken with these 
nozzle pressure ratios, chosen to match those conditions most frequently used in 
the full-scale experiments. After the initial series of experiments, during 
attempts to compare full-scale and model results, it was discovered that the 
full-scale conditions tc be matched were based on area averaged total pressures . 
Additional tests were conducted at the most critical conditions with the correct 
pressures. Non-dimen.sional results were found to have very little dependence 
on the total pressure. 

Integration of the pressure profiles in Figures 13 and 14 (assuming 
axial symmetry) gives a value of 36.5 N thrust for one fan-jet nozzle. However, 
this approach provides only an approximation of the thrust because the fan 
flow is not truly eudally symmetric (engine support pylon blockage on one side 
of the fan duct) and because this method of determining thrust does not account 
for pressure forces on the outer cowl of the core nozzle (centerbody) . 

A more accurate determination of nozzle thrust was made by measuring 
the impingement on a 0.914-meter square ground-plane plate that was instru- 
menteo with three strain beams. With baseline pressure conditions, the thrust 
produced by one nozzle was found to be 36.2 N. Maintaining pressure within 5% 
of the baseline conditions , the core exit temperature was raised from 24° to 
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Figxire 13 - Pitot pressure profile at fan nozzle exit 
baseline pressure operation. 
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168°C, As shown in Figure 15, there was a negligible variation in nozzle thrust 
over this range of core exit temperatures when minor variations in total pre- 
ssure at the core exit and the fan exit were taken into account. 


The thrust of a fan-jet nozzle was measured for a range of both core 

and fan nozzle total pressures. The thrust was found to be a Unear function of 

fan and core pressures (Figure 16) . This is true over a range of pressures 

exceeding those used in this study. An empirical fit of the data yields 

F = 1.55 q - + 0.719 q -0.445, in which q is in kilopascals and F 

^ fan, avg. ^ core ^ 

is in Newtons. As an indication of the origin of this equation, the following 

relationship was assumed: 


F=/K, +K_\ 2q^ A- + K„2q A 

\ 1 2 ) ^fan fan 3 ^core core 


( 2 ) 


where and represent the relative thrust efficiencies and K 2 refers to the 

drag of the core centerbody in the fan flow. For the nozzles used in the study, 

2 2 
A. =7.8 cm and A =3.37 cm 
fan core 

From integration of the core nozzle exit profiles, was found to be 0.0962. 
Combining the tenns relative to the core, + ^ 2 ^ x 2 A^^^ = 
we have the relationship 

F = K . X q- + 0.650 q (3) 

4 ^fan ^core. 

Solving for with room-temperature, baseline condition data yields = 1.46. 


Substituting, we find 

AF = 1.46 q^ + 0.650 q 

^fan, avg. ^core, avg 

This agrees quite closely with the empirical relation of Eq. (1) 


(4) 


Figure 17 shows free-jet total-pressure profiles that were taken (with 
both fan jet nozzles operating) at successive axial stations for a total temperature 
at the exit of 24°C. The lack of symmetry in each of the fan flows can be 
attributed to the core pylon blockage in the fan nozzles. The apparent 
convergence of the two flows, indicated by measurement on these profiles of 
center-to-center distance between the two flows at several axial stations, was 
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Figure 16 - Effect of fan and core 
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Figure 17 - Free- jet. pitot pressure profiles, baseline operation. 


ORIGINAL PAGE IS 
OF POOR QUAUTY 


not caused by actual convergence of eith'^r the fan or core nozzle hardware. 

The reason for this characteristic of the free-jet pitot pressure profiles is still 
under active inv'^estigation , but it is assumed to be due to the pylon-induced 
asymmetries. 

5.1.2 Ground Flow. - Pressures taken on the ground surface were recorded by 
a 30-channel A/D system connected to an HP- 1000 minicomputer as the ground 
height was varied. In addition, seme pressure data were taken at fixed ground 
heights 6is continuous profiles by displacing the ground plane in the y-direction. 
Figure 18 shows the variation of ground pressure along a line between the two 
jet impingement points for three body heights above ground. These ground 
pressure profiles were taken using baseline nozzle pressure conditions and 24°C 
core nozzle exit temperature. The maximum ground pressure in the stagnation 
line region (central section of the profiles) shows little variation with body 
height. The ground pressure in the impingement regions (outer sections of the 
profiles) shows an inward shift of the ground impingement center for each jet 
as ground height is increased, which is analogous to the convergence of the 
free jet centerlines illustrated by Fig^ure 17. 

Ground temperature measurements were taken at points around one jet 
impingement region for the same three heights above ground, using almost the 
same nozzle pressure conditions and a core exit temperature of 425°C. Figure 
19 shows the ground temperature radial variation from one jet impingement 
point for three different angular orientations ( 9 ) . A curve was faired through 
the round sjrmbols, corresponding to 9 = 90° (oriented in the wing spanwise 
direction). Note that, for 9 = 0°, the ground temperature was lower, and for 
9 = 180° the temperatures were higher, than this curve. The higher ground 
temperature under the wing may have been caused by blockage of air entrain- 
ment into the jet exhaust by the presence of the wing. 

Figure 20 shows representative wall-jet profiles of total pressure and 
total temperature that were obtained by traversing a probe normal to the ground 
at the ground location scaled to rake nvunber 6. Data were taken at the three 
gpround heights corresponding to those run during the full-scale tests. Data 
were obtained at this ground location and at two other locations f rakes number 
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8 and 10) . Figure 21 shews a tabulation of representative features of the pres- 
sure and temperature profiles taken at ail three ground locations and three 
ground heights. In this figure, the maximum wall-jet total pressures have been 
normalized by the maximum total pressure at the fan exit and the wall-jet thick- 
ness (at one quarter the maximum pirofile dynan»5c pressure) has been normal- 
ized by fan-nozzle outside exit diameter. 

The temperature data shown in Figure 20 were taken with a bare 
thermocouple that was located adjacent to the total pressure tube. These data 
must be corrected to account for radiation losses. The correction has been 
applied tc the data listed in Figure 21. 

5. 1.3 Up wash Flow . - Flow properties in the upwash were measured by traver- 
sing a probe assembly through the centerline of the upwash in the y-direction 
with the model in place. The total pressui-e, static pressure and total temper- 
ature profiles were measured along a liiio .56D below the model lower surface, 
using baseline nozzle pressure conditions and core exit temperatures of 24°C, 
128°C , and 192*^C . Data could not be taken at the lowest body height tested 
full-scale (1.07D) because the probe would have been in the upwash formation 
region. For the greatest full-scale body height (4.59D), the upwash was not 
recorded because the level was too low. Data were taken at the intermediate 
full-scale body height (1.54D), Figure 22, and at a somewhat higher body 
height (1.99D), Figure 23. 

Figures 22 and 23 show the variation of total pressure and of the 
difference between total and static pressures across the upwash for the lowest 
and highest core temperatures. It was found that the upwash pressure pro- 
files were well-centered around the midpoint between the nozzle centerlines and 
that they were essentially independent of core exit temperature. The tempera- 
ture across the upwash was almost constant and can be represented by the 
centerline values plotted in Figure 24. 

5.1.4 Model Forces . - Ground-induced interference forces were meastired for 
both the 15^ and 45^ strake angle configuration for several nozzle pressures and 
temperatures. The first example, varying fan total pressure with core temper- 
ature equal to fan temperature (approximately 24°C), is shown in Figure 25. 
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Figure 21 - Data from wall jet profiles resulting from fan 
jet impingement. 

For both the 15° and 45° strakes, changing fan total preosure with core total 
pressure held constant produced a progression of interference forces with no 
change in the qualitative behavior with height above ground. Note that, for 
all these cases, the fan total pressure was well above the core total pressure. 
The behavior may be changed if the core total pressure becomes equal to or 
larger than the fan total pressure because of changes in the impinging jet 
behavior. 

A more meaningful view of the effects of the fan total pressure can be 
found in the non-dimensional presentation with interference force normalized by 
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Figure 24 - Variation of centerline upwash temperature with ground height . 


39 



the nozzle thrust shown in Figure 26. The total pressure level has no significant 
effect on the non-dimensional interference forces. Full scale data (from Ref. 1) 
are shown on this plot for reference. 

At a fan average total pressure of 19.8 kPa gauge (and core total pres- 
sure, again, at 10.4 kPa gauge) five core total temperatures were run: 43°C, 
104°C, 254*^C, 322°C and 416°C. As shown in Figure 27, there is essentially no 
effect of temperature on the interference forces. Experiments were also con- 
ducted with the nozzle exit flow more closely matching the full-scale engine 
conditions. For a fan total pressure of 25 kPa gauge, model forces were mea- 
sured with core exit temperature of 421^C and 43^C. Figure 27 shows these 
results compared to the full-scale data (Ref. 1). Again, there is a very small 
effect of temperature on interference forces. The close agreement between the 
model test results and the full-scale data indicates that good simulation of air- 
craft ground interference forces can be obtained with a l/24th scale model even 
without core temperature simulation. 

5.1.5 Model Surface Pressures and Temperatures. - Pressure and temperature 
measurements were taken on the model undersurface with the baseline nozzle 
pressure conditions and core nozzle exit temperatures of 128°C and 192°C. 

Figure 28 shov/s dimensionless temperature distributions along the model under- 
side (x-direction) for three heights above ground. The temperature across the 
model surface (y-direction) showed little variation. Values of the surface 
temperature coefficient ((T-T ^/(T -T )) at station number 22 (0.54D aft of 

the naceUe centerline) are very close to vedues found with a probe on the up- 
wash centerline at 0.56D below the model surface (Figure 24). This comparison 
provides some verification for the model surface temperature data, which were 
taken by thermocouples attached to the inner surface of the metallic underside 
of the model as described in Subsection 4.7.2. Figi’re 29 shows the variation 
of model temperature and pressure with body height above ground. The temper- 
ature was measured at station number 22 and the pressure at the midpoint be- 
tween nozzles. 

Model pressure and temperature data were also taken for a higher 
core temperature and a higher fan pressure that more closely matched the full- 
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Figure 26 - Effect of fan nozzle pressure on interference forces 
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scale tests. Figures ^0 and 31 show pressure and temperature measurements 
taken at three heights above ground for a core exit temperature ri 427°C. The 
relatively low model surface temperatures shown at X/Dp = 0.5 in Figure 31 are 

not due to faulty thermocouples. Temperatures on both sides (y-direction) of 
this centerline thermocouple showed the same values, which verified that the 
temperature profile alonj the model had a local minimum slightly forward of the 
nacelle centerline. The magnitude of this local minimum varied with height 
above ground. This distortion did not appear in the data taken at lower core 
temperature (Figure 28) because this centerline thermocouple and one of the 
side thermocouples were inoperative during the earlier tests. 

Comparison of the data shown in Figure 31 with those shown in Figure 
28 indicates that the model surface temperature does not scale well with wide 
variations in core exit temperature. The model surface te.*nperature coefficient 
decreased by more than a factor of iwo when the core flow temperature was 
raised from 128°C to 421*^C. The same poor temperature scaling appeared in 
probe temperature measurements on the upwash centerUnc (Figure 24). In 
this case, the local temperature coefficient in the flow decreased by 15 percent 
when the core exit temperatxire was raised from 128°C *o 192°C . 

Probe temperature surveys were made in the area corresponding to the 
inlets on the fuU-scaie model (l.SD^ above the model upper surface, close to the 

outer surface of the fan nozzles). It was found that the temperature in this 
region was ver«' close to ambient, and fluctuating greatly. The maximvim 
temperature rise observed for a nozzle total temperature of 192®C was approxi- 
mately 2°C. The strakes app.u*ently turn the flow enough to prevent any 

significant amount of upflow from reaching the inlet area. The inlet suction 
may affect the flow field to alter this behavior, but this is not expected to 
happen . 

5.2 Open Circular Nozzles 

5.2.1 Nozzle Flow . - Free jet data were taken with the ground plane removed, 
using settling chamber temperatures from eirabienl to 232°C, and a settling 
chamber pressure of 27.6 kPa above ambient. Total pressure and temperature 
profiles were recorded for probe traverses across both jet centerlines for dis- 
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Figure 30a - Pressure distribution along model. 
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Figxire 30b - Pressure distribution uc/oss model. 
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Figure 31 - Model temperature distribution. 
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tances (Z') from the jet exists up to 15 diameters. Figure 32 shows pitot pressure 
profiles obtained with a settling chamber temperature of 227°C . Profiles taken 
close to the exit exhibit a flat-topped profile with a 1 percent ripple that was 
caused by fluctuations with time of the air supply pressure. Note that there is 
no inward shift of the maxirautr. pitot pressure points at large values of Z' as 
orjurred in the free-stream data obtained with fan- jet nozzles. The corres- 
^ondirg teniperature profiles are shown in Figure 33. These temperature data 
wore plotted directly from the thermocouple output and have not been corrected 
for radiation losses. (This correction was typically of the order of 10*^C.) The 
pressure profiles begin to merge about seven diameters from the nozzle exit. 

The temperature profiles show a greater radial spreading and begin to merge 
about five diameters from the nozzle exit. 

Figure 34 shows data taken from profiles that were obtained with 
different settling chamber temperatures. The data show the decay of pitot 
pressure and temperature along one jet centerline with two jets in operation. 

The centerline pitot pressure decay is slightly greater for heated than for un- 
heated jets. Several runs made with one nozzle blocked off showed no differ- 
ence between the single- and dual-jet results in the centerline values of both 
pitot pressure and temperature out to V - 15 diameters. Hence, the decay in 
centerline properties obtained from two jet operation can be used to represent 
single- jet operation. 

5.2.2 Ground Flow without Model . - Ground-pressure and wall-jet profiles 
were taken at three values of H/D corresponding to the full-scale test 
conditions. Ground pressure measurements were taken without the model to 
determine how closely the radial pressure distribution around one jet impinge- 
ment point matched existing single-jet impingement data. Figure 35 shows the 
data taken along radial paths at two orientations around one jet centerline. 

The solid curve, taken from Ref. 2, represents an empirical approximation for 
the ground pressure distribution for a single jet with H/D between 2.0 emd 
5.0. 

Wall- jet profiles (similar to those shown in Figure 20) were taken at 
the same values of H/D and the same ground locations used in the full-scale 
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Fi^re 35 - Radial ground pressure distribution 
around one jet impingement center. 


fan-jet tests. These measurements, summarized by plots of maximum total 
pressure aiid thickness in Figure 36, show little influence of temperature on the 
wall- jet pressure data. 

This figure suggests that the waU-jet flow was not axially symmetric 
around one jet impingement point. Wall-jet profiles taken at the same radius 
from the closest impingement center show that the flow directed toward the 
stagnation line has a higher peak total pressure and a greater thickness than 
the flow directed away from it. This trend occurs for all tested heights above 
the ground and all settling chamber temperatures. This apparent lack of 
symmetry in the wall jet data around one impinging jet is probably a real dis- 
tortion of the ground flow caused by the impingement of the second jet, but may 
have been caused by interference of the probe support with the upwash 
formation region when it was located at Station No. 10 (see sketch in Figure 36). 
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Figure 36 - Data from wall jet profiles, open jet impingement. 
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5.2.3 Up wash Flow without Model. - Measurements were taken in the upwash 
without the model to evaluate flow properties in the absence of model blockage 
effects. The probes were traversed through the upwash at the nozzle exit 
plane for several values of H/D and a range of settling chamber temperatures. 
The results are summarized by the centerline values of upwash properties at the 
nozzle exit plane that are plotted in Figure 37. The data obtained from the 
pressure probes were essentially unaffected by temperature. Temperature pro- 
files across the upwash were almost flat, and the exit plane temperatures show 
little variation with gpround height cut to H/D = 4. 

Additional upweish surveys were made at a fixed nozzle height above 
ground by traversing the probe through the upwash at various probe heights 
above ground. The results are summarized in Figure 38 which shows the 
variation of temperature and pressure edong the upwash centerline for H/D = 3. 
These pressure plots also show little influence of jet exit temperature. 

The effect of settling chamber pressure on upwash properties with jet 
temperature equal to ambient was investigated by conducting exit-plane probe 
traverses at various ground heights. Fig^ure 39 shows that the centerline 
properties at the nozzle exit plane, for ground heights of 2-1/2 nozzle dia- 
meters or greater, vary linearly with settling chamber pressure over a wide 
range. Furthermore, these measurements agree with past data taken at this 
jet spacing in our low-pressure facility. The lack of pressure sealing at low 
ground heights may be caused by the influence of the ground interaction flow 
on nozzle exit conditions. As bix)ught out in Ref. 3, the pressure distribution 
inside the nozzles close to the exit (for this jet spacing) becomes non-uniform 
when the ground is brought closer than 2-1/2 diameters to the nozzle exit. 

The magnitude of the upwash pressure data shown in Figure 39 
illustrates that the upwash flow at the nozzle exit plane can be treated as an 
incompressible flow for all of the nozzle exit conditions and all of the ground 
heights, even though the nozzle exit flow itself was in the compressible flow 
range. Flow with a dynamic pressure Ipss than about 5 kPa (corresponding to 
airflow velocity less than 90 m/sec) can be considered incompressible without 
a loss in accuracy. Below this value, the dynamic pressure can be found from 
the difference between total and static pressures. Dynamic pressure much 
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F'iguie 38 - Upwash properties between nozzle exit plane and ground, H /D = 3. 
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Fi^re 39 • Effect of settling chamber pressure on upwash properties 
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higher than this (somewhat arbitary) 5kPa level requ.res a compressible flow 
approach to determine dynanuc pressure, since it cannot be found by sub- 
tracting static pressure from total pressure. Nozzle exit total pressure was 
greater than 5 kPa for almost all of our experimental work. 

Therefore, for all of the data shown in Figure 39, the quantity (P.j,-Pg) 

represents the upwash dynamic pressure on the centerline at the nozzle exit 
plane, and the flow velocity corresponding to this dynamic pressure can be 
evaluated with incompressible flow relations. Total pressures were measured 
with a Kiel probe and static pressures with a static probe aligned with the 
mean flow direction (probe dimensions are given in Subsection 4.7.3). Tra- 
versing both probes through the flow simultaneously ( with fixed inter-probe 
distsince, Y’) provided direct plots of the total and static pressure pi-ofiles. 

When overlaid and shifted a distance, Y', one can trace a profile of (P.j,-Pg). 

Data shown as Figure 3'", 38 and 39 were all obtained from the 

peak value of such tracings. 

5.2.4 Effect of Model on Ground Flow. - Ground pressure measurements taken 
with and without the model in place showed that the model had little influence 
on the pressure close to the stagnation line but did change the pressure in the 
jet impingement region. Figure 40 shoivs the distribution of ground pressure 
in one of the jet impingement regions along a line connecting the two jet center- 
lines. The solid curve represents an empirical approximation for the ground 
pressure distribution for single jet impingement. Comparison with data taken 
at the same locations without the model (Figure 24) shows that the model re- 
duces the pressure in this part of the impingement region. Data in Figure 40 
for different jet temperatures show no significant temperature effect. 

5.2.5 Effect of Model on Upwash. - Probe traverses were taken across the up- 
wasa with and without the model in place to determine the influence of the 
model on upwash properties. The temperature profiles across the upwash pre- 
sented in Figure 41 show that the presence of the model raises the temperature 
of the upwash approaching the model. The increased upwash temperature with 
the m del in place indicates a recirculating flow entrainment caused by model 
blockage . 
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Figure 42 shows pressure profiles taken across the upwash with and 
without the model for the same nozzle height above ground. These profiles were 
taken far enough from the model undersurface so that the increased upwash 
dynamic pressure caused by the presence of the model was not a local flow dis- 
turbance involving flow around the model. The dynamic pressure (and hence 
flow velocity) throughout the entire upwash was affected by model blockage. 

This change is significant for flow modelling. The cause of this change is not 
understood at present. The data in Figure 42 were taken with a let exit 
temperature of 238°C. Data taken at temperatures of IIG^C and 24°C showed 
almost identical profiles with and without the model, again illustrating that jet 
exit temperature has little influence on the pressure measurements in the up- 
wash. 


5.2.6 Model Forces. - The thrust of one of the open circular nozzles was 
determined from the relation 


Thrust = 2 q r; A , 

where q is the dynamic pressure associated with the nozzle pressure ratio, A is 
the nozzle exit area and ri is a thrust efficiency factor that accounts for 
momentum loss in the nozzle boundary layer. Based on thrust messureracnts 
made earlier on this 5-cm diameter nozzle, we found r, = .90 for the open cir - 
cular nozzles. _ _ 


Using q =-^PM^ and = 5 
2 




/ Prr.'V 


- 1 


for air. 


the thrust developed by one circiilar nozzle can be found as a function of 
(settling chamber pressure): 

Thrust = 6.30 

Force data taken with the open circular nozzles were normalized by twice the 
value of nozzle thrust. 



Forces were measured on the model with 15° strakes for a range of 
settling-chamber temperatures and pressures. Figure 43 shows the variation of 
model interference force with height above ground for a settling chamber pres- 
sure of 27.6kPa gauge and jet temperatures of 24°C, 116°C and 240°C. These 
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three curves are quite similar, having small differences that appear to show a 
progression with temperature. Further investigation showed that differences 
in the force curves in Figure 43 are more likely to represent scatter attributable 
to slight differences in model alignment between runs. Subsequent measurements 
were taken at a fixed height above ground (H/D = 3.2) as the settling chamber 
temperature weis rsdsed from 24° to 221°C. Testing in this manner, essentially 
no change was found in interference force with jet exit temperature. The non- 
dimensional forces for the fan-jet simulation (Figure 27) are also shown and are 
much lower than those of the present section. 

Figure 44 shows the variation of model force with height above ground 
for different settUnj chamber pressures and a jet ex: 1 nperature of 24°C. 

While these curves are qualitatively the same, they - me variation with 

settling chamber pressure. Ihese data were all tr v at disturbing align- 

ment between tests. 
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Figxire 44 - Effect of nozzle stagnation pressure 
on ground interference forces. 



5.2.7 Model Surface Pressure and Temperature . - Model surface pressures 
were measured with and without the 15° strakes for 24°C nozzle exit temper- 
atures. At higher jet temperatures, the data were taken only with the 
strakes attached. Figure 45 shows the effect of the strakes on the pressure 
distribution along the fuselage undersurface. The pressure distribution across 
the model is shown in Figure 46. The presence of the strakes raises the sur- 
face pressure level and flattens the profile across the undersurface. Figure 
'*7 shows the variation with height above ground of model pressure at three 
centerUna taps for the model with strakes. The pressure on the fan nacelle 
centerline (from Figure 29) is shown for comparison, and is much lower. 

Figure 48 shows that the pressures along the model at a fixed height above 
ground are not significantly affected by jet temperature. 

The temperature distribution along the model is shown in Figure 49. 
For H/D = 1.54, the temperature coefficient at a point between the nozzles was 
approximately 0.8, which is considerably higher than values measured with 
probes (range of 0.5) on the upwash centerline without the model (see Figures 
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Figure 45 - Pressure distribution along fuselage underside. 
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Figure 47 - Variation of model pressure with height above ground. 
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37b and ?8a). The high model temperature appears to be caused by blockage 
f the entrainment flow brought about by the presence of the model. As 
discussed in Subsection 5.2.5, higher temperatures were measured in the 
upwash with the model in place than without it. 

The model temperature distributions in Figure 49 do not show a local 
temperature minimxir'. jvist forward of the nacelle centerline that was found using 
the fan jet nozzles (Figure 31). In addition, comparison of Figure 49 with 
Figures 31 and 28 shows that the model temperature coefficients appear to scale 
better with jet exit temperature using the open circular nozzles than using the 
fan-jet nozzles, where only the core flow was heated. Temperature data taken 
in the upwash using the open circular nozzles (Figure 37b and 38a) also 
showed better temperature scaling than the fan-jet data (Figure 24). 

5.3 Conclusions 

5.3. 1 Fan- Jet Simulation. - The primary conclusion of this experimental study 
is that a geometrically correct small-scale model testing program can predict 
V/STOL aircraft hover performance qmte well. Interference force results of 
the 1/ 24-scale experiments matched full scale results closely. Pressure ratio 
changes (in this case with subsonic jet exit conditions) produced a neglibible 
effect. Core nozzle exit temperature produced a small effect, easily compensated 
for b3’ empirical scaling. Aircraft lower surface temperatures did not scale well 
with core nozzle exit temperature. 

5.3.2 Open Circular Nozzles. - The normalized interference forces resulting 
from operation with open circular nozzles were much larger (approximately a 
factor of two) than those found with the fan jets. The qu.'ilitative behavior 
with height above ground was the same. As with the fan-jet nozzles, there was 
a small change in the force with changes in stagnation temperature, " here was 
also a small change in nonnalized force with total pressure in the region of 
highest forces (H/D ' li), which did not occur for the fan-jet case. Aircraft 
lower surface temperature did scale well with nozzle exit temperature, a 
different result than that obtained with thf' fan-jet nozzles. 

The presence of the model weis found to raise significantly the temper- 
ature in the upwash. Unexpectedly, model presence also raised the total 
pressure in the upwash. 
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6 - PREDICTION METH .iLOGY 


The establishment of a vertical jet impingement model is the basis for the 
prediction of the behavior of two jets impinging on .ound plane. A sketch of 
the flow problem is shown in Figure 50. Figure 50a -itiows the vertical plane 
containing the jet stagnation points and the ground plane. Each jet impinges on 
the ground plane and deflects to form radial-wall jets. The wall jets then in- 
teract and form an upwash deflection zone where ihe wall jets collide and are 
turned upward on leavinj the ground plane. The stagnation line lies on the 
ground in the vertical plane of symmetry between the two impinging jets. The 
maximum upwash stagnation pressure occurs at the midpoint of the Kne connect- 
ing the two jei stagnation points on the ground. The ground pressure then 
drops off with lateral distance along Ihe stagnation line. Figure 50’o shows the 
radial streamline pattern that has been observed (Ref. 5, 6 & 7) both in the 
ground plane and in the vertical upwash plane of symmetry between the two 
jets. This flow situation only exists wncp the jets are spaced far enough 
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Figure 50 - Two-jet interacting flows - negligible 
deflection zone interaction. 
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apart so that the jet impingement zont ^oes not interact or has a negligible 
affect on the upwash defle'ition zone. (Ref. 3). 

e. 1 Theoretical Models for Ncn-Isothermal Jet Impingement 

In Ref. 3, semi-empirical analytical models were formulated to simulate the 
global behavior of two incompressible jets impinging on a ground plane. Models 
were generated for the behavior of the free and wall jets and, finally, the up- 
wash sheet generated by the interaction of the opposing wall jets. The effect 
of jet proximity is considered in Ref. 3, leading to the development of upwash 
momentum models. 

These incompressible models are extended , in the present study , to include 
temperature or density effects. The flowfield is divided into three major regions 
where viscous or turbulent mixing effects dominate, being; 

• Free jet 

• W. jet 

• Upwash. 

Subdomains of these regions include the jet and upwash deflection zones, in 
these smaller regions, the flow changes from an inviscid behavior to a shear flow 
along with streamline deflections, causing a change in static pressure due to the 
stagnation and acceleration of the flow. 

6.2 General Temperature a:; d Velocity Equations 

The transport or diffusion of heat in a free jet is much like the transport 
of momentum. The viscous mixing and entrainment of ambient air causes a shear 
layer of momentum (velocity) and heat to occur. SchUcting (Ref. 8), in quoting 
the experimental and theoretical observations of Reichardt (Ref. S), implies that 
the temperature profile distribution of a two-dimensional jet behaves like the 
square root of the velocity distribution Tay .or's Free Turbulence Theons as 
quoted by Abramovich (Ref. 10) , arrives at this result for the thermal layer by 
a turbulence theory that is based on vorticity transfer rather than momentum. 

In Ref. 10, the temperature profile in a jet is obtained by a "Mew Prandtl- 
Gortler Theory of Turbulence”. The following relation ’aip is obtained between 
the te ■ ■--ature and velocity : 


77 




where 5 is some constant. If 5 =1, the dimensionless temperature and velocity 
profile coincide. It has been observed experimentally that this behavior does 
not hold. For 6 - 2 , Taylor's analytical result is obtained as well as Reichardt's 
experimental observations. 


Experimental observations indicate that the thermal layer spreads faster 
than the dynamic (velocity) layer. In order to establish an empirical value for 
5 in eq. (1), some non -isothermal profile data are plotted in Figure 51a for a 
circular jet. Similarity is not obtained if the temperature data are related to the 
velocity half widths as shown in Figure 51b. In addition, the data do not seem 
to confirm the square root of the velocity relationship. If, on the other hand. 
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a) Related ;o jet radius. 
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Figure 51 - Dimensionless velocity and temperature 
difference profile data. 


78 



t r 





b) Related to jet radius nor,-dimen$ionaiized by velocity half-width. 
R81-1622 040(2/3)0 


Figrire 51 - (Continued). 

the velocity and temperature difference data are plotted dimensionless to their 
respective ha^f- velocity and half- temperature difference lengths, the data fall 
on one dimensionless curve as illustrated in Figure 51c. This result was also 
observed in Ref. 11 and in Ref. 12 for wall jets. It seems to be clear from the 
available data that the thermal layer persists beyond the point where a measur- 
able velocity exists. The form of eq. (1) implies that the dynamic and thermal 
layers have the same scale length. If an exponential or infinite layer is 
assumed for the velocity profile and, hence, the thermal layer, then eq. (1) 
would yield results consistent with measured data. Indeed, ’f such a relation- 
ship as eq . ( 1) exists , the exponent may be a function of the nozzle thermal 
conditions. 

In the models to be developed , it is more convenient to assume a finite 
thickness to the various layers to facilitate the various integrations involved 
in computing the velocity and temperature decay rates. Equation (1) is not 
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c) to jet radius non-di'r>ensicnaJizet/ by veiocity and temperature haif*widths. 

R81-1622*O40(3/3)O 


Figure 51 - (Concluded). 


used in the 
instead: 


following models but the relationship depicted by Figure 51c is used 


T-T 



( 2 ) 


where T = r~ , 
*>T 



and f = 



The dimensionless forms (i.e., velocity and temperature each have tneir 
own scale lengths, and b^) of the velocity and temperature difference pro- 
files as depicted ’.n eq. (2) are assumed identical. The relatonship between 
the thermal and dynamic (velocity) layers can be simply expressed by the ratio 
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of the half widths of the two layers: 


X = >1 l'*'' 

Hence, the thermal and dynamic layers do not have the same scale. The velocity 
goes to zero while the temperature persists beyond a measurable velocity. 


6.3 Conservation Equations 


The basic governing equation to model the behavior of jots is the conser- 
vation of momentum equati« n. To account for non-isotherma! or temperature 
e^'fects, an additional equation is required. The conservation of excess heat 
content is used as the governing equation to account for heat diffusion or 
temperature effects. These two integral equations can be expressed as: 


Momentum - 


Heat Flux 


M, 


H, 





J" J 

N =//’oVATdA 


(4) 

(5) 


where AT„ = T„-T . , T„ is the nozzle temperature, and A„ is the nozzle area. 


Eouations ( 4 ) and (5) yield two expressions for the two unknowns, maximum 
centerline temperature and velocity. Unfortunately, more than two unknowns 
exist, namely, the profile shapes and growth rates of the various shear layers. 
These additional parameters must be given by empirical observations. The 
spreading models of Ref. 3 are used with the addition of eq. (5) and the further 
compUcation of the density or temperature occurring in eq. (4) which necessi- 
tates the simxUtaneous solution of eqs. (4) and (5). 

6. 4 Heated Free-Jct Model 

The nozzle exit conditions of a free jet can be specified by two parameters, 
the nozzle exit stagnation pressure ratio (NPR) and temperature ratio (NTR). 

The exit Mach number of the jet is purely a function of NPR. .4, compressibility 
correction is included in the incompressible models, based up.n an approximate 
Mach nurii, er computed using Bernoulli’s equation. 
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M 

To be consistent with an incompressible approximation, the nozzle exit static 
temperature is assumed equal to the nozzle stagnation temperature, 
or 

T„ . NTR 
N 

The free jet is subdivided into three regions, as shown in Figfure 52; 

• Potential core , Region I 

• Transition, Region II 

• Fully developed. Region III. 

In order to solve eqs. (4) and (5), some information about the half-velociiy 
boundary growth behavior must be assur ed to be known from experimental data. 
To be consistent with the models of Ref. 3, and for the sake of simplicity, the 
jet radius mcdel assumes a simple linear boundary growth in each region. 
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Figure 52 - Jet half-velocity wi:nh model. 
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T^'_ cc stants (aj, aj. as) governing the half-velocity boundary growth 
in the tru; \~i. and fully- developed regions are determined empirically. The 
poteuti'i rrjre length (Zp^) and the exponent oipj^ governing the shape of the 

shear profile in the fully developed region also must be specified. These con- 
stants vary slightly as functions of nozzle Mach number and temperature ratio. 
Suitable values have been determined as: 



/ 2 \] 

1 


as = .09 

1 + ln| 1 + 1 V 

V 

* 

1 

1 - 


where S3-.C.50. 

The potential core length is approximated by: 



where 3^~0.75, and is also slightly a function of the noz-^’e temperature ratio. 


The growth constant of the half-velocitjf boundary in the potential core 
region is determined by matching the boundaries at the end of the potential core, 
yielding. 


( b.-l ) 




In addition , r. linear decay in potential core radius , is assumed : 



The beginning of the fully developed region of flow, denoted by the length 
parameter Zp^, is determined by matchinp" the boundaries between the transition 

and fully developed regions. 



"“n (as-a,) 
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The momentum and heat flux equations become, for a circular jet. 


M 


2~p i-V^rdr 


J = 


H 


F =- 


2ttP 

i 


J v ^r 
— 

/ 


VAT rdr 


using = oRqT to eliminate density. 


( 6 ) 

(7) 


The velocity and temperature profiles are represented by the following 
dimensionless self-similar qufiuitities: 


V 






AT 

AT 


M ^ 


(n,j.) 


( 8 ) 

(9) 


Since finite profile functions will be used in the modeling, the velocity and 
temperature scale lengths, as represented by b^ and b,p, are not equal. 


Substituting eqs. (8) and (9) into eqs. (6) and (7), and rearramging. 


yields , 


Momentum 


N \''n/ \ N,' / |-f^ 


Heat 
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M 


•N 


= 2 


aYAvVAIm / VtAvAv 

L J 


( 10 ) 


(11) 


where T„ 



Unfortunately, bringring the density into the integrands makes the profile 
integrals functions of the maximum centerline temperature ratio. Hence, simple 
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geometrical similarity does not exist for non-isothermal jets. The velocity scale 
length is used in the integration since there can be no contribution to eq. 

(10) or eq. (11) beyond the dynamic (velocity) layer. The integrals in the 
above equations do not lend themselves necessarily to analytical integration , 
and are defined as: 


C (T ) - /* 

Vo |~f.. + 


M 




/■ Vl 
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fTT^-pHydny 

(1-f.J )1 
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Equations (10) and (11) can be rewritten simply as 
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= 21 
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(12) 

(13) 

(14) 

(15) 


The values of the integrals and C.p are also functions of the profile 
characteristics. In general, it is assumed that 


— J=V'V = 

M 



(16) 
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tj(n.j.) 


(17) 


where 
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C/by- 


8S 



Since the dimensionless velocity and temperature profile functions have the 
same form, then 


n 


T = 



(18) 


In general, it is assumed that 




Hence : 




(19) 

( 20 ) 


and the temperature scale length can be directly related to the scale length 
of the dynamic layer. 

6.4.1 Potential core region - In the potential core region, the maximum 
temperature and velocity are known and are equal to the nozzle values. The 
momentum and heat equations, (14) and (15), can be reduced to. 


-I 


XH 


Cm (a, X,j.) = 0 



V 


^(a) 



(ct,X^) - 0 


where Ry 



(1-n ) 
c 


1 



( 21 ) 


( 22 ) 


Equations (21) and (22) are solved simultaneously for the exponent a and the 
thermal-to-dynamic layer ratio For simplicity, the thermal and dynamic 

cores are assumed identical. 

6.4.2 Transition and fully-developed regions - The conservation equations can 
now be written as: 
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The above equations are solved simultaneously for the dimensionless max- 

T V 

M M 

imum temperature ratio, , and maximum velocitv ratio, , given values for 

the exponent a and Suitable values based on empirical data are assigned 
to these parameters in the fully developed regions as 

X,j, = 1 + 0.185 where B^^O.50, 

and a linear variation from potential core values to the fully developed ones is 
assumed. 


Some typical results are shown in Figures 53 through 57 and compared to 
the test data of Corrsin & Uberoi (Ref. 13) for low Mach number or basically 
incompressible jets. Figure 53 shows the relationship between the dimensionless 
velocity' and temperature decay rates for a slightly-heated jet (NTR = 1.05) and 
a significantly-heated jet (NTR = 2.0). Both the velocity and temperature decay 
rates are enhanced due to an increase in heat content of the jet. The tempera- 
ture decays faster than the velocity. Figure 54 shows the dynamic pressure 
or tots’ pressure decay for the two different temperature ratios. The total 
pressure decays faster for the hotter jet. Figure 55 shows the relationship 
between the dynamic pressure and total temperature decay. Unlike the velocity 


and temperature, the pressure and temperature curves cross each other. The 


dynamic pressure is computed as Figure 56a and b shows tne 

relationship between the half dynamic pressure (q), half velocity (V../V ) and 

[VI N 


halt dimensionless temperature (T^^-T^ /Tj^-T^) boundaries for the two nozzle 
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(Corrsin & Uberoi, Hef. 13) 
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Figure 55 - Comparison of computed free-jet dimensionless dynamic 

pressure and temperature difference decay with test data 


temperature ratios. The half dynamic pressure boundary exhibits a contraction 
in the potential core region (Zp^ ~ 4D). The half velocity boundary lies between 

the half dynamic pressure and temperature boundaries. The half temperature 
boundary being significantly wider than the half velocity boundary as the 
nozzle temperature is increased. Figure 57 shows a comparison of the half q 
and half temperature boundaries for the two nozzle temperature ratios. Hence, 
heating the jets causes an increase in the jet half boundaries, which is con- 
sistent with the higher decay rate exhibited by the hotter jet in the previous 
figures. 

6.5 Jet Deflection Region 

A schematic of the wall jet impingement model is shown in Figure 58. The 
jet begins to stagnate ai some ground effect height given by Aj. Values for 

Aj are given in Ref. 14. The ground stagnation pressure and temperature is 

then determined by the jet properties at the ground 'ffect height denoted by 
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Figure 57 - Heated free- jet boundary growth characteristics. 


the subscript (g), or: 



V /V-, 
g N 



T /T 




T 


A 


(25) 


The ground pressure distribution in the impingement region is obtained in 
the same fashion as for an isothermal jet. A momentum balance is performed with 
the integral of the ground pressure equal to the momentum or thrust of the 
incident jet. 

r 1/2 

M ° (26) 

—r=S / i Prdrd 

4 ^ 

00 


The ground pressure distribution is assured to have the form: 


where 




f f ^ ^ 

= ig ^ = 




(27) 
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Equation (26) can be rewritten as, 


1 = 




sp 


'V 



( 28 ) 


where dr,^ 

0 ‘ 


Given the impingement radius r^, equation (28) can be solved for the 
exponent ot governing the shape of the ground pressure distribution. The 
deflection radius r^ is estimated empirically and is based upon the incident 


half velocity radius (b,.,,:) of the jet: 

VH 


2.9 < 3.6 ; 1.2 <_< 3.0 

^VHi 


h 
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r h „ _ 

o _ , fi . — > 3.0 

Shi'-® • “ 


6.6 Heated Wall- Jet Transition Model 

Figure 58 shows a schematic of the wall- jet model. The wall- jet transition 
model developed in Ref. 3 consists of three subregions: 

i) Inviscid deflection region where the effects of viscosity are assumed 
to be negligible except in a region close to the wall and near the edge 
of the deflected flow. In this region, the inner boundary layer is 
established as governed by axially-symmetric stagnation flow. 

ii) Transition region where the effects of viscosity are beginning to 

dominate and the inner boundary layer and outer shear flow changes to 
the fully-developed turbulent wall-jet. 

iii) Fully- developed flow, where the effects of turbulent viscosity dominate 
and the static pressure through the wall layer is considered ambient , 
or fully recovered. In this region, the nearly similar wall- jet develops. 

The addition of heat to the wall-jet requires the introduction of a radial 
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heat flux equation into the model of the isothermal wall-jet. The two con- 
servation equations (assuming negligible losses) applied to the wall jet are: 


Momentum - 


Heat - 


2 

Mp = J (pV + AP) rdz 


0 

Hp = f pVATrdz 


( 29 ) 


(30) 


where is the height or thickness of the dynamic (velocity) layer. 

Before being able to apply eqs. (29) and (30) to the wall layer, the 
velocity, pressure and temperature profile behavior must be approximated. 

Two sublayers are assumed, an inner boundary layer and an outer shear layer. 
Adiabatic wall conditions are assumed and the thermal and dynamic boundary 
layers are also assumed to coincide. 


Inner Boundary Layer 
1 



where 7 4 N ^ 15. 
Outer Shear Layer 


where ^ 


Z 

b* 
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K 



(31) 


(constant pressure through boundary layer) 


(adiabatic wall, constant temperature 
through boundary layer) 







(32) 
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z 

where =^— and is the height or thickness of the thermal layer. The 


dynamic and thermal layers are assur ad to be related by; 



( 33 ) 


The introduction of equations (30) and (31) into eqs. (28) and (29) lead 
to the development of several pTOfile integral expressions. 

The momentvun integral has a velocity squared and a pressure term. The 
pressure term contributes to the momentum until pressure recovery occurs. 
Substitution of the profile expressions into the velocity term of the integral 
yields the following integrals for the two layers : 



(34) 


where ~ ^ 

Ta 

The pressure term leads to a simpler expression : 


1 



^6 


which can be integrated analytically. 

The heat equation leads to integrals similar to eq. (34): 


(35) 
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(36) 
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T 




pT(^r) 


The pressure distribution along the wall was defined previously as: 


^^S^g (^w) 


(37) 


Equations (29) and (30) can now be expressed as 


Momentum - 


'm (^ w) 


VIS 


(r) =1 
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(38) 


N 




Heat Flux - 



(39) 


The terms on the left side of eqs. (38) and (39), and F,j, , reflect the 


gain in radial momentum or heat flux as the jet deflects and becomes parallel to 
the wall. These functions are made proportional to the pressure recovery 
function, or 


^m(^w) ^t(''n) (^w) 


(40) 


since little or no empirical data are available in the deflection region. Hence, 
the radial momentum and heat flux reach a maximum at pressure recovery, and 
for 


n 


W 



1 


The term F^j^lr) reflects the loss of momentum due to frictional loses at the 
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wall, and 


(41) 


(I""’' 

Equations (38) and (39) represent two simultaneous equations for the two 
unknowns, and The values of the integrals and C.^,, which arise 

from the profile functions, are determined numerically and are, in general, a 
function of the maximum temperature ratio (Tjyj ). 


The inviscid deflection region is assumed to occur under the half velocity 
width (byjjj) of the incident jet profile. Constant total pressure and total 

temperature is assumed throughout the inviscid layer to be equal to the incident 
jet values. The wall-jet properties at the beginning of the wall- jet transition 
region are obtained from the inviscid values. The initial velocity at the end 
of the inviscid region is defined from the static pressure and Bernoulli’s 
equation. 

Equations ( 38) and ( 39) are used to determine the initial values of the 
exponent a , governing the wall- jet profiles in the shear layer , and the initial 
ratio. \^ . of the thermal-to-dynamic layer thickness. The wall layer is now 

completely initialized upon prescription of the initial boundary layer 
characteristics. Fully developed values for the exponent a and the ratio \,j, 

are now prescribed and a linear variation between the initial values at the be- 
ginning of the transition to the fully developed region is assumed. Equations 
(38) and (39) can then be solved simultaneously for the maximum temperature 
and velocity throughout the wall- jet. given a prescribed behavior for the half- 
velocity thickness. 


For the fully developed wall-jet region, the boundary and shear layers are 
assumed to have a simple linear behavior: 





(42) 




(43) 
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v/here 






= 0.0175 

3 

and M is the stagnation Mach number of the jet. The wall jet thickness in the 
s 

transition region is obtained by a linear variation between the initial and fully 
developed values. Some total pressure and temperature decay data is shown 
in Figure 59 for single circular jet impingement into a radial wall- jet. The 
wall-jet pressure decay is only slightly influenced by the nozzle temperature 
ratio for the range tested. The temperature decay shows a greater dependence 
on the nozzle temperature. Both the temperature and pressure decays more 
rapidly with an increase in nozzle temperature. 


6.7 Two- Jet Interaction Model 

Figure 60 shows a sketch of the ground plane coordinate systems for the 
two jet impingement interaction problem. The jets stagnate on the ground, de- 
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Figure 59 - Heated wall-jet dimensionless total pressure and 
temperature difference dacay with ground radius. 
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Figure 60 - Two- jet impingement without deflection zone interaction. 


fleet into radial wall jets and collide along their plane of symmetry. The 
formation of an upwash stagnation zone is a result of the upward deflection of 
the wall jet flow . The point , on the line joining the jet centers , is the only 

true stagnation point where the flow comes to rest and deflects vertically. At 
points off of the center, the flow is deflected at an angle such that the pressure 
is representative of the stagnation of the normal component of the maximum wall 
jet velocity. 

6.7.1 Maximum pressure and temperature distribution along the upwash stagna- 
tion line - To estimate the upwash stagnation line ground pressure and 
temperature distribution, the maximum velocity of the wail jet normal to the 
stagnation line is computed: 
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The maximum pressure along the stagnation line is computed from the static 
phis the contribution of the normal velocity of the wall jet. 

Pressure - 

* MW 

(45) 


1 P 

AP = AP +"“-M55L V ^ cos^ 

^ umg 2R„T.„., MW u 


G MW 


Dividing by the nozzle pressure given by 
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Temperature 
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umg _ MW 
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(48) 


The maximum temperature along the stagnation line is then just equal to 

3 

the maximum temperature of the wall- jet at the stagnation line location, X = -^ 

The pressure along the stagnation line given by eq. (47) uses the maximum 
temperature, velocity, and static pressure of the wall-jet evaluated at the 
stagnation line location. 

6.7.2 Upwash momentum models - The upwash momentum model, including the 
effects of close jet spacing, remains unchanged in principle due to temperature 
effects and are those described in Ref. 3. The addition of temperature effects 
has a slight effect on the ground pressure distribution since the addition of 
heat to the flow somewhat affects the overall decay rates of the various regions. 
The overall upwash momentum can then be related to the integral of the upwash 
ground pressure distribution as described in detail in Ref. 3. 


ICO 



6.7.3 Heated upwash decay model - Figure 61 shows a schematic of the upwash 
sheet model and the characteristic scaling parameters. The upwash sheet is 
assumed to be a reflection of the wall-jet flow into the vertical plane of symmetry 
lying between the jets. The radial streamline pattern of the wall layer is 
assumed to continue into the upwash sheet and, to a first approximation, be un- 
perturbed by the turning region. The wall-jet flow is assumed to enter the up- 
wash deflection region with a characteristic profile, half velocity width and 
maximum velocity or total pressure. The characteristic length scale in the up- 
wash sheet is the half velocity width of the incident wall-jet layer estimated at 
the wall location. The pressure recovery region in the upwash sheet is 
assumed to be approximately three times the half velocity width of the incident 
wall jet profile. The magnitude of the momentum flux in the resulting upwash 
streamline is assumed equal to that of the incident wall-jet. Hence, the upwash 
sheet is treated in a similar fashion as the wall-jet with a few exceptions. 

Due to the high turbulence levels typically measured in the upwash sheet , 
the flow is considered to be fully turbulent and similar. The upwash velocity 
profile is taken to be that of a free shear profile with no internal momentum 
defect due to the waD layer. The upwash sheet growth rate is assumed to be 
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Figure 61 - Characteristic scaling parameters for upwash model. 
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constant without any initial transition phase. It is also assumed that the 
magnitude of the momentum flux distribution in the upwash sheet is pven by 
the previously discussed momentum models. Expressions for the momentum and 
heat flux similar to those for the wall jet can be derived for the upwash flow : 



( 49 ) 


where 
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The pressure terms in eqs. (49) and (50) only have a contribution in the 
upwash deflection regions. The above integrals, as with the wall-jet and jet 
flows, are carried out over the w'dth of the dynamic layer given by b^. The 

thermal layer is considered to be somewhat wider and is given by: 


buT = 


The growth rate of the upwash flow is assumed to be approximately three 
times the rate of the wall-jet due to the higher turbulence. Equations (49) and 
(50) are then solved simultaneously for the temperature and velocity. Initial 
application of this model to the heated upwash led to the results shown in Figure 
62 for the pressure and temperature decay. The pressure correlates well with 
the test data but the temperature decays too rapidly. The wall- jet and upwash 
temperature behavior for two interacting impinging jets is significantly different 
from a single isolated jet impinging to form a wall- jet. Recirculation effects 
dominate for multi-jet impingement. These effects are due to several aspects of 
the two- jet problem. Hot upv/ash flow tends to recirculate back into the free jet 
and wall jet regions, altering their effective temperature decay. Confinement 
of the region between the free jet and upwash flow tends to heat the ambient 
flow or restrict the influx of cooler ambient air into this region to be entrained 
by the wall jet and vipwash flow. As a result of these effects, in order to be 
able to estimate the temperature behavior of wall-jets and upwash flow , a 
recirculation model must be established. 


6.7.4 Recirculation model - As a preliminary attempt at estimating recirculation 
effects, a model was generated for both the wall-jet and upwash flows. The 
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Figure 62 - Heated upwash centerline decay characteristics with 
recirculation effects. 


basic principle of such a model is demonstrated in Figure 63. The model takes 
into account that the wail jet and upwash streamlines do not effectively er. train 
air at the reference ambient condition. This situation is most extreme along 
the inner wall- jet and upwash streamline (<^=0°) and is enhanced as the jet 
spacing is reduced. The recirculation effect is also assumed to be negligible 
along the outermost ray In reality, the outermost ray may also be 

affected if the free jet decay has been altered. This is neglected in the present 
study. In order to quantify the recirculation effect, it is assumed that it 
results in a local change in ambient conditions represented by: 


"Aeff , ^ 2 

:= 1 + cos 
1 . rec 

A 


(-f) 


AT 


(54> 


where 


A 


rec 


g 


(^’ 2 ) 


rec 


:i ^ 1.30 
rec 


and AT is the ground stagnation temperature of the incident jet. 

g 
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Figure 63 - Recirculation effects. 


Essentially, the amount of excess heat flux in any given layer is then 
altered by the effective local ambient condition: 



^Aeff 


Hence . if the effective or local ambient temperature is equal to the nozzle 
temperature, the excess heat flux will be zero ar.d no temperature decay will 
take place as a result of the heat equation. 


If the incident radius (r ) of the impinging jet is equal to the half spacing 
of the jets, then the effective ambient temperature becomes equal to the 



incident temperature ratio of the jet along the stream e f> =0" and no further 
temperature decay will occur. Figures 64a and 64b demonstrate the effect of 
this model on the wall jet decay rates for two jets with S/D = 3.0 and H/D = 3.0, 
and two different nozzle temperature ratios. In this figure, the decay is 
computed until the wall jet interacts vith the stagnation line. For the 

wall-jet decay is computed to a fixed radius of three diameters. The least 
temperature decay occurs along the inner (i=0°) streamline. The wall-jet 
decays more rapidly as the outermost ray is approached and. at i = -. the decay 
of the wall- jet becomes equal to the isolated impingement value. 
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ORIGIN, 1 P V - r3 
OF POOR QOALnY 

A slight overshoot of the values occur at the end of the inviscid region 
due tc a n’!iniatch either in the radius of the inviscid region or in the thickness 
of the wall layer. Even though an effective change in ambient temperature 
conditions affects the temperature decay rate dramatically, it only affects the 
total pressure decay slightly as indicated by the figures. Figure 65 more clearly 
demonstrates the effect of the recirculation model as a function of jet spacing, 
the decay of the innermost streamline (>>=0) is plotted as a function of jet spacing. 
.As the spacing increases, the wall- jet temperature decays more rapidly and, as 

. the decay rate of the i.aeracting wall jet approaches the isolated or single 
jet impingement behavior. The recirculation effect on the wail- jet feeds into the 
upwash as an alteration of the initial upwash temperature as well as changing the 
up wash decay rate. Figure 66 shows the improved correlation in ihe upwash 
temperature decay with the implementation of the recirculation model in comparison 
with Figure 62. Obviously, more detailed correlation is required xeally to be able 
to quantify this effect . It is also expected that the presence of a confining plate 
or vehicle will significantly alter or enhance the recirculation effects. 



Tomp^riturc 

QO 


Pr«4sur« 


Rai-l622-056(l/2)0 

Figure 65 - Heated-wall jet decay characteristics ( y? = 0 °) vs. jet spacing. 
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7 - CONCLUSIONS 


For the configuration considered in this investigation , small-scale models 
can be used to predict full-scale induced forces. No scale effects were found, 
and jet temperature and pressure did not affect the nondimensionalized induced 
lift. The induced lift in- ground-effect was found to be fiigher for the uniform 
circular jets than for the simulated fan jets. 

■Modifications to an existing wall-jet transition model adequately predict 
the trends with height above ground of upwash temperatures and pressures. 

The addition of a recirculation model is necessary to predict upwash temper- 
atures in the presence of an aircraft. 
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APPENDIX 


COMPUTER PROGRAM DESCRIPTION 


PROGRAM NAME : GRUMHOT2 

PURPOSE: VTOL TWO -JET IMPINGEMENT INTERACTION PROBLEM 

This program is designed to estimate the dynamic and thermal flow 
characteristics associated with two vertically-impinging and equal-strength hot 
jets. This program is specifically oriented towards the problem of closely- 
spaced jet interaction, where the deflection regions interact until eventually jet 
coalescence occurs. 

Aside from the basic flow characteristics, the program assumes the 
symmetrical placement of a slender fuselage in the upwash flow . The upwash 
lift force is then computed for a cylindrical fuselage of constant cross-sectional 
shape. The body parameters, in terms of width and depth, do not vary 
longitudinally. The program estimates the force on the basis of a rectangular 
and a circular cross section. Two planes of symmetry are assumed and all out- 
put appUes to one quadrant of the flow field (ie. equal jets and nozzles located 
at the midpoint of the fuselage) . Jet entrainment effects may be significant but 
are neglected in this program. Hence, the force is only that due to upwash 
impingement. Some residual programming exists in the code for a parabolic 
body of revolution. These cards have been commented out but may be used if 
desired. 


INPUT DESCRIPTION 

Note: All input parameters are in terms of nozzle diameters. 


Card No . 
1 


Code Names 
NPR, NTR 


Name Defin ition 

NPR Nozzle pressure ratio 

NTR Nozzle temperature ratio 


Format 

2F10.5 


110 



Card No. 


Code Names 


Format 


2 


HD, SD, ZPLD, DZPL, ZFINAL 


5F10.5 


Name D efinition 

HD Nozzle height above ground 

SD Nozzle spacing 

Note : The program will compute one or several positions of the fuselage 
relative to the ground for a fixed nozzle height above ground. 

ZPLD Initial fuselage height above ground 

DZPL Increment in fuselage AZ above ground 

ZFINAL Final Z coordinate of fuselage height relative to ground 

plane. 


Card No. 
3 


Code Names 
XL2, WCON, ZCON 


Names 

XL2 

WCON 

ZCON 


Definition 
Fuselage Length 
Fuselage Width 
Fuselage Depth 


Form at 

3F10.5 


Note: ZCON determines position of fuselage underside relative to its ZPLD 
location. Bottom of fuselage wiU be located at ZPLD- ZCON at first 
computed location. The upwash sheet properties are also computed at 
this Z location. 


Card No . 
4 


Code Names 
IPBAR 


Format 

II 


Name Definition 

IPBAR Integer controlling the output of ground pressure 
pattern: IPBAR = 0, no pressure pattern output 
IPBAR = I, pressure pattern output is desired 

Note: Card No. 4 is not required if IPBAR = 0. 
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Card No 


Code Name 


Format 


5 NU 

Name Definition 

NU Number of pressure values to be input for computa- 
tion of ground isobar pattern NU < 25 

Note: Card No. 5 is repeated NU times. 

Card No. Code Name 

6 PU 


Name Definition 

PU Nondiraensional pressure for isobar pattern 

0.0 < PU < 1.0 

^igure 67 shows a typical input set. 
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Format 

F10.5 


rof: 


INPUT DATA 




t.007 2.00 

4.00 3.000 2.00 1.00 4.00 

8.00 l.OO 0.0 

1 

10 
.9 
.8 
,7 
.60 
.50 
.4 
.3 
^ 2 
.10 
.05 


EOF; 


R8I'1622<062D 


Figure 67 Typical program input 



PRINTED OUTPUT DESCRIPTION 


Figure 68 shows a typical computer printout. Most of the geometrical 
output quantities are nondimensionalized by the nozzle radius. All velocities 
and pressure are initially nondimensionalized by the nozzle exit velocity and 
stagnation pressure. All pressures are relative to ambient conditions. 

Note: RN and D refer to nozzle exit radius and diameter. VN refers to nozzle 
exit velocity. 


Input Parameters 

The first set of output echoes the input parameters. 


Output Titles 
NPR 
NTR 

Mach No. 

TN 

H/D 

S/D 

Z/D 

DZ/D 

ZFINAL/D 

L/D 

W/D 

ZB/D 


Definitions 

Nozzle pressure ratio 

Nozzle temperature ratio 

Mach number computed from NPR and NTR 

Nozzle temperature ratio (NTR) relative to ambient 

or TN/TA 

Nozzle height above ground 

Nozzle spacing, distance between jet centerlines 

Initial fuselage height 

Increment in fuselage position 

Final location of fuselage 

Body Length 

Body Width 

Location of underside of body relative to specified 
fuselage location 


Note : If ZB « 0, underside location is coincident with specified fuselage 
position . 


Jet Decay Region 

Output Titles 
DELG/D 
ZPC /RN 
Z7D/RN 
Z/RN 


Definitions 

Jet ground effect height relative to ground plane 
Length of potential core 

Length of potential core and transition regions 
Jet axial location measured from nozzle exit 
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ORiG:.- ^ 
CF PCCK 


««• INPUT PAPAHCTCRS M« 


HPR- 1.0070 NTR- 2.0000 

HACH N0.« 0.1000 TN/TA« 2.0000 

H/0- 4.00000 S/D« 3.00000 Z/P* 2.00000 

Z/0- 2.00000 PZ/D» 1.00000 ZFINAI /D* 4.00000 


iODY LENGTH L/D* 0.00000 HIPTH U/0- 1.00000 DEPTH ZO/P' 0.0 


t iittf JET DECAY REGION ttttt 


JET DEFLECTION HEIGHT DELO/D • l.lf340 

ZPC/RN- 6.016t3 ZFD/RN* 12.78302 


Z/RN 

RJH/RN 

RJ/RN 

ALPV 

VJ/Vtt 

RC/RN 

ALPT 

RJHT/RN 

TH/TH 

TH/TA 

DTN/DTN 

RJHG/RN 

0.0 

1.0000 

1.0000 

0.0 

1.0000 

1.0000 

1.3000 

I .0000 

1.0000 

2.0000 

1.0000 

1.0000 

1.7016 

1 .0181 

1.0673 

6.0144 

1.0000 

0.8009 

4.0144 

1 .0181 

1.0000 

2.0000 

1.0000 

1.0047 

3.4033 

1.0363 

1.1363 

6.3209 

i.OOOO 

0.3691 

6.3209 

1.0363 

1.0000 

2.0000 

1.0000 

1.0088 

3.1049 

1.0344 

1.2083 

6.5791 

1.0000 

0.3033 

4.3791 

1.0344 

1.0000 

2.0000 

l.OOOC 

1.0120 

6.8066 

1.0726 

1.2830 

6.7837 

1.0000 

0.0019 

6.7837 

1.0726 

1.0000 

2.0000 

1.0000 

1.0138 


4tt$i JET DEFLECTION REGION 9tttt 

ROH/RN- 1.07237 RG/RN- 1.28301 RO/RN- 3.8A124 UG/VN* 1.00000 
STAGNATION PRESSURE* DPS/DPTJ- 1.00000 GROUND NAX. TEHP.t T08/TA« 2.00000 

STAGNATION HACH NUH8CR OF JET- 0 .10000 

AlVg- 1.50462 

SINGLE JET GROUND PRESSURES 


R/ROH 

DP8/DPTQ 

R/RN 

DPS/PPTJ 

0.0 

1.00000 

0.0 

1.00000 

0.13000 

0.96689 

0.16088 

0.96609 

0.30000 

0.90821 

0.32177 

0.90821 

0.43000 

0.83608 

0.48263 

0.83608 

0.60000 

0.73619 

0.64334 

0.73619 

0.73000 

0.67237 

0.80442 

0.67257 

0.9000C 

0.38832 

0.96331 

0.58832 

1 .OSCOO 

0.30392 

1.12619 

0.30592 

1.20000 

0.42734 

1.28708 

0.42734 

1.33000 

0.33409 

1.44796 

0.33409 

1.30000 

0.28731 

1.60883 

0.28731 

1.63000 

0.22773 

1 .76973 

0.22773 

1.80000 

0.I/3S6 

1.93062 

0.17386 

1 .93000 

0.13173 

2.09130 

0.13173 

2.10000 

0*09328 

2.23239 

0.09328 

2.23000 

0.06606 

2.41327 

0.06606 

2.40000 

0.04330 

2.37413 

0.04330 

2.35000 

0.02683 

2.73304 

0.02683 

7.70000 

0.01324 

2.89392 

0.01324 

2.85000 

0.00771 

3.05681 

0,00771 

3.000^0 

0.00331 

3.21769 

0.00331 

3.14999 

0.00110 

3.37838 

0.00110 

3.29999 

0.00023 

3.33946 

0.00023 

3.44999 

0.00001 

3.70033 

0.00001 

3.39999 

0.00000 

3.86133 

0.00000 


wall jet region «t4il« 


STAGNATION POINT DOUNDAOY LAYER THICKNESS* DELS/RN- 0.03733 
PHI/FA-- 1.00373 (PTTOT/PA-1 . ) • 0.00700 (T0T/TA-1.>- 1 . Ov/000 

THI/T.4- 2.00000 VNI/UH- 0.68210 

VELOCITY AT START OF TURiULENT UJ^LL JET* VH/VO- 0.68210 

START OF WALL JET REGION* ALPUO- 0.60408 DHOH/RN- 0.34307 8U0/RN- 2.38600 

R81*1622*063O(19 Sheets} 

Figure 68 Typical program printout. 
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GJ/GN 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 



ORIGINAl. F.^ 

OF POOR QUALiT 

ISOLATED WALL JET FROFERriES 


R/RN 

UH/VO 

VH/VN 

iUH/RN 

8U/RN 

DEL DL/RN 

DELFJ 

KDEL 

Th/TH 

TH/TA 8UHT/RN DTH/DTN 

0.0 

0.0 

0.0 i 

0.3431 

2.3868 

0.0374 

1 .0000 

0.0136 

1 .0000 

2.0000 

0.3383 

1 *0000 

0.2300 

0.2312 

0.2312 

0.3431 

2.3868 

0.0374 

1 .0000 

0.0136 

1.0000 

2.0000 

0.3383 

1 .0000 

0.3000 

0.4!30 

0.4138 

0.3431 

2.3868 

0.0374 

1.0000 

0.0136 

1.0000 

2.0000 

0.1383 

1.0000 

0.7300 

0.3433 

0.3455 

0.3431 

2.3868 

0.0374 

1.0000 

0.0136 

1.0000 

2.0000 

0.3383 

1 .0000 

1 .0000 

0.6342 

0.6342 

0.3431 

2.3868 

0.0374 

1.0000 

0.0136 

1.0000 

2.0000 

0.3383 

1 .0000 

1 .2SC0 

0.7238 

0.7238 

0.3483 

2.0313 

0.0370 

0.9818 

0.0182 

0.9844 

1.9687 

0.3479 

0.^687 

1.3000 

0.7622 

0.7622 

0.3329 

1.6938 

0.0366 

0.9376 

0.0216 

0.9610 

1 .9219 

0.3613 

0.9219 

1.7300 

0.7764 

0.7764 

0.3373 

1 .4680 

0.0362 

0.8799 

0.0247 

0.9357 

1.8713 

0.3730 

0.8713 

2.000'» 

0.7741 

0.7741 

0.3620 

1.3087 

0.0338 

0.8153 

0.0273 

0.9096 

1.8193 

0.3889 

0.019? 

2.2300 

0.7377 

0.7397 

0.3666 

1.1916 

0.0*133 

0.7492 

0.0297 

0.6836 

1.7673 

0.4031 

0.7673 

2.5000 

0.7367 

0.7367 

0.3712 

1.1028 

0.0349 

0.6833 

0.0316 

0.6386 

1 .7173 

0.4174 

0.7173 

2.7300 

0.7081 

0.7081 

0.3738 

1.0338 

O.0343 

0.6239 

0.0333 

0.B3S2 

1.670* 

0.4320 

0*6704 

3.0000 

0.6764 

0.6764 

0.3803 

0.9789 

0.0340 

0.3722 

0.0348 

0.8138 

1.6276 

0.4469 

0.6276 

3.2300 

0.6439 

0.6439 

0.3849 

0.9347 

0.0336 

0.3243 

0.0360 

0.7943 

1.3889 

0.4619 

0.3809 

3.3000 

0.6123 

0.6123 

0.3893 

0.8983 

0.0332 

0.4827 

0.0369 

0.7773 

1.3343 

0.4772 

0.3343 

3.7300 

0.3830 

0.3830 

0*3941 

0.8683 

0.0328 

0.4461 

0.0377 

0.7620 

1.3239 

0.4927 

0.3239 

4 . 0000 

0.3476 

0.3476 

0.4103 

0.8860 

0.0330 

0.4023 

0.0393 

0.7434 

1.4907 

0.3177 

0.4907 

4.:soo 

0.3071 

0.3071 

0.4360 

0.9386 

0.0394 

0.3320 

0.0419 

0.7306 

1.4612 

0.5300 

0.4612 

4.3000 

0.4718 

0.4718 

0.4616 

0.9912 

0.0437 

0.3102 

0.0441 

0.7176 

1.4352 

0.5824 

0.4332 

4.7300 

0.4409 

0.4409 

0.4873 

1.0438 

0.0481 

0.2733 

0.0461 

0.7060 

1.4120 

0.6148 

0.4120 

3.0000 

0.4133 

0.4133 

0.3129 

1.0964 

0.0323 

0.2458 

0.0479 

0.6937 

1.3913 

0.6471 

0.3913 

5.2300 

0.3892 

0.3892 

0.3386 

1.1490 

0.0369 

0.2207 

0.0493 

0.6864 

1.3727 

0.6793 

0.3727 

3.3000 

0.3674 

0.3674 

0.3642 

1.2017 

0.0613 

0.1991 

0.0310 

0.6779 

1.3338 

0.7U8 

0.3338 

3.7300 

0.3477 

0.3477 

0.3899 

1.2343 

0.0636 

0.1804 

0.0323 

0.6703 

1.3405 

0.7442 

0.3403 

6.0000 

0.3300 

0.3300 

0.6133 

1.3069 

0.0700 

0*1642 

0.0316 

0.6633 

1.3263 

0.7>63 

0.1263 


•«$•$ TUO-JCT INTERACTION PRODLEN tittt 


ttt RECIRCULATION EFFECT ON NALL JETS tt$ 


PMI« 0.0 DEGS. 


TNEFF- 2.00000 TAEFF- I . 4283A 


R 

UN/UN 

TH/TN 

TH/TA 

DFH/DPN 

OTH/DTN 

0.0 

0.0 

1.00000 

2.00000 

1 .00000 

1.00000 

0.12300 

0.14993 

1.00000 

2.00000 

1.00000 

1.00000 

0.25000 

0.25116 

1.00000 

2.00000 

1.00000 

1.00000 

0.37300 

0.33728 

1.00000 

2.00000 

1.00000 

1.00000 

0.50000 

0.41383 

1.00000 

2.00000 

1.00000 

1.00000 

0.62500 

0.48281 

1.00000 

2.00000 

1.00000 

1.00000 

0.75000 

0.54348 

1.00000 

2.00000 

1.00000 

1.00000 

0.87300 

0.60243 

1.00000 

2.00000 

1.00000 

1.00000 

1.00000 

0.63416 

1.00000 

2.00000 

1.00000 

1.00000 

1.12300 

0.69741 

0*99740 

1.99320 

0.99369 

0,99320 

1.23000 

0.72787 

0.99183 

1.98369 

0.98077 

0.98369 

1.37300 

0.73098 

0.98387 

1.97173 

0.96012 

0.97173 

1 .30000 

0.76783 

0.97962 

1.93925 

0.93487 

0,93923 

1.62300 

0.77938 

0.97312 

1.94623 

0.90638 

0.94623 

1.73000 

0.78629 

0.96641 

1.93261 

0.87539 

0.73281 

1.87300 

0.78918 

0.93936 

1.91913 

0.84284 

0.91913 

2.000C0 

0.78839 

0.93263 

1*90330 

0.80933 

0.90530 

2.12300 

0.78498 

0.94574 

1.89147 

0.77334 

0.89147 

2.23000 

0, >7877 

0.93888 

1 .97776 

0.74203 

0.87776 

2.37300 

0.77036 

0.93213 

1 .86425 

0.70933 

0.86425 

2.300C0 

0.76008 

0.92353 

1 .83106 

0.67744 

0.85106 

2.^':?300 

0.74829 

0.91913 

1.83823 

0.64696 

0.63823 

2.75000 

0.73328 

0.91294 

1 .82589 

0,6177’ 

0.82589 

2.87300 

0.72133 

0.90701 

1 .61402 

0.39023 

0.81402 

3.00000 

0.70679 

0.90133 

1 .60266 

0.36423 

0*80266 


Figure 68 - Continued 
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oRicmra. 

OF POOR QOkLi. . 


PHI> 45.0000 0E0<^ 


THEFF- 2*00000 

TACFF* 1 

*34541 



R 

UN/VN 

TH/TM 

TN/TA 

DPH/OPN 

OTM/DTM 

0.0 

0*0 

1.00000 

2.00000 

1.00000 

1.00000 

0.17078 

0.19441 

1.00000 

2.00000 

1.00000 

1 .00000 

0.3S3S9 

0.32347 

1.00000 

2.00000 

1.00000 

l.OOOCO 

0.S3033 

0.43124 

1.00000 

2.00000 

1 .00000 

1.00000 

0.70711 

0*32440 

1.00000 

2.00000 

1 .00000 

1 .00000 

0 *88389 

0.40428 

1.00000 

2.00000 

1.00000 

1.00000 

1 .OAOAA 

0.47739 

1.00000 

2.00000 

1.00000 

1.00000 

1*23744 

0.72490 

0*99149 

1*98298 

0.9B242 

0.98290 

1*41421 

0*73432 

0*98195 

1.94390 

0.95298 

0.9439O 

1.5900? 

0.77378 

0.97180 

1.94340 

0.91485 

0.94340 

1*74774 

0*78334 

0.94117 

1.92233 

0.87140 

0.92233 

1 *94454 

0.78744 

0.95024 

1.90051 

0.82494 

0.90051 

2.12132 

0.78287 

0.93928 

1.87857 

0.77744 

0.87057 

2.29609 

0*77309 

0.92844 

1.85488 

0.73037 

0.05486 

2*47467 

0.75920 

0*91/88 

1.83577 

0.48473 

0.83577 

2*45144 

0*74213 

0*90773 

1.81550 

0.44144 

0.81550 

2.82842 

0.72284 

0*69813 

1.79424 

0.40132 

0.79424 

3.00519 

0.70210 

0.88909 

1.77810 

0.54415 

0.77818 

3.18197 

0*48048 

0.88044 

1.74132 

0*53000 

0.74132 

3.35873 

0.43927 

0.87283 

1.74570 

0.49919 

0.74570 

3*33552 

0.43844 

0*84344 

1.73128 

0*47101 

0.73128 

3.71230 

0.41841 

0.85898 

1.71/94 

0.44550 

0.71794 

3.88907 

0.59812 

0.85242 

1.70404 

0.4194V 

0.70484 

4*04585 

0*54844 

0.84472 

1.40944 

0.38781 

0.40944 

4.24243 

0.54030 

PHI* 89 

0.83824 
.9999 PECS. 

1 .47410 

0.34826 

0.47448 


TNEFF- 2.00000 

TAEFF* 

1.21417 



8 

VH/VN 

TH/TM 

TM/TA 

DPH/DPN 

DTM/DTM 

0.0 

0.0 

1.00000 

2.00000 

1.00000 

t .00000 

0.25000 

0.25114 

1.00000 

2.00000 

1.00000 

1.00000 

0.50002 

0.41383 

1.00000 

2.00000 

1.00000 

1.00000 

0.75000 

0.54548 

1.00000 

2.00000 

1.00000 

1.00000 

1.00000 

0.45414 

1.00000 

2.00000 

1.00000 

1 .00000 

1.25000 

0.72490 

0.98828 

1.97454 

0.98127 

0.97454 

1.50000 

0.74518 

0.97072 

1.94144 

0.93417 

0.94144 

1.75000 

0.78159 

0.95177 

1.90354 

0.87750 

0.90354 

2.00000 

0.78170 

0.93210 

1 .84419 

0.81215 

0.84419 

2.25000 

0.74971 

0.91244 

1.82472 

0.74534 

0.82492 

2.50000 

0.74898 

0.89349 

1.78499 

0.40107 

0.78499 

2.75000 

0.72238 

0.B7S47 

1.75135 

0.42130 

0.73135 

3.00000 

0.49239 

0.85929 

1.71858 

0.54789 

0.71858 

3.23000 

0.44122 

0.84449 

1.48897 

0.32032 

0.46897 

3.50000 

0.43070 

0.83124 

1.44249 

0.47881 

0.44249 

3.75000 

0.40214 

0.81944 

1.43888 

0.44249 

0.43888 

4.00000 

0.54724 

0.80454 

1.41307 

0.39897 

0.41307 

4.25000 

0.52458 

0.79499 

1.58998 

0.34880 

0.58998 

4.50000 

0.49105 

0.78474 

1.54952 

0.30724 

0.54932 

4,75000 

0.45974 

0.77544 

1.35128 

0.27252 

0.53128 

5.00000 

0.43202 

0.74744 

1.53491 

0.24319 

0.53491 

5 25000 

0.40724 

0.74007 

1.52014 

0.21822 

0.32014 

5.50000 

0.38504 

0.75334 

1.50473 

0.19479 

0.50473 

5.75000 

0.34499 

0.74725 

1.49450 

0.17828 

0.49450 

4.00000 

0.34482 

0.74144 

1.40331 

0.14219 

0.48331 


Figure 68 - Continued 


116 



OF POOR 


PHI-134. 99f? DEQS. 


THEFF- 2.0000e TAEfF- 1.04273 


8 

VH/OH 

TM/TN 

TH/T4 

DPN/PPN 

0.0 

0.0 

1.00000 

2.00000 

1 .00000 

0.23000 

0*23114 

1.00000 

2.00000 

1 .00000 

0.30000 

0.41363 

1 .00000 

2.00000 

1 .00' . 

0.73000 

0.34348 

1.00000 

2.00000 

1.00000 

1 .00000 

0.45414 

1 .00000 

2.00000 

1.00000 

1.23000 

0.72414 

0.98333 

1.97UI 

0.96144 

1 .50000 

0.74313 

0.94392 

1.92784 

0.93718 

1.75000 

0.77798 

0.94040 

1.88119 

0.87914 

2.00000 

0.77441 

0.91442 

1.83285 

0.81434 

2.25000 

0.74272 

0.89233 

1.78470 

0.74799 

2.30000 

0.74041 

0.84914 

1.73833 

0.48397 

2.75000 

0.71242 

0.84743 

1.49484 

0.42450 

3.00000 

0.46128 

0.82750 

1.45501 

0.57087 

3.23000 

9.44917 

0.80954 

1.41908 

0.32317 

3.30000 

0.41792 

0.79331 

1.5B702 

0.48147 

3.75000 

0.38883 

0.77923 

1.35850 

0.44494 

4.00000 

0.35341 

0.74373 

1.32750 

0.40126 

4.25000 

0.51303 

0.74993 

1.49987 

0.35097 

4.50000 

0*47748 

0.73774 

1.47547 

0.30930 

4.75000 

0.44443 

0.72488 

1.45374 

0.27443 

5.00000 

0.41913 

0.71714 

1.43433 

0.24497 

5.25000 

0.39447 

0.70840 

1.41481 

0.21988 

5.50000 

0.37274 

0.70047 

1.40093 

0. 19835 

3.75000 

0.3529* 

0.49325 

1.30451 

0.1797^ 

4.00000 

0.33312 

0.48444 

1.37331 

0.14355 


PHI 

*179.9999 PEGS. 




THEFF 

2.00000 

TAEFF* 

1.00000 


R 

VH'VH 

TH/TN 

TH/^h 

PPH/PPM 

0.0 

0.0 

1.00000 

2.00000 

l.OCOOO 

0.25000 

0.25114 

1.00000 

2.00000 

1.00000 

0.30000 

0.41383 

1.00000 

2.00000 

1 .00000 

0.73000 

0.54548 

1.00000 

2.00000 

1 .00000 

1.00000 

0.45414 

1.00000 

2.0000C 

1.00000 

1.25000 

0.72383 

0.98437 

1.96874 

0.981B3 

1.50000 

0.74224 

0.94097 

1 ,92194 

0.937.2 

1,75000 

0.77441 

0.93375 

1.87149 

0.97987 

2.00000 

0.77409 

0.90943 

1.81925 

0.81533 

2.25000 

0.75947 

0.08344 

1.74727 

0.74917 

2.50000 

0.73447 

0.03843 

1.71724 

0.48527 

2.75000 

0.70807 

0.83522 

1.47044 

0.42565 

3.00000 

0.47441 

0.81370 

1 .42735 

0.5/222 

3.23000 

0.44390 

0.79447 

1.58893 

0.52447 

3.50000 

0.61233 

0.77725 

1.55450 

0.46248 

3.75000 

0,58299 

0.76196 

1.52392 

0*44404 

4.00000 

0.54743 

0.74536 

1.49073 

0.40235 

4.25000 

0.50710 

0.73061 

1.44121 

0.35197 

4.50000 

0.47J83 

0.71759 

1.43518 

0.31024 

4.73000 

0.44068 

0.70602 

1.41203 

0.27531 

5.00000 

0.<1352 

0.69567 

1.39135 

0.24580 

5.25000 

0.38916 

0.68636 

1.37271 

0.22044 

5.50000 

0.34736 

0.67792 

: .35503 

0.19907 

3.75000 

0.34774 

0.67026 

1.34051 

0.18042 

4.00000 

0.33000 

0.66325 

1.32651 

0.1641V 


Figure 68 - Continued 


oth/ptn 

I .00000 
1.00000 
1.00000 
1 .00000 
1 .00000 
0.97111 
0.927%4 
0.80119 
0.83203 
0.*>e470 
0.73833 
0.49404 
0.4SS01 
0 41903 
0.SP702 
0. 33830 
0.52730 
0.49987 
0.4734/ 
0.43374 
0.42433 
0'.'41481 
0.40093 
0.38431 
0.37331 


DTH/DTN 

1 .00000 
1 . 00000 
1.00000 
t. 00000 
1.00000 
0.94874 
0.92194 
0.87149 
0.81925 
0.74727 
0.71726 
0.47044 
0.42753 
0.33893 
0.33430 
0.52392 
0. 49073 
0.44121 
0.43318 
0.41203 
0.39135 
0.37271 
0.3;>585 
0.34051 
0.32iSl 


117 




OF FCOr. Cc 


rAXimin GROJ:'*? r-RtSSURe and TCHP^RATURE ALOKO UFbASHaROUND STAGNATIOK LINE 


XU/RN 

XU/ S 

OPNAX/OPJ 

PHAX 'PHAXO 

DTNAX/DU 

tkax/tnaxu 

0.0 

0.0 

0.36423 

^ . 00000 

0.S0266 

1.00000 

0.17473 

C. 02912 

C. 56128 

0.99478 

0.80208 

0 99927 

0.35063 

0.03844 

0.55263 

0.97944 

0 80071 

0.99707 

0.S269G 

O.OS8M 

0.53952 

0.95443 

0.79734 

C-99337 

0.71101 

C. 1«850 

0,51902 

0.919B7 

0./9313 

0.98815 

0.89814 

0.14949 

0.4918? 

0.87701 

0.78771 

0.98137 

1.09171 

0.18198 

0.^.6421 

0.32629 

C. 78096 

C. 97296 

1 .294C7 

0.21368 

0. 4338 J 

0.76884 

0.77286 

0.96287 

1.50663 

0.23111 

9.39862 

0.706-«8 

0.76333 

0.95104 

1.73203 

o.r&a6v 

0.36127 

0.64028 

0.75241 

0.93739 

i .97313 

0.32835 

0.32275 

0.57202 

C . 73994 

C. 92186 

2.23341 

0.37224 

0.28402 

0.50337 

0.72591 

0 90438 

2.31770 

0.41955 

0.24282 

0.45036 

0.70073 

0.88298 

2.83035 

0.47172 

0.19617 

0.34768 

0.68737 

0.85636 

3.17981 

0.52997 

0.15246 

0.27198 

0.66516 

0.8286’> 

3.5^523 

0.59568 

0.11645 

0.20639 

C. 64118 

0.79882 

4.02969 

0.67162 

0.08533 

0.15127 

0.61545 

0.76673 

4.36127 

0.7602; 

0.06009 

0.10630 

0.58/96 

0.73251 

3.19614 

9.86602 

0.04024 

0.07149 

0.55374 

0.6961! 

5.97346 

0.99558 

0.0255 

0.04531 

0.52T81 

0.65757 

6.95473 

1 . 15912 

0.01506 

0.02673 

0.4952C 

0.61 694 

8.24235 

!. 37372 

0.0C811 

0.01438 

0.46293 

0.57425 

12.02036 

i. 67009 

0.00385 

0.00663 

0.42503 

0.52752 

12.63775 

2.10962 

0.00153 

0.00272 

0.38750 

0.<8277 

17.01334 

2.83536 

0.00046 

0. 00062 

0.34829 

0.43392 
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Definitions 


Output Titles 

RJH/RN Half-velocity radius of jet measured from jet center- 
line 

R J /RN Half width of jet 

ALPV Exponent of velocity profile 

VJ/VN Centerline velocity 

RC /RN Potential core radius 

ALPT Exponent of temperature profile 
RJNT/RN Half temperature radius 
TM/TN Center line temperature non-dimensionalized by nozzle 
temperature 

TM/TA Centerline temperature nondimensionalized by ambient 
temperature 

DTM/DTN (TM-TA)/(TN-TA) 

RJHQ/RN Half dynamic pressure radius 
QJ/QN Centerline dynaraic pressure 

Note: RN , VN . TN , and QN denote nozzle radius . velocity , 

temperature, and dynamic pressure. TA denotes ambient 
temperature. 


Jet Deflection Ret 


Single jet impingement characteristics are printed. 


Output Titles Definitions 

RGH/RN Hsdf- velocity radius of jet at ground-effect height 
RG/RN Jet half width at ground-effect height 
RO/RN Deflection zone or pressure recovery radius 
VG/VN Square root of jet g^round stagnation pressure 
DPS/DPTJ Jet ground stagnation pressure nondimensionalized by 
nozzle stagnation pressure 

TGS/TA Stagflation temperature of incident jet at gfound-effect 
height . 

ALPG Exponent of gpround pres.sure distribution 
R/RGH Radial location in deflection region relative to ground 

stagflation point. Nondimensionalized bjr ground-effect 
half-velocity radius of jet 





DPS/DPTG Ground pressure nondimensionalized by ground 
stagfnation pressure 

R/RN Radial location in deflection region nondimensionalized 

by nozzle radius 

DPS/DPTJ Ground pressure ncndimensionalizad by nozzle 

stagnation pressure 

Wall Jet Region 

Isolated wall jet properties are printed in this set of output. 


Output T itles 
DELS/RN 
VMI/VG \ 
VMI/VN i 
ALPWO 
BV/OH/RN 
BWO/RN 
PMI/PR 
TMI/TA 

(PTTOT/PR-l.) 

(TOT/TA-1.) 

R/RN 

VM/VG 

VM /VN 

BWH/RN 

BW/RN 

DELBL/RN 

DELPJ 

KDEL 

TM/TN ♦ 
TM/TA I 

BWNT/RN 

DTM/DTN 


Definitions 

Boundai^ layer thickness at stagnation point 

In viscid maxinmn) velocity at jet half width (RGH) radial 

location or start of transition to turbulent wall jet 

Exponent of velocity profile at start of wall jet 

Half velocity thickness at start of wall jet 

Initial thickness of wall layer 

Initial static pressxire in wall layer 

Initial temperatiure of wall layer 

Initial total pressure of wall layer 

Initial total temperature of wall layer 

Radial location in wall jet measured from ground 

stagnation point 

Maximum velocity in wall layer nondimensionalized by 
square root of groxind stagnation pressure 
Maximum velocity in wall layer referenced to nozzle 
exit velocity. 

Half velocity thickness of wall layer 

Thickness of wall layer 

Boundary layer thickness 

Maximum stagnation pressure in wall layer 

Ratio of boundary layer to total wall layer thickness 

Maximum temperature in wall layer profile 

Half temperature thickness of wall layer 
(TM-TA)/(TN-TA) 
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TWO-JET IMPINGEMENT IN TERAC T ION OUTPUT 


Recirculation Effect on Wall Jets 

Note: This set of output recomputes the wall jet properties accounting 
for the interaction of two hot jets. 


Output Titles 
PHI 

TNEFF ( 
TAEFF i 
R 

VM/VN 

TM/TN ( 
TM/TA j 

DPM/DPN 

DTM/DTN 


Definitions 

Azimuthal angle in jet-centered ground polar coordinate 
system 

Effective nozzle and ambient temperature ratios with 
recirculation effects taken into account. 

Radial location in wall jet , measured from ground 
stagnation point 
Maximum velocity in wall layer 
Maximum temperature in wall layer 

(PM-PA)/(PN-PA) or (PM-1.)/(NPR-1.) 
(TM-TA)/(TN-TA) or (TM-1. ) /(NTR-1. ) 


Maximum Ground Pressures along Upwash Ground Stagnation Line 


Output Titles 
XW/RN ^ 

Definitions 

Coordinate along stagnation iine measured from upwash 

XW/S \ 

stagnation point 

DPMAX/DPJ 

Stagnation line pressure nondimensionalized by nozzle 
stagnation pressure 

PMAX/PMAXO 

Stagnation line pressure nondimensionalized by upwash 
stagnation point pressure 

DTMAX/DTJ 

Stagnation line temperaxure nondimensionalized by nozzle 
stagnation temperature 

TMAX/TMAXO 

Stagnation line temperature nondimensionalized by 
upwash stagnation point temperature 

Upwash Momentum Function 

Output Titles 

Definitions 

RG/RN 

Radius of jet at ground effect height 

AGON 

Constant in momentum function 



XMQMZ 


Total vertical monientum In up wash sheet non-dimen- 
sionalized by the optimum value of Mj/2tt 

PH 10 Coalescence angle used in upwash momentum model 

If PKIO = 0 jets are not coalesced 
If PHIO > 0 jets have begun to coalesce 

Note: Depending on the value of PHIO (i.e. zero or non-zero) the constant 
ACON applies to the appropriate upwash momentum model. 

CU Nondimensional upwash deflection region width constant 

Note: There are two possible outputs that can occur at this point. If the 
i»ts are spaced far enough apart the comment: 

JET AND UPWASH DEFLECTION REGIONS DO NOT 

INTERACT. CU ESTIMATE IS CORRECT 

In this case the value of CU is correct and the jet impingement and 
upwash deflection regions are independent. The perturbation 
parameters are then defined as: 

EPS = 0.0 
SIG =1.0 
PHIUO = 0.0 
PHIO = 0.0 

PMIN =0.0 along upwash line 
ALPUG = 1.50 

If the nozzle spacing and height above ground are such that the deflec- 
tion regions interact, the subroutine INTERG will be called and 

♦♦♦CALL INTERG ♦♦♦ 

will be displayed. The following output will be printed. 

Note: The ground pressure distributions are computed without temperature 
effects. 

Output Titles Definitions 

Iteration cycle Number of iterations required to find solution using 
Newton's method 
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EPS or EPI Perturbation parameter for jet impingement pressure 
distribution 

SIGMA or SIG Perturbation parameter for upwash deflection region 
pressure distribution 

CU or cm L'pwash width estimate prior to iterative solution 
PHIU'O Angles defining the intersection of jet and upwash 

PH 10 deflection regions 

Note : The following values apply along the line connecting the jet stag- 
nation points on the ground (i.e. x - 0) . 

DELPWO Pressure at upwash stagnation point nondimensionalized 

by nozzle stagnation pressure 

PMIN Minimum pressure between jet and upwash deflection 

regions nondimensionalized by nozzle stagnation 
pressure 

ALPUG Exponent of upwash ground pressure distribution 

function 

(CU) X 

(SIGMA) Final value of upwash thickness constant 

Note: If IPBAR - 0, the following output will occur. 

Computation of Two-Jet Ground Isobar Pattern 

Output Titles Definitions 

PBAR Input values of pressure, nondimensionalized by nozzle 

stagflation pressures, will be echoed. 

Note : The following will be repeated NU times. 

Output Titles Definitions 

IJET Number of points on jet impingement region isobar 

lU Number of points on upwash deflection region isobar 

If IJET = 0 or lU = 0, the specified value of PBAR was not found 
in ground pressure distribution. 

Note : The following coordinates are referenced to the jet stagnation point 
on the ground. 

Output Titles Definitions 

XISOJ X coordinate of jet isobar 

YISOJ Y coordinate of jet isobar 
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XISOU X coordinate upwash region isobar 

YISOU Y coordinate of upwash region isobar 

The final set of output in this section is the upwash deflection zone 
line . If the deflection zones do not hiteract , ambient conditions exist 
along this line. 

XUP X coordinate of upwash line 

YUP Y coordinate of upwash line 

COMPUTATION OF UPWASH FLOW FIELD 


Upwash Streamline Properties 

Note: Two streamlines are printed 


Output Titles 
PH ID 

INV. VEL. 

TMI/TA 

DELPWJ 

BWGH/RN 
ZOU /RN 
BUO/RN 
BUOHN /RN 


Definitions 

Azimuthal angle of upwash streamline referenced to jet 

ground coordinate system 

Inviscid turning region maximum velocity 

Inviscid turning region maximum temperature 

Maximum pressure on upwash stagnation line on the 

ground where upwash streamline orginated 

Half velocity width of incident wall jet streamline 

Upwash turning region height above ground 

Initial upwash width 

Initial upwash half velocity width 


Upwash Streamline Decay Properties 


Output Titles 

Z/RN 

VMU/VN 

TM/TN \ 
TM/TA j 

BUH/RN 

BU/RN 

DPU/DPJ 


Definitions 

Upwash Streamline coordinate (ZS) 

Upwash maximum or centerline velocity 
Upv;ash maximum or centerline temperature 


Half velocity width 
Half width of upwash 

Upwash maximum or centerline total pressure non- 
dimensionali?ed by nozzle stagnation pressure or 
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Upwash Properties Computed At Z Location of Underside of Body 


If ZB = 0, the output will yield the upwash properties at Z = constant 
plane above ground. 


Output Titles 
ZU/RN 
RW/RN 

PHIB 

X/RN 

VMU/VN 

DPU/PN 


BUH/RN 

DPU/PUO 


DPU /PJS 

TM/TA 

DTSM/DTN 

PBS/PN 


Definitions 

Upwash coordinate measured from groiuid plane 
Radial coordinate from jet ground stagnation point to 
upwash stagnation line 

Azimuthal location of upwash streamline referenced to 
jet coordinate system 

X coordinate in upwash sheet or X coordinate on 
fuselage 

Maximum or centerline upwash velocity 
Maximum or centerline upwash total pressure 
nondimensionalized by nozzle stagnation pressure 




Half-velocity width of upwash 
Upwash total pressure nondimensionalized by total 
pressure on streamline originating from upwash stag- 
nation point (X = 0, Y = 0) 

Upwash total pressure nondimensionalized by jet 

ground stagnation pressure 

Maximum or centerline upwash tempeiature 

(TM-TA)/(TN-TA) 

Stagnation pressure on underside of fuselage placed 
in upwash sneet 


Upwash Lift Force 

Two values are printed; 


FLAT BOTTOM VEHICLE WITH SHARP CORNERS 
BODY WITH CIRCULAR CROSS SECTION 


Output Title Definition 

LU/2TJ Upwash lift force nondimensionalized by the total 

thrust of the two jets 

COMPUTER PROGRAM LISTING 

Figure 69 sliows a Fortran listing of the computer program 
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U U u CJ 


^1. pi\:h 




c jt:ooo:o 

C VTOl TUO'JET inFlMCEHEMT INTERACTION COMPUTER PROGRAM JCT00030 

C FOR CLOSELY SPACED JETS JET00040 

C BY H*J.SiCLARl JETOOOSO 

C 5U 575-2207 JET00060 

C GRUNHAN AEROSPACE CORPORATION JET00070 

C B£TKPAGE»NEU YORK I17M JET00080 

C JET00090 

C JETOOUO 

C PROGRAM GRUHH0T2 JE TOO I 20 

C JET00130 

REAL KDELFDtKDELOfNFDfNOtNiKDEL JETOOMO 

REAL HACHfNPRtNTR JETOOI50 

COHHON/UALL/DELSfNOiKr«ELO*ALPG»ALPUO>ALPUFDtROfRGH>UGtBMHO JETOOUO 

COMHON/UALLT/ TLAHl »DELPS»DELTSDt PSDAR JET00170 

COHNON/MOZZLE/NPR«NTR»TNfGANfMACHfACHSfTAEFFrTNEFF JET00180 

CONHON/ JE T/SD t RCU JE TOOT 90 

COHMON/HEIGHT/ HDtIPBAR JET00200 


PROGRAM GRUHH0T2 


REAL KDELFDtKDELOfNFOfNOtNiKDEL JETOOUO 

REAL HACHfNPRtNTR JETOOT50 

COHHON/UALL/DELSfNOiKr«ELO*ALPG»ALPUO>ALPUFDtROfRGH>UGtBMHO JETOOTOO 

COMHON/UALLT/ TLAHl »DELPS»DELTSDt PSDAR JET00170 

COHMON/MOZZLE/NPR«NTR»TNiGANfMACHfACHSfTAEFFrTNEFF JET00180 

COMMON/ JE T/SD t RCU JE TOOT 90 

COMMON/HEIGHT/ HDtIPBAR JET00200 

FCTAUD(ETAtA)-ETA-(4./(AfT. ))IETAM(AfT. )T(6./(2.yATT.)>4ETA4t(2.IJET002T0 
TAFT.)-<4./(3.>ATT. )>IETAl4(3.lAfT.)fa./i4.«ATT. ))4ETASy(4.4AH.) JET00220 

JET00230 

BODY SHAPE FUNCTIONS FOR FORCE PREDICTION IET00240 

y2...HALF-UIDTH IN NOZZLE RADIT JET00250 

ZDODY...HALF DEPTH IN NOZZLE RADII JET002^0 

DZDXD... AXIAL SLOPES JET00270 

THSOSE..,NObE ANGLE OF PARADOLTC BODY JET002S0 


C parabolic body UlTH CONSTANT CROSS SECTIONAL SHAPE 
C 

C W2<XB)^.5tTAN(THMflSE)4<XL2»l2-XBtl2)/XL2 

C ZBODY(XB)-.54TAH(TnNOSE)|(XL24l2-XB4t2;/XL2 

I DZDXB ( XB ) TAN< THNOSE ) « X B/XL2 

C 

C CYLINDRICAL BODY UITH CONSTANT CROSS SECTIONAL SHAPE 
C 

U2aB)=UC0N 

ZBODY<XB)=ZCON 

DZDXB(XB)=7SC0N 

C 

C CYLINDRICAL BODY PARAMETERS IK TERMS OF NOZZLE DIAMTERS 


ZSCON^O.O 

YHN08E^14.03A2«PI/180. 

CtRCULAA CROSS SECTION 

POINT Of FLOy SEPARimON 

PI>3.141S924S 

QAM>T«4 

FHISEP-135.4PI/TB0. 

READ (5>T> NPRtNTR 
READ (5>T) H0*8D»ZPLDiDZPLfZFlNAL 
T FORMAT (5F10.5) 

READ (5rl> XL2rMC0NtZC0N 
WRITE (Af2U 

21 FORMAT ( ' 1 ' »///30X» 'Fit INPUT PARAMETERS ttt' 
M^CM•0nRT^2. t(MPR-l . I/OAM) 
rN*N?R 

WRITE NPRfNTR 

17 FORMAT <20X» 'MPR»' »F8.4»2X» 'WTR«* ,F8.4) 

WRITE (6i27> MACHfTN 

27 FORMAT <20X»'HACH HO . • ' i F 1 0 . 4 . 2X t ' TN/T A* ' »F TO . 4 > 
WRITE <A»9) HDfSD.ZFLD 

9 FORMAT « 20’'# 'H/0» ' #FJ0.5»2X » S/0*' »Kld.3 »2Xf ' 7/0* 
WRITE CA.907) ZPLD*DZPL*ZFtmaL 

R8 1-1622-0640(30 Sh««tJ) 


Figure 69 Program FORTRAN Listing 
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OF p - . :,.;ALrTY 


c 

c 

c 

c 


c 

c 

C 

C 

C 

C 

C 

C 


907 FORHAT ( /20X » ' 2/D» ' t F 1 0 . 5 r 2X » ' DZ/D* ' » F 10 . 3 i 2X f ' ZF INAL/0* ' »FtO .3/ ) JET00700 
UKire (6»908) XL2>UCON»2CON J1T00710 

008 FORMAT (/20Xf'|i0DY LENGTH L /D« ' » F lO . 3 » 2X f ' WIDTH U/0« ' »F 10 *3 » 2X • ' DC JET00720 


IfTH ZD/D*' »F10.5/> 

ZPLANE^2.iZPL0 

CQHPUTE JET standoff DISTANCE AND EFFECTIVE JET DECAY HEIGHT 
(GIRALT ET AL) 

IF iHD.LE.&.8> DELG«0. 133«7,8 
IF (HD.OT.6.8) DEI.G«0. 1S3«( 1 . fHD> 

ZEND‘(HD>.30tDEL0)«3. 

WRITE (6»22) 

22 FORMAT < ///30X t ' Ft JET DECAY REGION 
WRITE .6*7) DELG 

7 FORMAT (lOX.'JET DEFLECTION HEIGHT DELG/D ='fFl0.3) 

COMPUTE JET DECAY 

ALl DIMENSIONS NONDIMENSI AL IZED $Y NOZZLE RADIUS 

SET EMPIRICAL CONSTANTS FOR HALF-WIDTH BEHAVIOR AND LENGTH OF 
POTENTIAL CORE (HRYCAK»ET AL> 

A2*.040 
•2*. 000 
FH»^l. 

A3->«09$(1 .4AL0C(2./< 1 . ^1 . /TN«9 .3 > ) ) t ( 1 . 0- . I AtFHAMACH ) 


JET00730 

JET00740 

JET00730 

JET00740 

JIT00770 

JET007B0 

JET00790 

JETOOaoO 

JETOCaiO 

JCT00820 

JET00830 

JETOOB40 

JET00830 

JETOOBAO 

JET00870 

JET00880 

JET00890 

JET00900 

JET00910 

JET00920 

JET00930 

JET00940 

JET00950 

JET009&0 

JET00970 


C 

r 

c 


B3»0. 

ZCCQ-8.55 

FM«0.0 

IRC* .3$ZPC0«( 1 *41 ./(TN94.73) ) t( 1 .04 . 1A*FH9HACH> 

Al-A24(B2-1 . )/ZPC 

ZFD«B2/<A3-A2> 

ALPFD«1.S 

WRITE (6*23) ZPCfZFD 

23 FORMAT M OX * ' ZPC/RN- ' , F 1 0 . 5 * 2X r ' ZFD/RN« ' *F 10 . 3 ) 

WRITE (4*12) 

12 FORMAi <//3X. ' Z/RN' *3X* 'RJH/RN' *3X» 'RJ/RM' t2X* ' ALPV' 
. # RC/RN' *3Xi ALPT' *3X. ‘R JHT/RN ' *2X* ^TH/TN' r2X» TM/TA 
2f 2X* 'RJMO/RN' .2X» 'OJ/ON' ) 

REGICN T. . .POTENTIAL CORE REGION OF JET.. 

NJPC»3 

0MPC*NJPC-1 

IF CEND.GE.ZPC) DZPC*ZPC/DNPC 
IF (ZEND.LT.ZPO PZPC«ZENP/DNPC 
Z=0. 

DO 100 I I .NJFC 

IF ( I .EQ. I ' 60 TO 102 

GO TO 103 

102 VJ»l. 

RJH^l . 
ftJHT*! . 

RJ-1 . 

ALPVC»0. 

ALPTC-1 .3 

CV2*.3 

RJPC«1. 

TMal. 

DELFMN*! . 

THB-TN 
DELTMN»1 .0 
RC«I. 

RJHQ-1 . 

GO TO 104 

103 ETAC«l , -Z/ZPC 
RJM»AI*241 . 

CALL JETPCT ( TN*ETAC tRJH* ALPVC * TLAMf CV2»CT ,RV»RT ) 
TLAMPC"T1 :.r. 

RJ*RJM/Rv 
R JMT »TLAMtRJM 
TM-l . 

TMB*TN 
DELTMN-I . 

DELFMN*! . 


JET00980 
JET00990 
JETOIOOQ 
JETOiOlO 
JET01C2Q 
JET01030 
JET01040 
JET01030 
JET01060 
JET01070 
t4X* 'VJ/VN' *4XJET01080 
' *4X* 'DTH/DTN'JEf 01090 
JETOllOO 
JETOlllO 
JET0U20 
JETOnSO 
JET01140 
JET01150 
JET01160 
JETOU70 
JET01180 
JET01190 
JET01200 
JET01210 
JET01220 
JET01230 
JET01240 
JET01230 
JET01260 
JET01270 
JET0128Q 
JET01290 
JET01300 
JET01310 
JET01320 
JET01330 
JET01340 
JET01300 
JETOnoo 
JET01370 
JET01360 
JET01390 
JET01400 
JET01410 
JET01420 
JET01430 
JET01440 
JET01430 
JETOl AAC 


Figure 69 - Continued 
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U O u U CJ (J 


Ci V* ■■ 


\S 


( ^ » «*fV 


RJF‘C»#tJ 

RJHTPC*RJHT 

VJ-l. 

RC-ETACtRJ 

AtRTC-ALRMC 

CALI QHALF <TH9*ETAC»ALPVC»7LAHtRa) 

RJHQ«RJ4RQ 

M0«1. 

T0«t. 

T«STAO»rHt*VJRVjt(KTR-TM) 

DttTA0-(TH8TAQ-l . )/ (NTR-i. > 

t04 URITE U»tl> ZtRJH»RJrALPMC»VJ»RC»ALF'TC*RJHT rTHt rnBfDCLTt1N»RJMQ> 
lOCLPNN 

11 FORHAT (IXr4(F7.4»lX>»A<F7.4«tX> i4<FBM*U>> 

Z-240ZPC 

100 continue 

TRANSITION RCOION II ANO FULCY DCUCLQFEK REGION OF JET 

IF (ZEND.LE.ZFC) 00 TO 

NJFO^ZS 

0NF1V«NJFD-1 

ALPVFD*1.S 

ALPTFD«1.50 

TLAHFD>1.4. tes«TN«i.5 

RC0N-<2.-S0RT<2.))/'. 

I RVFD«RC0N«4(i*/ALPMF0) 

RJF0*(A3$ZFD4»3)/RUF0 

RTFD«RVrI> 

RJHTFD«TLAH» I •(AS4ZFDFB3) 

DZ«(ZEND-ZPO/ONFO ,, 

Z-ZPC4DZ 

00 200 I*=2fNJFD 

IF <Z.0T.ZFD) 00 TO 201 

RJ-RJPC4( (RJFP-RJPC)/(ZFO-ZPC> >t(Z-ZPC) 

RJH«A24Z4B2 

RV-RJH/RJ 

ALPV«ALOO ( RCON ) / ALOO (RV ) 

AlPT«AtPV 

TLAH«TLAHPC4UTLAHFD-TLANPC)/<ZFD-ZPC))4(Z-ZPC) 

00 TO 203 

201 ALPV«ALPVFD 
ALPT«ALPTFP 
RJH-A34Z4B3 
TLAN»TLA«FD 

202 CALL OTSOLV(TNtRJHiALPVrTLAN»UJ»TH»RU»RT) 

RJ«RJH/RU 

RJMT*TLAH$RJH 

THR-THRTN 

IF <A»S<TN-1 * ) .LT. 1 .E-S) DELTHN'l.O 
IF (ABt(TN-l . ) .OT. 1 ,E-5> PELTRN-l THB-t . >/(TM-l . ^ 

0EtPHN*VJ4VJ/TH 

CALL OMALF ( TflO » 0 . # ALF V p TL AHt Rfll 

RJHQ»RQ«RJ 

TG'TH 

VO-Uj 

TH8TA0«Tl1ifVJ4UJ*(NTR-TM) 

PTSTAO’^dHSTAO-! . )/(MTR-I . ) 

WRITE (4»U) Z»RJH*RJ»ALPV» VJf RCf ALF T rR JHT f TM. THRp DELTHNpRJHO 
1 f OELPnN 
Z-Z40Z 

200 CONTINUE 

JET DEFLECTION REGION 
JET IHPINGEHENT PRESSURES 


•F V01470 

JET01480 

JET01490 

JETOISOO 

JETOlStO 

JET01520 

JCT01930 

JET01940 

JET01930 

JETOISAO 

JCT0137C 

JET015BO 

JET01500 

JETOIAOO 

JETOlOlO 

JE701620 

JETOUJO 

JET01A40 

JET01630 

JETOIAOO 

JETC1A70 

JC701A80 

JETOlAfO 

JE TO 1700 

JET01710 

JET01720 

JETOr/30 

JET01740 

JET01730 

JET01760 

JET01770 

JET01780 

JETOk7?0 

JETOIBOO 

JET01810 

JET01820 

JET01630 

JET01840 

JETOieSO 

JET01B40 

JET01870 

JET01880 

JET01890 

JETOIOOO 

JETOIOlO 

JET01920 

JET01930 

JET01940 

JET01950 

JET01960 

JET01970 

JET01980 

JET01940 

JET02000 

JET02010 

JET02020 

JET02030 

IET02040 

JET02050 

JET02040 

JET02070 

JET02080 

JET02090 

JET02100 

JET02110 


998 A0C0N«-.03202A 
POCON^. 3354248 
C0C0N«2.8 

IF <HD.6T.3.» A0*3.4 
A0SLH*(3«4-2.9)/(3.-l .?) 

IF iMD.LE. 3.) A0*2.9+A0SLHt(H[»-l .2 > 

WRITE (4t79) 

79 FORHAT ( / //30X » ‘ tit** JET DEFLECTION REGION 

80 FORHAT ( I OX • ' RQH/RH* ' i F 1 0 . 5 1 2X # ' RG/RN« ' » F iO . 5 f 2X i 
I t 'VO/VN ' ' tFlO.S/ ) 


JET02120 
JET02130 
JET02140 
JET02130 
JCT02160 
JET02170 
JET02XBO 
JET02190 
ttttt'/) JET02200 

PO/RN* ' tFl0.3»2XJET022l0 
JET02220 
JET02230 


Figure 69 -• Continued 
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n nri nnnr>r»n 


ORIGINAL PACti 5S 
OF POOR QUALITY 


SET CHARACIERISTIC SCALE LENGTHS FOR IHF INOCNENT REGION 

RQH-RJH 

RO-RJ 

RO'AOtRGH 

ETAUOsRGH/RO 

IP (2.RHD.LE.ZPC) 00 TO 

deterhine stagnation point pressure 

434 WRITE <4»80> RGHtROiROtVO 

STAGNATION PRESSURE NONDINENSIONALIZED BY JET TOTAL PRESSURE 

ROU*RO ' 

OEtPS«<VGR«2>/TO 
PSBAR-1 .4DELPS«<NPR*I . ) 

ACHS*SQRT(2«t(PSBAR>l . >/CAH> 

tob-tortn 

T6BS*TOB4(NTR>TN>tVOBVO 
WRITE U>ei) DELPSpTGBS 

8t PORNAT (lOXt'STAONATlON pressure* DPS/PPT J- ' * F 10 . 3 »2X * 2X » ' GROUND 
lAX. TEHP.t TQS/TA«‘ tPlO.S) 

WRITE (4#eOU ACHS 

001 FORHAT </10X> 'STAGNATION MACH NUHBER OF JET- ' ' F tO . 3/ > 

CALL OPRES T I RO r UO* TO * ALPO f R * CU: ) 

WRITE Ut02) ALPS 
82 format <iOXr 'ALP0«* rFlO.S) 

NG-2S ^ 

RN0M«N0-i 

DR-RO/(RNGH$ROH) 

DCTA«L ./CRNOH) 

ETA-0. 

R«0* 

WRITE (Ai83) 

63 FORMAT (30Xf'SIN0LE JET GROUND PRESSURES ' //23X »' R/ROH '» 6X •' DPS/DP 
10* »8X? 'R/RN' i7X» 'DPS/DPTJV) 

DO 300 I-liNG 
DPT0-(l*-£TA«FALP0)t«4. 

DPTJ-DPTORDEl PS 
RRN-RfRQH 

WRITE (4*84> R»DPTOtRRN«DPTJ 
84 FORNAT < 20X t 4 ( F 1 0 . 3 * 5X ) > 

R-RIDR 

ETA-ETA4DCTA 
300 CONTINUE 

initialize wall JE ..EOION 

AW3-.09ta *>AL0e<2./< 1 .11 ./TM*$.5> > »»< 1 . - , lARACKS) 

TLAMF0*1.4.se3tTHf4.5 

ALPWFD-1.3 

KDCLFP-;0P42 

ALPUO-l.S 

NFD-7. 

CALL SlKM<KD£LFDrHfDiALPWPDiR4TW.CSV2UfCPW) 

CALL SIMU0<ALPUG#RATUG*CSPU9> 

CU-2 . «AW3tCSU2U/ ( RATUtCSPUG > 

ASl*. 013438 
AS2-. 128033 
PELS-2. tA51t(HD)$tAS2 
WRITE (6>2«) 

24 FORMAT (//y30Xt ^ f tt«t WALL JET REGION ttttt'//) 

WRITE (4*84) DELS 

84 FORMAT CSX* 'STAGNATION POINT BOUNDARY LAYER THICKNESS* DELS/RN-'* 
IFIO.3) 

NO-14 

KDEL0»DELS/(2*4D£L0) 

FETAO-U .>CTAWG«9ALP0)t«4 
PHBl-l . ♦(MPR-1 , )*OELPS*FETAn 
OELPSB»-(NPR-l . >*DELPS 
ACMI*SORT(2.4<PSBAR/P«BI-l . )/GAI1) 

TMBI-TGBS 

VMI- (ACM1/MACM)*50RT (TMBI/TM) 

OELTSB^-iGBS-l .0 

WRITE PMII .DELPSBfDELTSB 


JET02240 
JET0225O 
JET02260 
JET02270 
JET02200 
JET02290 
JET02300 
JET02310 
JCT02320 
JET02330 
JET02340 
JET023SC 
JCT02340 
JET02370 
JET02300 
JET0239Q 
JET02400 
JET02410 
JETQ2420 
JET02430 
JET02440 
MJET02450 
JET02440 
JET02470 
JET02480 
JET02490 
JET02500 
JET02310 
JCT02320 
JET02330 
JET02S40 
JET02330 
JET02S40 
JET02570 
JET02380 
TJE^02390 
JET. 00 
JET02e 
JCr02420 
JET02430 
JET02440 
JET02430 
JET02440 
JET02470 
JET02480 
JET02490 
JET02700 
JET02710 
JET02720 
JET02730 
JET02740 
JET027S0 
JET02740 
JET02770 
JET02760 
JET02790 
JET02900 
JET028IO 
JET02820 
JET0203Q 
JCT02840 
JET02850 
jrT02840 
JET02870 
JET02880 
JET02890 
JET02900 
JET02910 
JET02920 
JET02930 
JET02940 
JET02950 
JET02940 
JET02970 
JCT029&0 
JET02990 
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89 FORrtAT < IX. 'PHI/PA*^' .F10.5.2X. ' (PTTOT/PA-l . .F10.5.2X. ' (TOT/TA 
' n.)a'.F10.3) 

URITE <6.92» T - 

92 FOFMAT (SX.'THl i A« ' »F 10. 5 » 2X . ' V«I 'VH* ' »F 10 . 5) 

TH1-TH8I/TN .t 

CALL Gt1ATT( TN.F0H.DELPS.FN8I »Tf1 1 fVn I . ALPG.KDCLO.no. FTAMO. DELS I 
lALPUO.R.DUOiDUHO.TLAHI > 

RATUO«R 

DCLSH»8U0«KDCL0 
OH01«VH1/V0 
WRITE (4. ft) VKOl 

ft FORHAT (5X, 'VELOCITY AT START OF TURBULENT WALL JET. VH/V0«'.F10. 

n 

C WRITE (4.87) TLAHFO 

87 rORHAT (5X. 'TLAHFO-'. FIO. 3) 

WRITE (4.85) ALPWO.OWHO.BOO 

85 FORNAT (lOX. 'START OF WALL JET REOION. ALPW0-'.F10.3.2X.'iW0H/RN- 
1 .F10.S.2X. 'fPUO/RN> ' .FI 0.5) 

HAXIHUH VELOCITY DISTRIBUTION IN DEFLECTION REOION 

NWPTr23 

XNUPT-NWPT~1 

RCN0»2t$D 

OR-REND/XNWPT 

iU3HO«AW14RO 

DELiRH- ( OW3HO-OWHO > / ( R0-R6H > 

DELNR-(NFO*NC)/(RO-ROH) 

DELAR- ( ALPUFD-ALPWO ) / ( RO-ROH ) 

A0ELiT-.0i7S 

ADELiL*(ADELBT«(R0-2. )-DCLS)/<RO-ROH) 

R-0. 

TLAH8L«<TLAHF0>TLAHI)/(R0~RGH) 

WRITE (4.107) 

107 F0RHAT<//40X» ' ISOLATED WALL JET PROPERTIES'//) 

WRITE (4.97) 

97 FQRHAT ( //5X . ' R/RH ' »3X . ' VH/VG ' . 5X» ' VN/VN' . 3X . tfWH/RH ' . 2X » 

1 '8W/RN' »3X. 'DELIL/PN' .3X. 'DELPJ' .3X . ' KPEL ' . 3X ' TH/TN ' 

2f tX« 'TN/TA' .IX. 'BWHT/PN' . 1 X . ' DTM/DTN ' / ) 

DO 400 t»l>NUfT 

IF (R.OC.ROH) TLAM>TLAN!f TLAnSL«(R-ROH) 

IF (R.LT.ROH) TLAH^^TLAHI 

IF (R.LE.RO) FCTA^( 1 .'(R/R0>4fALPG)9f4 

IF (R.OT.RO) FETA^O.O 

IF (R.OT.RO) TLAH-TLAHFD 

PH8-1 , f (NPR>1 . )«DELPS9FCTA 

IF (R.OT.RGH) 00 TO 402 

ACH«S0RT(2*t(P9iAR/PHD-l . l/OAH) 

TN8«T0»S 

VHN-(ACH/HACH)$SQRT(THB/TN) 

TH-THD/TN 

OCLPJ-DELPS 

DELTA-KDELOtBWO 

N«M 

KDEL-KDELO 

»W-»WO 

DUH-RATU04DW0 
VHO*VHN/VO 
4WHT-TLAH199WH 
DCLTHN*(THB-1 . )/( TN-l.) 

TH8TAC«’TH»iVHNtVHNt(HTR‘*TN) 

DTSTAOXTHSTAO-l . )/<MTR-l. ) 

00 TO AOl 

402 IF (R.OT.RO) 00 TO 4)3 
l>UH*BUH04n£LiRHt(R-R0H) 

ALPU-ALPW040ELAR*(R‘R0M» 

DFLBL ‘PCL54ADELr..t(R>ROH) 

N-N04DCLNRt(R-R0M) 

00 TO 403 

403 DUN-AW3IR 
DELDL^ADtLBT4(R-2. > 

N-HFD 

ALPU-ALPWFD 
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-JET03000 
•ET03010 
JET03020 
JET03030 
JET03040 
JET03050 
JET03040 
JET03070 
JET03060 
JET03090 
JET03100 
5JET031 10 
JCT03I30 
JET03130 
JET03I40 
JET03130 
' JET03140 
JET03170 
JCT03180 
JFT03190 
JET03200 
JET03210 
JET03220 
JET03230 
JET03240 
JCT03250 
JET03240 
JET03270 
JET03280 
JET03290 
JfT03300 
JET033I0 
JET03320 
JCT03330 
JET03340 
JET033S0 
JCT03340 
JCT03370 
JCT03380 
JET03390 
JET0340O 
JET034 
JET03^ 0 
JET03430 
JET0J440 
JET03450 
JET03440 
JET03470 
JCT03480 
JET03490 
JCT03300 
JET03510 
JET03320 
JCT03330 
JE703340 
JEr035S4 
JET03340 
JET03570 
JCT03S80 
JET03S90 
JET03400 
JET03410 
JCT03420 
JET03430 
JET03440 
JET03650 
JCT03440 
JET03470 
JCT03480 
JE703690 
JCT03700 
JET03710 
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^05 ALAH«t (2.-SQRT(2. ) )/2. >•»<! ./ALAW) 
i;DCt-»CLBt/(OELBLf(|»UH-DEL!iL)/ALA«) 

IF (R.LC.4.) FVIS^l.O 

IF (R.0T.4.) FVIS^l )»^.24 

FH0M*l .-FCTA 

»ELFM^PFLPS»FtTA 

TAEF^^l. 

7HFFF»7N 

CALL VTUALL (R »|iUH»FVlS*F«OH»DELPH»TN» TLAM* AHa.KDEL f M» AUFM* 
lVH«*7HrTNEFFf TAFfF) 

VHO>VnN/VG 

OELFJ' . >/<NFR-l . ) 

RC0N»(?.-SQRT(2. > )/2. 

RAt>^KliEL7<l.-KtiEL)*RCQNt*< I ,/ALFtf) 

RU«EiWH/RAT 

r>CL7A-KI»fL4l»y 

T«»«Trt*TN 

f^UHT«7LAH(l>UH 

t>EL7HH^<THR-i . )/( 7N-1. > 

tHSTAG»THIf*VHNfVHH$(MTR'TN) 

DTSTAG-< THSTAG-l . )/(NTR-l.) 

AOl MRITE (6t9S) R i V^HG f VHNf E^UM » RUt I«EtTA ' DELF J»KDEL f TH» ms 
1 f i^UHTttiELTHN 
F-RfDR 

95 F0RrtA7 <lX»5<F8,4t2X)»8(F7.4»lX)) 

400 CONTINUE 

UPUASH GROUND NAXlNUN PRESSURE DISTRIDUTION 
IIRXTE <6«24) 

24 FORMAT ( / / //30 < , * t * • 4 t TUO-JET INTERACTION PROBLEM 


JET03730 
JET03740 
JET03750 
JCT03760 
JET03770 
JET03780 
JET037V0 
JET03800 
JET0T810 
JCT03820 
JET0383C 
JET03840 
JET03BS0 
Jer03860 
JET03870 
JCT03880 
JETO309O 
JET03700 
JET03910 
JET03920 
JET03930 
JET03940 
JET03950 
JET«3940 
JET03970 
JET03980 
JET03990 
JET04000 
JET04010 
JET04020 
ttIttV JET04030 


C 

C 

c 


* //) 

UALL JET REC*: - JL7AI0H EFFECTS 


IRAY-5 
AREC=1 . 

OnEG»( TGDS-1 . (RG/SD>«9AREC> 

IF (RO.GT.SD) OHEG-'TGPS-l . 

PRAY =IRAr-l 

PPHI"P1/DRAY 

PMI«0, 

WRITE (A«7t0) 

710 FORHAT (//30Xr'*t# RECIRCULATION EFFECT OH HALL JETS ••$'///) 
DO 750 XW«IflRAY 

PHID-PH1M80./PI 
WRITE (4t711) PHID 

711 FORMAT </50Xf 'PHI* ' fF8. 4# ' DEOS-V) 

CPMI2 COS( .59PH1 ) 

TAEFF*1 .F0HEG*CPHI2fCPMI2 
TNEFF*TH*(TAEFF-l . >»( TOPS* TNI/ ( f08S*l . I 
WRITE Ui 740) TNEFF»TAErF 

740 FORMAT ( /40X * ' TNEFF« ' *FlO , P . 2X.' ' TAEFF* ' .FlO .5/ ) 

IF <PH!.0E..5tPI) GO TO 701 
XH*SMTAN(PHI ) 

RU«SQRT(XU««2FSD9«3) 

IF (RW.6T .2.9SD) RWO.fSD 
00 TO 702 

701 RW*2.*SD 

702 NHPT*25 
XMUPT«NWPT-1 
PR«RW/XNWPT 
RrO. 

WRITE <«f7t.T> 

713 format (/37Xf *R' t 7X* 'VM/VN' i7Xt 'Th/TN' # 7X » ' TM/TA' i 7X» 'DPM/DPM' r 
I 'OTM/PTN'/) 

DO 703 IR^WNWPT 

CALL WALLJT <R*UMN*DHHf DELPW.TM»TMDiPM8f DCLTHNrDTSTAG) 

WRITE tAf704) R.UMN,TM,TMDfDELPH,DTS7AG 
704 FORMAT ( 30X f A < F 1 0 . 5 r 2* > ) 

K> RFDR 

703 CONTINUE 
rMI»FHHDPHI 

750 CONTINUE 

WRITE <At4AC) 


JET04040 
JET04059 
JE7O40A0 
JET0407G 
JET04080 
JET04C90 
JET0410C 
JET04U0 
JET04I20 
JET04130 
JET04I40 
JET041S0 
JET0416O 
JET04170 
JCT04180 
JET04190 
JET04200 
J 704210 
JLi04220 
JET04230 
JET04240 
JET04250 
JET042A0 
JET04270 
JET04280 
JET04290 
JCT04300 
JCT043IO 
JET04320 
JETC4330 
JET04340 
JCT043‘J0 
JET043AO 
5X. JET04370 
JET04^BQ 
JET04390 
JET04400 
JET04410 
JET04420 
JET04430 
J-T04440 
JET04450 
JET04460 
JET04470 
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460 FORHAT (//3Xf HhX;hu« GROUND FRESSURE AND TEhPEftATURE ALONG UPUASH.ET04 480 
IGROJHD STAGNATION LINE'//) j£T044P0 

WRITE (6.470) JET04300 

470 FORhAT ( 14X , XV/RN' .6X .‘ 'XU/S ' ,6X . ' DPHAX/DPJ' . 3X , ' F HAX /f NAXO ' . 3X t ' D JET04S 10 
ITHAX/DTJ' .3 x. ' TMAX/TNaXO' / > JCT04330 

FHI^O. JET04530 

FHIH*80.fPI/ie0. JETC4340 

DFHl =PMI‘1/24 . JET04550 

DO 450 Ii*S.25 JETD4560 

XU^SD4TAN(PHt > JET04570 

RM^S0RT(XU«$24SD««2) JET04580 

XUS=XM/(2.tSC* JET04590 

CPHI2-CC8< .5»PMU JET04600 

T«FF-l.40HE0»CPHI2»CPMI2 JETO^ilO 

TMEFF*TM4(TAEFF-1 . )t(TaPS-TN)/<TOIS-l . ) JCT04620 

CALL WALLJT ( Ry • UHN - »WH f DELPb t TN > TH» . PKB a*ELTNN t DT9TA0) JETC4630 

PnAX0«(PNB*l* )/(HPR'l. }fPNB«VHNVVHN«COS<PHl)$COS<PHI)/Th JCT04640 

TSO=*THB+(NfR-TH)tVHN»VHN#COS(PHI >$COS(PKI) 4ET04650 

^«AXG^(TSG-l . >/(NTT>l.) JCT04660 

ir (IL.EQ.l) DELPWO^PNAXG JET04670 

IF (IL.EO.I) r^*.TtfO-THAXG JCT04680 

PHAXHD-PHAXG/D..LF80 JET04690 

TKAXMO: Tf^AXQ/DELTUO JET04700 

WRITE (6.4SS) XU.XUS.PNAXG.PHAXND.THAXQ.TNAXND JET04710 

4SS FORHAT ( lOX » B ( F 10 . 5 . 2X ) ) JCT04720 

PMI*PHI4DPHI JET04730 

430 CONTINUE JE704740 

AC0N«1 .*(R0/S0>*T2 JET047S0 

PHI0D«0.^ JCT04760 

IF (ACON.UE.O. ) PHIC*0.0 JET04770 

IF (ACON.GE.O.) GO TO 667 j£T04780 

^'HIO-ARCOSC (SP/R0)*t.20) JET04790 

PHIOD=.PHIOtiaO./PI JET04BOO 

AC0K«2./(S1M .5tPI“PHIC)tt2) JE’ 04810 

XHOHZ»ACOF«( 1 .O4SlH<PHIO)t(SIN(PHI0)-2,O) )/3. JET048'0 

GO TO 668 JET04II30 

667 XNOHZ«ACeN42.t(l .'AC0N>/3. JCT04640 

668 WRITE (6»3j3> JETC48S0 

553 FORHAT < //5X » ' UPWAf ; HCHENTUH FUNCTION; JET RAOlUSf INTERACTION JETO406O 

ICONSTANTf VERTICAL UPPASH 'iHEET HOHENTUH. AND COALESCENCE ANGLE V ) JE704870 


WRITE <6. 551) ROiACQNfXHOHZ.PHIOD 
551 FORHAT < i5Xr 'RO/RH*' .Fi0.5i2X. 'ACON-' »F 10 . 3 t 2X » ' XHOHZ* ' tF10.St2Xi 
IPMI0(0EGS)«' rFlO.5/) 

CPMI?*!. 

TAEFF=1,40HEG«CPHI2*CPHI2 
T.4EFF»TNf(TAEFF-l . )»(T0eS-TN)/(TG9S-l .) 

CALL WALLJT<5D.VHH»DWH.DEL^b.TH>TKR.PHB»DELTHN.DTSTAG) 

CALL P3U>14T(CJ.DELpS.PELPWrRQ.SDfALP0f ALFUG) 

CALL SIHWiKDL..Fr . NFDt ALPUFD . RATU . CSV2W . CF U ) 

CALL S1HUC( ALPUG.RATUG.CSPliG) 

CU*2. tAW3FCSV2W/(RATU»C£PUtf ) 

WRITE (6*557) CJ 

357 FORHAT (25X.'UFwASH WIDTH CONSTANT ESTIHATE 
READ '5.66n IPDAR 
661 FORHAT (ID 


CU*' .FiO.5) 


TnST FOR IMTERACTiJH OF OELFLECTTOH REGIONS 


IF CUfSDrRO.LE.SD) 


TO 222 


c confutation of deflection KEGION interactjou 
C 

CALL lNTERQ(XHOHZ.PELPS*SD.RO.ALPCi£PS*CUf SIO) 


GO TO 223 
222 EFS^O.O 
nIOM.O 
FMIUO=0. 

PHI0»0.0 

WRITE (6.28) EP« 3IG.FHIU0.FHI0 

28 FORMAT (V20X-'JET AND UFu;,SH DEFLECTION REGIONS PC NOT INTERACT 
I CU ESTIMATE IS CORRECT'/ 

I /35X . ':,FS^ .Fio. 5. 2Xf SIG- ' » F 1 0 . 5 . 2X . ' TH I UO ^ . F 1 0 . 5 » 2X . ' FHI 0* ' 

I . 5> 

WRITE (d.:-> AlPUG 


JET048S0 
JET04890 
JCT04900 
JET04910 
JET04920 
JET04930 
JET04940 
JET049S0 
JET04960 
JCT04970 
JET04980 
JET049?0 
JETC3000 
JET05010 
JET05020 
JET05030 
JET05040 
JETC5030 
jer03060 
JCT05070 
JETOSOaO 
JET05090 
JET05100 
JET05UO 
JCT05120 
JET03130 
JET05140 
JET03130 
JET05160 
JET05170 
JET051S0 
JET05190 
FI0JET03200 
JET05210 
JET05220 
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^ PCoa 


29 

FORHAT <30Y. 'FMIN*0.0 ALONG UFUftSH tlNf. ALF »’G ^ » F8 . 4 . 1 X » ' IH 

UPWJET05230 

1 

lASH riCFLECTlON REGION'/) 

JET05240 


IF tIPliAR.NE.O) CALL GfLOT { AL PO f DELF S » S 0 . RO . EPS ; CU » S 1 0 r PHI UO » 

PHIQJET05250 

1 

1 ' 

JET05260 



JET05270 



JET052B0 

223 

URITE (6*111) 

JET05290 

III 

FORHAT t//30X» COMPUTATION OF UPUASH FLOW FIELD 

»tBBJET03300 


1 • ' / / ) 

JET0S310 


WRITE (6*11?) 

JET05320 

U2 

FORMAT (//35X» 'COMPUTATION OF UPUASH STREAMLINE PROPERTIES'/) 

JET0S330 


4PI/180. 

JET03340 


r MI»0. 

JET05330 

* • t 

FORMAT <//30X* 'PHlIi=' .FB.4» ' DECS'//) 

JET05360 


DO 600 l=l»? 

JET0S370 


PHlD = PMIMeO./PI 

JCTO330O 


WRITE (6*114) PHID 

JET05390 


ZEND^ZPLANE/COSiPHI) 

JET03400 


Y=SPtTAH(PHI) 

JET05410 


kU-SQRT'SD$t?FY9t?) 

JXT03420 


IF (RU.LE.RO) DUH=>BUHO(DELDPH«<RU-RGH> 

JET0S430 


IF (RU.GT.RO) DWN'AU3«RU 

JCT05440 


CPHI?-COS( .54PhI) 

JET05450 


TAEFF 1 .^OMEGtCPHI?tCFHi:; 

JET0S46O 


TMEFF=Tn+( TAEFF-1 . )«<TGBS-TMv/ (TGDS-1 . ) 

JET05470 


CALL UALLJT(RU»UHN»DUHG*DELPySiTM* TH8»FH8>DELTHH*DELTUS) 

JET03480 


IUR*1 

JET05490 


CALL UPWAST tROfRU*PHI^ZEHD,PMHG»D£LPS.DELruS»DELTUS*VUN*DELPUJ 

JET0S50O 


1 * TM*TM8*DELTMN,DSTAG*DUH*20bS» lUP) 

JET05510 


pmi^phhdfhi 

JET05520 

600 

CONTINUE 

JET05530 


HURT -25 

JET0554C 


XNUPTsNUFT-l 

JET05554 


A| PU=1 ,5 

JET0S56O 



JET 05570 

CIRCULAR CYLINDER DRAG FUNCTION 

JCT0SS8O 



JET05590 


DRAQsSlN(PHlSeP)4( t .*(4./l. )t<SlN<FHlS£P)«t2)) 

J£T0567e 


DO 700 J=l*23 

JET 056 10 


SUMF-0.0 

JET05620 


SUMFC^O.O 

JCT0S63O 


X8«0, 

JET0S64O 


DX8*XL2/XNIiPT 

JETO56S0 


WRITE (6*116) 

JCT05460 

U6 

FORMAT (/////30X*'UPWASH PROPERTIES COMPUTED AT I LOCATION OF 

UNDeJCT05670 


1R8IDE OF iOPY'//) 

JET0S660 


WRITE (6*118) ZPLANE 

JCTe5690 

lie 

FORMAT (/35X*'80DY CENTERLINE LOCATION AT Z/RH- ' * F 10 . 5// ) 

JCT09700 


WRITE (6*78) 

JET03710 

78 

FORMAT {//IX* '2U/RN' t5X» 'RU/RN' *5Xi 'PHI8' iSX* 'X/RM' *5X# 'VMU/VN 

'» JET05720 


1 3X* 'DPU/PN' f 4Xf '8UH/RN' »2X r ' DPU/PUO' *1X * 'DP'J/PJS' *3S* 'TM/TA' *3X* 'CJCT05730 


2TSM/DTH' *1X, 'P8S/PN'/) 

JET03740 


IWR»0 

JET05730 


DO 500 r 1*NUPT 

JET05760 


FACT*2.0 

JCTC5770 


IF (I.EQ.l) FACT*i,0 

JETQ3780 


IF ( 1 EO.NUPT > FACT^l , 

JET0S79Q 


Z8*Z80DY(X8> 

JET0S800 


PH18xATAN(XB/(SD4ZPLANE-ZD>) 

JCT0S810 


XW-SDiTANCPHli) 

JET0582C 


RM«SORT(SD««2FXW»«2) 

JET0583C 


CPHI2*=C0S( .StPHIp) 

JET03E40 


TA£FF*1.+OHEO*CPHI2ICPHI2 

JCT05850 


TNEfF»7M4(TAEFF-l , )»(TG8S-TM)/(TGIiS-l . ) 

JET0S840 


CALL WALUJT<RW*UMNr8UHG«D£LPUS»Tn*TMBrPMB>DELTMNiDELTWS) 

JET0Sa?0 


ZU»2PLANE-ZB 

JC70S880 


ZEND«ZU/C0S(PH1P) 

JET05890 


CALL UPWAST ( RO* RUr B * ZEND * 8UH0 * DELPS * DCLPWS * DELTWS * VUN • 

JET05900 


lDELPUJ»T)1*TMB*&ELTNNf DTSTAOtBUH.ZOUSr lUR) 

JET03910 


THEr#>^ATAN(DZDXB(XB) > 

JET05920 


ETAWB»W2( XB)/ ( 3.66BUH) 

JET05930 


IF (ETAWB.6T. 1 . ) ETAW8-1 . 

JET05940 


FB'FETAWBCETAUBfALFU) 

JET05930 


DCYL-^W2(XB)»0RAG 

JET0S94O 


DFLAT*3.6tBUHBFI 

JET05970 


SUMF*SUMF+ .54DXBtFACTf DFLAT«(UUN»*2)*(C0S(PHIBFTH£TB)tt2). TM 

JET03980 


Fignre 69 - Cominueri 


14 




SUHFC«3UNFC^.S>DX»«FACT«DCYt«(VUN$*2»«tC0S<PHIfifTHCTB)t»2)/TM 

ZOU«ZOUStCOS(FHlB) 

PHID«FHIB»ieO./FI 

IF (ABS(FHIB) .LT. 1 .£-4) D£LPH«D£t.PUJ 
IF <1»£LPH.M£.0. ) 0ELPB^D£LFUJ/D£LPH 
IF (P£LPH.£Q.O. ) D£LFB«I*£tPUJ 
l•£LPUS*P£LPUJ/^£LPS 
P|iS«< (VUN«CaS<PHlB4TN£TB> >>:«2)/TH 

• VfkITE ZUtRM»PHIO»X>»VUN»D£LPUJ>BUHf D£LPafI£LPUS*ThB>DTSTftOf 

IPBS 

77 FQ^^HAT UXt 1 1 (F8 .4i |^) rF6 .1) 

XB*Xt4DXB 
500 CONTINUE 
C 

C FLAT BOTTOM LIFT FORCF 
C 

XLUT»< 1 ./PI )#SUNF 
UNITE (Af778) 

778 FORMAT (//30X»'FLAT BCTTOM VEHICLE UlTH SHARP CORNERS'/) 

UNITE <4>771) XLUr 

771 FORMAT (35X»'UPUASH LIFT FORCE ' *5X . 'LU/2TJ* ' t FI0.5//) 

C 

C UPUASH LIFT FORCE FOR CIRCULAR CROSS SECTION 
C 

XLUT-f l./PI )*SUMFC 
UNITE «4»773) 

773 FORMAT <//30X#'BODY UlTH CIRCULAR CROSS SECTION'/) 

UNITE (At77l) XLUT 
2PLANE-ZPLANE42.tOZPL 
IF (ZPLANE.6T.7. fZFINAL) STOP 
700 CONTINUE 
EHB 


SUBROUTINE VTSOLV < TNr BVH* ALPU* TLAH»Vn» TMtRVrRT ) 

NC0M»(2.-S0RT(2. ) )/2. 

C UNITE <6>200} TM. BVH*ACPV. ALPT 

200 FCNNAT < IX» ' TN« ' *Fl0.5#2X. 'BVH«' tF|0.3t2X» 'ALFV=' t Fi0.5»2X* ' ALPT-' 
ltFlO.3) 

BTN-TN-l, 

REL*-3 

BELMlN-l.C-5 
RU«<RCGN)4t< I . .'ALPV> 

NTcRV 
ALPT-P.LPV 
VH^l . 

TH*1. 

BO 100 1-1*300 

CALL StMTVlTNfTMrALPU*TLAMtCH*CT*CNTN*CTTN> 

C UNITE i0*203> CHtCTfCMTMrCTTM 

203 format iiX» •CM-'#F10.3*2X* 'CT-' * Fl0.5*2X* CMTM* ' »F10.3* 2X* 

! 'CTTM-' .F10.3) 

F=RV«NVtTH-2.tVM7v;i55yHtByMBCM 
6*DTN4TMtNV$RV-2.<yMU<TM>TN>I . )tfVHtBVHtCT 
FTM-TM-2.^VMtVM9BVM«ByhtCNTM 
) VM»- 4 .«VMfByH* 8 'i:f 5 Crt 

GTM »PTNfRV$NV-2. •VM4< TMtTN-1 . > t BVH«»VHtCTTM- 2 . FUMt BUHt Ry>!»CT«TN 
fiVM^-2. t'TMtTN-l . ) BVHtBVHtCT 

bET*FVMt6TH-6VMTFTP 

RE*.yM.(GtFtM-FtGTM)/OET 

DCLTH-(FtGVH>6fFyM>/9iT 

VM-UMfRELtDELUM 

TM-TM4NCL4PELli; 

C UNITE 40-505) Vn«DEL^Mt?N>tELTM*F*6 

303 format UX- UM- ' -FIO.S-IX- BELW- ' -F 10 . 5* 1 X * ' TM« ' *F 1 0.S . 2X * ' DTM* ' ,F 
ll0.5.IX.'F^**FI0.3*lXf'8*'tFl0.3) 

IF 4ABSCF) .LT.rtELHIN.AND.ABSCC) .LT.DELMIN) 00 TO 101 
100 CONTINUE 
STOP 

to I RETURN 
END 


Figure 69 - Continued 


148 



ORIGINAL PAGE IS 
OF POOR O’^AUTY 


SU»ROUTINe UALLJT <R>VNN»»UH»PCLPJ«rNf Tn»»PH)»DCLTHNtMSTAO) 

REAL KDELO»HO*ADELFD>NFD>Ni ADCLfNPRrNTRtHACH 

COKHO'I/UALL / DCLS • NO P KDCLO p ALP 6 1 ALPNO p ALPWFD p RO t ROHp VO p »UH0 

CONNON/UAIL r/TLANf p DELPS* OELTSi «F'St AR 

COHHON /NQZZUE/ HPR p NTR p THt OAH p HACH t ACl'IS p TaCFF p TNEFF 

RC0H«<3. -SQRTC:. ) )/2. 

NFD*7*0 

KPEtFI>«l./f. 

»U3HO>^Ay3tRO 

AU:'*.09f (l.fALGG<2./( 1 .FI ./TNtt.S) ) ) • ( I . - . 1 6«ACHS) 
r J!'3«AU3tRQ 

•ELRRH* < BU3H0-pyHQ ) / ( RQ-RGH) 

PELHR>4HFp-NO)/(RO>RGH> 

DEL AR* ( ALFUFD'ALPVO) / ( RO-RGH > 

ADELBT«.0175 

ADEtBL*(ADELBTt(RQ-2. >-D£LS >/(RO>ROH> 

TLANFD»1 .F.185$TN«t.5 

TLANSL»(TLA«FII~TLAHI)/<RO-RCH) 

iF (R.0E.R6H) TLAN-TLAHIfTLAHSL«(R-ROH) 

IF (R.LT.ROH) TLAn«TLAIIl 

IF (R.LE.RO) FCTA«( I .-<R/R0>t*ALPG)ft4 

IF (R.OT.ROl FETA»0.0 

IF <R.OT.RO> TLAh«TLAHFO 

F«B»1 .F<NFR~1 . )$DELPS«FETa 

IF (R.GT.RGHI GO TO 402 

Ttt»«l.iD£L7S» 

TH«THB/NTR 

ACHI*SQRT(2.t(PSBAR/FHB-l . >/ qaM> 

V«H» < ACHI /HACH) tSORT ( TMB/MTR ) 

OELPJ«PELPS 

BELTA^KPELOtByO 

N=>14 

khel^kdelo 

BU*BUO 

BUH«RATUOfBUO 
VltG>VHN/VO 
BUHT>TLAHI4BUH 
DELTHN«(THB-1. )/(HTR-l.) 

TNSTAC>TNB 

r-TSTAG«(THSTAQ-l . >/INTR-l .) 

60 TO 401 

402 IF IR.OT.ROl 00 TO 403 
iUH«BWHOfDCLBRHt<R-RON> 

ALPy-ALPW04DELARt(R-RGH) 

BELBL«DELSFABELBL«CR>RGM) 

N«NOtPELNRf (R-RGK) 

GO TO 403 

403 BUH«AV3fR 
DELBL«ABFlBT4(R-3.) , 

N«NFO 

ALPU-ALPyFD 

405 ALAN* C (2. -SORT ( 2. >)/2. )••<!. /ALPy) 
ROEL»P£LBL/(BELBL*(ByH-DELBU/ALAN) 

IF (R.LC.4.) rviS«1.0 
IF (R.GT.4.) FVlS»l./<R/4.)t4.24 
FMOH-1 .-FETA 
DEtPH»D£LPS*FETA 

CALL VTWALL (RpBUHfFUlSpFnON.DELPHpTHpTLANtPNBpKBELpNpALPyp 
iVriNpTNpTNEFFpTAEFF) 

V«C*VHM/VO 

RC0M«O.-S0RT<2. ) )/2. 

RAT *KDELFU . -KD£L )»RC0N4* U . /ALPN ) 

BW»ByH/RAT 

B£LTA«KBCLtBy 

0£LPJ*(PNB-I . )/<WPR-i . >fPrtB4VHIf*VHH/TH 

THB-TN9NTR 

ByNT-TLAHtiMH 

B£LTNH*( TNB-1 , )/(NTR-l . > 

THSTAG'TNB 

BTSTA6*<TNSTA6-1 . )/<NTR-l . ) 

C 401 VRITE <4p95> R pVHGpVHKp BUH- BV » DELTA * OELPJ pKDEL p TH p THB 
C IpBllHTf DELTHMpDTSTAO 
401 CONTINUE 
R*RFPR 

93 FOPNAT ( 1 X f3<F0. 4 p 2X > ,6?F7. 4f IX) ) 

RETURN 

END 


JETOA7AO 

JET04770 

JET0A760 

JET00790 

JET06800 

JETOA810 

JET0A820 

JET04830 

JET06840 

JET04630 

JET0A840 

JET04870 

JET04880 

4ET04890 

4ET04f00 

JET04910 

JET04B20 

JET04930 

JET04940 

JET04950 

4ET04940 

JET04970 

JET04980 

JET04990 

JET07000 

JET07010 

JET07020 

JET07030 

JET07040 

JET07050 

JET07040 

JET07070 

JET07080 

JET07090 

JET07200 

JET07HO 

JET07120 

JET07130 

JET07140 

JET07130 

JETO7140 

JET07170 

JET07180 

JET07190 

JET07200 

JET07210 

JET07220 

JET07230 

JE^07240 

JET07250 

JET07240 

JET0727O 

JET07280 

JET07290 

JET07300 

JET07310 

JET07320 

JET07330 

JET07340 

JET07330 

JET073A0 

JET07370 

JET07380 

JET07390 

JET07400 

JCT07410 

JET07420 

JET07430 

JCT07440 

JET07430 

JET07440 

JET07470 

JET07480 

JCT07490 

JET07500 

JET07510 

JET07520 
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JET07930 

JET07540 

JET07550 

JET07360 

SU9R0UTIMC VTUALL <R i RUH i TVIS f FHCH* DPM* TMt TLAH » RH»rKDCL » Ht ALf‘V» JCT07370 
iTHfTNEFFf TAEFF) JET07380 

REAL KDEL>N JET073^0 

VHsi, JET07AOQ 

TN«i. JET07A10 

DELH1N*1.E-3 JCT07A20 

REL-.3 JET07430 

DO too I«l*300 JCT07440 

TK»«TH$TN JCT07430 

CALL SIHUTV (rNtTH*TLARrEf1ifK0CL»N»ALPV«CV2»CPtCT*CV3TH>CTTH» JCT07440 

IRV) JCT07470 

FnRV«FH0N«FVIS-Rt»UH«(2.4PNi>VH$VHtCV2/THfDPHtCP> JET07460 

8aRM$FH0nt(TNEFF<TAEFF)tTN-2.tPHBtv*H«<TH»-TAEFF)«R4tUHtCT JCT07490 

FVH«-4.tRttyH4(P?lD$VnfCV2/Tfl) jeT07700 

FTH*2.iRfBI»HtPHDtVH«VH8<CV2/(TNtTM>-CV2TN/TH) JCT07710 

Ovn«-2 . tPH»« < TPD-TAEFF )4R»DyH4CT JCT07720 

6TH«RW$FH0Ht<THEFF-TA£FF)-2.tPRD4VM$RttyMt(CTtTM4<THi*TAEFF)»CTT«J JET07730 
DET*FVHtGTH-OV««FTN JET07740 

DCLVH>=(0tFTN'-rt6TH)/DET JCT07730 

DEL TH« ( F f 3Vn-GtFvn ) /DCT JET07740 

VN'^VH4RCL«DELVN JCT07770 

Tn«THfRELfDELTH JET07780 

) WRITE (4>70S) FtG*VHfTH JET07790 

703 FORMAT CIX * * F» ' * F 10 . 3 » 2X » ' G= ' *FI0.5» IX# 'WN» ’ .FlO.Sf IX» 'TM»' rFtO.S) JET07800 


'ELMIN.AND.ADS<G> .LT.DELNIN) GO TO lOt 


IF (ADS<F> 

100 CONTINUE 
STOP 

101 continue 

RETURN 

END 


SUiROUTINE UPUAST (ROr RWALL . PHI # ZEND* DUGH# DELPS» DELPyStDELTUSt 
1VNUN*DELPU3»TH*THD*DELTHN*DTSTA0#R|;N#Z0U# IVR) 

REAL NPRtNTRrhOHU 
COMMON /JET/ SD«RSH 

COHMON/NOZZLE /NPR#NTRtTN# GAM* MACH >ACHS#TAEFF#TNEFF 
REAL MACH 

THIS ROUTINE COMPUTES THE UPUASH STREAMLINE PROPERTIES 
A0U«3« 

PI«3. 14139243 

ALP0P*1.3 

ALPUH-l.S 

ACON> 1 . > ( RGH/SO ) tt2 

IF lACON.OE.O.) GO TO 907 

PHI0«ARCQSUSD/ROK)t«»20) 

AC0N-2*/<SlN( .34Pl*PHI0)«t2) 

PH10D»PHI0S180./Pt 

IF (PHI.OE.PHIO) M0MU«AC0NtSIN(PHI-FHI0)9«2 
IF (PHI .LT.PHIO) n0MU»0. 

GO TO 909 

907 MONU«ACON42.t< 1 .-AC0N)iSlN(PHl)t»2 
909 CONTINUE 

773 FORMAT < 1 X « ' ACON» ' *F 10 . 3« 2X t ' HOMO- ' *F 10 *3> 

AU3-.30 

ALU1-.3 

ALU2-.0 

Z0U-AOU9DU6H 

VRITE (4*108) SD*RGH*ZOU 
108 FORMAT ( IX* ' JET ' * 3F10.3) 

ZV>ZOLtCOS(PHl) 

ALU*. 5 
ALPU-t .3 

CALL SINU( ALPUtRAT*CV2U) 
yRITE (6*111) CU2U 
111 FORMAT ( IX* 'SINU' tFlO.3) 


determine OEPATURE FROM INVISCID DEFLECTION 


JET07810 
JET07820 
JET07830 
JET07840 
JET07aS0 
JET07840 
JET07870 
JET07880 
JET07890 
JET 07900 
JET07910 
JET07920 
JET07930 
JET07940 
JET079S0 
JET07940 
JCT07970 
JET07980 
JET07990 
JET08000 
JET08010 
JET08020 
JET08030 
JET08040 
JETO80S0 
JET08040 
JET08070 
JET08080 
JET08090 
JET08100 
JET08110 
JET08120 
JET08130 
JET08I40 
JCTC8130 
JET08140 
JET0B170 
JET08I80 
JET06190 
JET08200 
JET08210 
JETO022O 
JCT08230 
JET08240 
JET08230 
JCT08244 
JET08270 
JET08280 
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DETERHINE IHITIAL UPUASH UIDTH 


DELPMJ-DELFUS 
ZETAG«»UGH/ZOU 
FETAGP»U .-ZEfAGt*ALPUP)»»< 

C WRITE (6tll7> ZETAGtFETAOP 

117 F3RHAT UX i ' ZETAG- ' iFlO . 5i 2X» 'FETACP* ' iFlO .3 ) 

FETACN«U .*ZETAG««ALPUH)tt4 
FROHU^HOHU«< 1 * -FETAGfO 
P«BI*1 .F(NPR-l. >*DELPWJfFETAGP 
DELPW»-(NPR-1 . liDELFUJ 
PSG«l.4tiEUPMP 
DELTCPXMTR-1 , )«PELTWS 
C WRITE (A»i07> DELPWRtDELTSP 

107 FORHAT U X . ' DTP« ' . F I 0 . 5 . 2X i ' OTT- ' . f I 0 . 5 ) 
ACHl-S0RT(3.#<PS»/PHRl-l. )/CAH) 

THBI-l.+PELTSR 

WHNUl-<ACHl/nACH)«SQRT(TMBI/NTR) 

THl-THil/MTR 

IF <IUR.NE.O> WRITE (A»33> WtINUI r THB ! » DELPU J » DELTUS * BUGH 


JET06290 

JET08300 

JET0B310 

JET08320 

JET08330 

JEY08340 

JETO033O 

JET083A0 

JET08370 

JET083G0 

JET08390 

JET08400 

JE708410 

JET08420 

JET08430 

JET08440 

JET08430 

JETQ84AQ 

JET08470 

JET08480 


53 FORMAT (IXi ' IHV . VEL . « ' # F 1 0 . 3 » 3X t ' TMI /TA« ' t FI 0 . 3 1 2X i ^BELPWJ-^ «F10.SJET08490 


l»2X. 'OELTWJ-' »Fl0.3f2X» ' BUGH/RM= ' iF 10 .3) 

RAY-RUALL4BUQH 

CALL SIHUM (TN»THI .TLAMItFHBIrALPUfFETAOPiDELPWJtVHHUI >BUONrRAY» 
IRAT»FM0MU»TMEFF»TAEFF) 

BUOHH^PATtBUON 

IF <ABS(nOHU).LT.l.E-S) BU0N«l.E-3 
IF (ABS(M0MU).LT.l,E-5) BUOHM-l.E-3 
IF (IWR.NE.O) WRITE (A»10) ZOUt BUONt BUOHN* TLAMI 
10 FORHAT ( |3X» 'ZOU/RN*' >F10.5>2X» 'BUO/RN*^ tFlO.3 
1 »2X» BUOHM/RN=' rFlO,3#2X» 'TLAHI*' »Fl0.5//> 

NUPT-25 

XNUPT*NJPT-1 

0Z»<2EMD-BWGH)/XNUPT 

Z*BWGH 

IF <1WR*NE.0) WRITE <4»t2> 

12 FORHAT </X* 'Z/RN' »7Xt 'WHU/VM' *7X# 'TH/TM' f7X 
IX» 'BO/RN' lAX* 'DPU/DPJ' »AX» DTH/DTW'/) 

TLAHFB«1 . F . l03tTHf $ . 3 
TLAHSL>(TLAHFD-TLAHI )/ ( ZOU-BWGH) 

TM»THI 
UHUH-WHHUI 
DO 100 I«1»NUPT 

BUH«BUOHNf AU3»< Z-BWGH>f (COSrPHi )«$ALUl > 

BU«BUH/RAT 
ZETA-Z/ZOU 

IF (Z.GE.ZOU) ZETA«1. 

IF (Z.OT.ZOU> TLAM-TLAHFD 
IF (Z.LE.ZOU) TLAH-TtAHlFTLAHSLt(Z-BWGH) 

FETAP-( I .-ZETA»#ALPUP)»»4 
FETAH«( l.-ZETAi«ALPUH) tt4 
RAY-RUALLfZ 
FnCH«HOMUt( 1 .-FETAN) 

DPH«0£'CPWJ»FETAP 
PHB«l.F<HPR-l.>tOPM 


JCT08S00 
JET03S10 
JET08S20 
JET08S30 
JETO0S4O 
JET08S30 
JET08SA0 
JET08370 
JCT08S80 
JET08390 
JET08400 
JET08410 
JET06620 
JET0BA30 
JET08A40 

TH/TA' »7Xf 'BUH/RH' »6 JET08A30 
JCT08660 
JET08A70 
JET08680 
JET08690 
JET08700 
JET08710 
JET08720 
JET08730 
JET08740 
JET087SO 
JET087A0 
JET08770 
JET08780 
JET00790 
JETOBBOO 
JET08910 
JET08820 
JET08830 


CALL WTUPW (RAYfBUHtFHOHiDPHf TNfTLAnfPHBrALPUiVNUNrTHiTNCFFiTAEFF)JET0884Q 


905 DELPllJ«(PMB-l . )/(HPR-l . ) FPMBfVHUHfL'HUH/TH 
THB»THtNTR 

DELT«H»( THB-1 . )/(NTR-l . ) 

*^*FSTA0*TI1» 

tfTSTAO«(THSTAO-l . )/(HTR-l . > 

IF (IWR.EO.O) DELPUJ»DFLFUJ/P£LPS 

IF (IWR*N£.0> WRITE (4*11) Z * OfiUN i TH > THB • BUH * BU* DELPU J * PEL THH 
11 FORHAT <2X»9<F10.3*2X)* 

Z»ZFPZ 

100 COMTIHUE 
RETURN 
END 


JET08830 

JCT08864 

JET08870 

JET08880 

JETO089O 

JETOB900 

JET0B910 

JET08920 

JET08930 

JET09940 

JET0B9S0 

JETO09oO 

JET08970 

JET08980 

JET08990 

JET09000 
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SUiROUTXNE VTU»>U (R i f UH . FHOH f DFHt TNi TL AM » PM6 » ALPV. UH * TM* TNEFF r JET09010 

ITAEFF) JET09020 

DELMlN»l.E-5 JET09030 

REL-.5 JET09040 

X»0 100 I'ltSOO JET090S0 

THi>TH»TN JET09040 

CALL SlMUTV (THf TMf TLAMfPMB»ALPV»CV/2iCF t CT »CV2TM • CTTM • RM ) JET09070 

F»RV$FM0f1-RtBUH»(2.«PMii«VH*VH$CV2/THiDFM$CP) JET09080 

0’*R^/fFMOM ' I . )tTM-2*»FMB»VM«<TMB-l . )*RtBUH»CT JET09090 

C URITE <4»lli) FfQfTM.VM JET09100 

111 FORMAT <IX» 'F^' »Fl0.Si2X,'G*'.Fl0.5»2Xi 'TM*' »F10.Si2X, VM«'#Fl0.i)JET09ll0 
FVH^-4 . «RtBUHt (PHB$VM»CV2>'TM> JET09120 

FTM«2. tR*BUHIPMBfVH*VH$<CV2/(TM«TM)-CV2TM/TM) JFT09UO 


C 


c 

c 

c 

c 

c 

c 


c 

c 

c 


GVM*'2. tPMBf( T hB -1 , )«R*BOH*Cr JET09I40 

GTM^RVtFMOMt(Th-l, )-2.*FMB»VM*RBBUM*<CTtTMf(THB-l . )tCTTM) JET09130 

DET^FVMtOTM'OVMtFTM JET09U0 

BELMM-(3«FTM-FtGTM)/DET JET09170 

DELTM-(FtCVM-OtFVM)/MT JCT09180 

VM*OM+RCL«DCLVh JCT09190 

Tn-TM+RELtDELTM JET09200 

URITE <A»70S) FiGiUH»TH JET09310 

705 FORMAT (XXr ' F» ' #FlO . 3 »2Xt G- ' rFlO .5 . IX r 'Mfl« ' »F10 . 5f IX * ' TM« ' tFlO . 3 ) JCT09320 

IF (ABS(F) .LT.DELMIN.ANO.ABS(O).LT.DELMIN) 60 TO 101 JCr09230 

100 CONTINUE JET09240 

STOP .ICT09230 

101 CONTINUE JET09260 

RETURN JET09270 

end JET09280 


JET09290 


JFT09300 

JET09310 

JCT09320 

subroutine IHTCRO (XLAMfBELPSiSOtROf ALPOfEPStCUrSIO) JCT09330 

JET09340 

MAIN routine for JOHPUTIHG GROUND PRESSURE DISTRIBUTION JET09330 

THIS ROUTINE FINDS THE SOLUTION FOR THE UPUASH THICKNESS AND JET093AO 

JET IMPINGEMENT PERTURBATION PARAMETCRl SIGMA AND EPS BY JET09370 

MATCHING PRESSURE INTEGRALS JET09380 

JET09390 

COMMON /HEIGHT/ HDiIPBAR JET09400 


DIMENSION EP(200)»C<200)iPHIU{200) »PHK2C0) »CSJ< 20C ) tFUI < 4 ) »F JI< 4 > JETC<>4 10 
CSPJ(ETA»A)«.5*ETA*»2-(4./(AT2.))tETA*«(A' )*(3./(A*l . ) ) tETAt* ( 2 . JET09420 

ltA+2. )'(4,/(3.»AT2. ) U£TA««C3.*AT2. )♦< i4.tAT2. > ) «ET A«$ < 4 . BA42 . > JCT09430 


FETA<ETArA)><l .-£TA«$A>«S4 

EP<1)*.0 

URITE (6»777) 


JET09440 

JET09450 

JET094A0 


777 FORMAT <///|X»'*tt CALL INTERG «•»'//) JET39470 

IFLAG-0 JET09480 

C<1)-1. JET09490 

PI»3. 14159243 JCT09300 

SIC1»1 .0 JCT09310 

CUl-CU JET09320 

ALU-l .0 JET09330 

PARH*.01 JET09340 

DO 900 J«l>100 JET09S30 

EPI«EP(J) JET09340 

SIOI-C( J) JET09370 

DO 1000 XTR*1.4 JCT09S80 

IF UTR.EQ.2) EPI»EP(J)4PARM JET09S90 

IF (ITR.EQ.3) EPI«EP(J) JET09400 

IF nTR.EO.4) S101-C( J)4PARM JET09410 

JET09420 

SOLUE FOR INTERSECTION OF PRESSURE BOUNDARIES JETC9430 


JET09440 

DO 100 1»1»99 JCT094S0 

IF (I.EO.l) PHIU<1)«1.42 JET09400 

'"J>SOtTAN(PMIU<in JCT09470 

XUP«SD/(COS(PHIU(I) )t$2> JET09480 

TU»CUItSX01tSD/(C0S(PHtU( I > )ttALU) JET09690 

YUP»'CUIiSlG!4SD«ALU/(C08(PMIU< I ) >t«(ALU41 . > ) JET09700 

RUSa-XU»«24(SD-YU>4*2 JET09710 

RU-SORT(RUSG) JET09720 

ROP-RO+EPIt(l .♦( (5D-YU1/RU) ) JET09730 

R0PS0*R0P«I2 JET09740 

F-RUS0>R0PS0 JCT09730 

FP«2.*xu*xup-2.t<9P“YU)trap-2.frpr4ROP*rso-ro)t jcto97ao 

1 < YUP/RU*(2.4XU*XUF-2'.t<SD*VU>tYUP>/RUS0) JET09770 

PH1U( IM ) *PHIU« 1 ) -F/FP JET09780 

IF (ABS<PHIU(Hl)-PMIUU)).LT.l.E-5) 00 TO lOl JET09790 
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100 CONTINUE 
STOP 

101 PHIUO«PNlU<Zfn 
PHlUOD»PHIUO«ieO./PI 
YU0*CUI»SICIfSD/(C0S(PHIU0)$9ALU> 

PHI0-ATAN3<S0«TAN(PH1U0) iSO-YUO) 

PHl0D«PH10tie0./PI 

IF UFLAG.EQ.i) URITC <A»10) EPl »CUI i PH1U0D»PHI0D 
10 FORMAT (13X» 'EPI=' fFl0.5»2X» 'CUI«* >F10.3»2X* 'PHIUO*' fFtO.3 
l»2Xf PHIO«' *Fl0.5*3Xt 'PEGREES' ) 

DETERHINE HOriENTUH INTEGRAL FOR IHPINGEHENT REGION 

NUPTS-99 

IF <J.GT.20> NUPTS^198 
XPTS*NUPTS-1 . 

DFHIU^PHIUO/XPTS 
PHIU<1 >^0. 

SUHU1»0. 

00 200 IJ«l#NUPTS 
FAC1»2.0 

IF (IJ.EQ.l) FACT*1.0 
IF (IJ.ED.NUPTS) FACT*l.O 
RU*SO/COS(PHIU(IJ)> 

CALL UALLJ (RU *UMNtDUH2>DUH3»DUN4»0ELPU) 

IF < IJ.EO.l. AND. IFLA6. CO. 1) WRITE (6tlll> PELPW 
111 FORMAT (33X* 'PHI*0.0» DELPUO« ' »FlO.S) 
rU-CUX»SIGl«SD/COS<FHlUUJ) ) 

XU«SPtTAN<PHIU< IJ)> 

PHI(IJ?^ATAN2(XU*SD-YU) 

RU«<Srt-YU)/COS<FHX < IJ) > 

R0P«R0FCPI«(1 .FCOStPHI (IJ) > ) 

ETAN*^RU/ROP 

IF lETAN.GT.l.) FETAGsO.O 

IF (Etan.gt.1.) go to ton 
FETAG^FETA(ETAMf ALPG) 
toil PNlN^DELPStFETAS 

IF < IJ.CQ« i >AND. IFLAC.EQ. 1 ) WRITE (6*773) PHIN 
775 FORMAT ( 33X » 'PHI -0 .0* F MIN* ' »F 10.3> 

PHAX*PELPM-VHMtA2F(VNHtC0S<PHIU< IJ) ) )f *2 
CALL F HATCH ( PMZW# PNAX r ALPUO) 

IF <IJ.EQ.l.ANP.IFLAG.EO*n WRITE (6tlOSS) ALPU6 
1055 rORMAT (33Xf FHI-O.Of ALPUO- ' . F 1 0 . 5 ) 

CALL SinuOt ALPUO»RATU»CSFUO> 

P«AX-nELPU-VMN#f2f <\;HN*C0S<FHI0< IJ) ) )$t2 
FU=0ELPS*FETACFCSPU6»(FHAX-PELFS$FETAG) 
FUT'*5I61tFU/<C0S(FHIU( IJ) )**3> 

SUNL»l-*SUMJl F .5tOFHIU*FACTfFUT 
IF (ETAM.GT»t.) ETAM^^l.O 
CSJ( IJ) CSF J(ETAN* ALFG) 

PHlUdJF^ =MIU( I J)FPPMIU 
200 CONTINUE 

XIlUG>CUltSUNUl«SPtf2 

SUH1»0. 

NUPTSH«NUPTS*1 
PO 300 1J=«1»HUFTSM 
ROPP*ROFEPm 1 .FCOS(PHIUJFl))) 

R0FFSQ*R0PPtt2 

ROPM«ROFEPI«( 1 . FC0S(PHI ( I J) ) ) 

ROPftSO^POPMM2 
Fn>ROPHSQ«CSJ( IJ) 

FP»ROPPSQ«CSJ( IJFI ) 

PPMI*PMI ( 1 JFl )-PHI ( IJ) 

SO«!*SUMl+*3*PPHIt(F«6FP) 

300 CONTINUE 

Xll JO«^PELPStSUMl 

CALL SIH(0 ,0«ALPG*R« . j*CSPJO) 

XI2JG«CSPJ6tDCLPS«aPI-PNI0>t<R0«(R0F2*>EPZ )F1 .S9EPltt2) 
1-2. tEPI«<ROFEPl)«SlN(PHIO>- .25iSlH(2.9PMlO)tEPl**2) 

XJl*xn JG+XI3J0 

FJi'( ITR)*XJI-PI 

PPHIU«( .5*F I-PH1U0)/XPT5 

SUMU2*0.0 

PH»PM1U0 

XI2U0»SIGIt< 1 .-SIM<PHIU0)> 

FU: ( ITR )»Xr 1UGFXI2UG-XLAM 

Figure 69 - Continued 


JET09G00 
JET09810 
JET09620 
JET09830 
JET09B40 
JET09850 
JCT098A0 
JET09870 
JET09880 
JET09890 
JET09900 
JET09910 
JET09920 
JET09930 
JET099AO 
JET09950 
JET09960 
JET09970 
JET09980 
JET09990 
JETIOOOO 
JETIOOIO 
JET10020 
JETI0030 
JET 10040 
JET 10050 
JET 10060 
JET10070 
jETiooao 

J£Tt009C 
JETlOlOO 
JETIOUO 
JET10120 
JET10130 
JET10140 
JET10150 
JETlOlOO 
JCT10170 
JET10180 
JET10190 
JET10200 
JET10210 
JET 10220 
JETI0230 
JET10240 
JET10250 
JET10240 
JET1027O 
JET 1 0280 
JET10290 
JET10300 
JCT10310 
j£Tt0320 
JET10330 
JET10340 
JET103SO 
JET10360 
JET10370 
JET10380 
JET10390 
JET10400 
JETI0410 
JET10420 
JET10430 
JET10440 
JET10450 
JET10460 
JET1047O 
JET1048 . 
JCT10490 
JET10300 
JET103IO 
JET10520 
JET10330 
JET10S40 
JCT10350 
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IF <IFLAO»EO.n 00 TO 907 JET10560 

IF (iTR.EQ.U GO TO 1000 JCT10S70 

RELAX*. S JET10500 

OEFI-FARH JCTI0390 

DCU1*FA9H JETI0600 

IF (1TR.E0.2) FUIE*(FUI< ITR)-FUI ( ITR-1))/DEP1 JETIOAIO 

IF 4TTR.EQ.A) FUIC* < FUK I TR ) -FUI < I TR- 1 ) ) /DCUI JET10A20 

IF <UR.EQ.2> FJ!E»<FJI(1TR)-FJI<ITR-J))/DEFI JET10A30 

IF (ITR.EQ.4) F JIC*<FJIUTR)-FJ1UTR-1))/PER1 JE110640 

1000 CONTINUE JETIOASO 

DEr*FUlF«FJTC>F JlEtFUIC JETlOA^O 

P£LE«(FJI< 1 )tFUIC-f UK 1 )iFJlC>/DEV JET10A70 

OELC«<FUK DtFJIE'FJK t)9FUlE)/0ET JET10A80 

EP( J4n-CP( J)4RCLAX90ELE JCTS0A90 

C<J4t)>C< JKRELAX90ELC JET10700 

90S FORHAT (ISXf ITERATION CYCLE* ' r I 3t 2Xt ' EPS* ' > F tO .3* 2X> ' SI6HA* ^ f FIO. JET10710 

13> JCT1072V 

IF (APS(rELC).LT.l.E-3UND.ASS(0ELC) *LT.l .E-S> IFLAO*! JCT10730 

IF (IFLAG.EO.ll URITE (6»909) JET1O7A0 

909 FC.'^HAT ( 20X f ' SOLUTION OF GROUND PRESSURE DISTRIDUTION NAS »£EH FOUJCT10730 

1ND*> JET107A0 

IF (IFLAG.EO.n URITE (Af903) J* EP ( J U > f C ( J41 ) Jfn0770 

900 CONTINUE JCTt07EQ 

STOP JET10790 

907 $IO*C(J) JETIOUOO 

EPS«EP( J) JCT108A0 

CUU-SIOtCUl JET 10820 

WRITE (0rl3) cUU JCT10830 

13 FORHAT <20X»'UPUASH THICKNESS CONSTANTt <CU) X (SfOHA ) » * »FIC*3) JET10840 

IF <1P8AR.ME.0) CALL 6PL0T t ALPS » DELPS » SD> RO* EPS »CUI fSI6tPHIUO»PHIOJCTlO«SO 
1) JET108A0 

RETURN JET10870 

END JET 10880 

JCT10890 

JC710900 

JET10910 

JET10920 

SU8R0UTINE GPLOT ( ALPOt DELPSt SD»RO»EPS r CUI r SI6I t PHIUOrPHIO > JET10930 

MHEHSION PU(23) rRUALLI 102) »CRRa02) »XPL07U<202) r YPLOru( JET10940 

1202) tXPL0Tj(202> »YPLOTJi202>tXAXlS<30) «YAXIS(30) >XDATA{10)tYDATA<lJET10930 
20) JET10900 

COHHON /HEIGHT/ HDrIPDAR JCT10970 

JET10980 
JCT10990 
JETllOOQ 
JCrilOiO 
JCT11020 
JETU030 
JETU040 
JETllOSQ 
JETI1060 
JET11070 
JET11080 

CONFUTATION OF 7U0-JET OROUND IS08AR PATTEJETU090 


THIS ROUTINE COMPUTES THE GROUND ISOBAR PATTERN FOR THE TWO- JET 
IHPINGENENT FLOW FIELD 

FETA<ETA»A)-( 1 . -ETAttA) t94 

INPUT NUMBER OF ISOBAR VALUES 


READ(3*5t) NU 
31 format <I2) 

URITE (At33) NU 
33 FORHAT K K r //30X » ' « 9f 


IRN 4i$'///32XtI2r ' 
2 //) 


VALUES OF PRESSURE SPECIFIED FOR PATTERN' 


INPUT NU ISOBAR VALUES FOR GROUND PATTERN 

READOfl) (PU(t> f !*1 >NU) 

I FORMAT (F10.3) 

WRITE (4tU2) (PU(I) tKltNU) 


JETlllOO 

JCTIUIO 

JET11120 

JETU130 

JET11140 

JETlllSO 

JET11160 

JET11170 


112 FORMAT (lOXi 'PBAR*' rB<Fl0.3t2X)) 
URITE (6>121) 

121 FORMAT (//23X» 'GROUND PATTERN IN JET 
i> 

HOPTS-29 

C CALL PL0T(4. «0*t-3) 

XLMAX*ROfSD 
C X8C-XLMAX/8. 

XSC-I.O 

XHIN-0. 

DXH|N«-K 
DO 801 IP-IMOO 
IF (XMIN.LT.'RO) 00 TO 802 
XnlN*XMlN4DXHIN 
80 1 CONTINUE 


JCT11180 

JETU190 

CENTERED COORDINATE SYSTEM' //JE711200 

JCT11210 
JET11220 
JET11230 
JET11240 
JET 11230 
JCT11200 
JET11270 
JCTU280 
JET11290 
JETU300 
JCTIITIO 
JET11320 


Figure 69 - Continued 


154 



U o 


ORIG'NAL PACE tS 
OF PC JR QUALITY 


802 XHAX>SP 

0XAXlS«S0-XHli4 
NAnS»DXAXlS4l 
XAXIS( 1 >-XHtH/XSC 
YAX1S< 1 >»0. 

DO 809 1P*2>NAPTS 
XAXIS(IP>«XAX1S< IP-1>41 ./XSC 
YAXIS< IP)»YAX15< IP-1) 

809 CONTlHUe 

DYAX-l./XSC i 

XAXIS<NAPTS4| )>SD/XSC 
YAXIS(NAPTS4l)«0. 

DO 0OS 1P-2»MAPTS 
XAXIS<MAPTS4IP)*SD/X5C 
YAXIS(MAPTS41P)=YAXIS(NAPTS4IP-1>4PYAX 
605 CONTINUE 

MPL0T-24NAPTS 

CALL LINE (XAXtSrYAXI$«NPLOT»lf l»-3»l».l) 

CALL SYHDL4 (-3, r8. . .2. .Of S/D-2 . 00' » 0 . • 18 ) 
Pl»3* M159245 
XPTS*NUPTS-I . 

RELAX*. 5 
DO 500 I-l.NU 
WRITE <6f66) PU(I> 

66 fORHAT <///45Xf 'FDAR*' ,F10.5) 

IFLAC»I 

XPUF»0. 

IJ€T«1 
IU-1 
RU«SD 
ALPUQ*1 .S 

CALL UALLJ ( RU » ^HN» DUH2 » DUri3 • PUN4 1 DELPU) 
PNAXO>*DCLPU 

IF (PHAXO.LT.PU( I ) ) IFLAG«0 
RU«SD/COS (PHIUO) 

CALL UALLJ f RUtUMN f DUN2 i DUH3 > DUH4 f DELPU) 

PnAXF'0ELPW-UHN492f (UHNtCUS(FMIUO) )9«2 

IF (IFLAQ.eO.O) CO TO 502 

RUaLLU ) -SD 

RUALU2)^SD*.l 

DO 400 K>^l»100 

CALL WALL J( RUALL < K ) » OHH* DUN2 t MJri3 » DUH4 f DELPU ) 
XU**SORT<RUALL<K)«f 2-5Df >2) 

PHlU«ATAN2<XUrSD} 

PHAX»DELPU-VHN««24(VNN«COS(PHXU> ) 4t2 

ERR<K ^ =PU( I ) -PNAX 

IF (K.EQ.l) GO TO 400 

IF <ADS<€RR<IO).LT.l.E-5) CO TO SOI 

PRDE«<RUALL(K>-RUALL(K-| ) )/(ERR(K)-ERR(K-n ) 

RUALL(R41)«RWALL(K)-RELAX«DRDEFERR(X) 

400 CONTINUE 
S’OP 

501 RU*KUALL(K) 

TPUF = SORT(RW*t2-St«<»2) 

PUF»ATAN2(XPUFf SP> 

YUF*CUI*SIGItSP/COS<FUF) 

PUFP»PUF*180./PI 

502 OPMIU-PHIUO/XPTS 

IF (PHIUO.EQ.A. ^ GO TO 201 
PHIU*0. 

DO 200 IJ*1»NUPTS 
RW«SD/COS<PMlU> ' 

CALL UALLJ (RUtVHNtpUH2tPUH3fPUN4tPELFU) 
YU*^CU14SI0I*SD/C0S(FH1U) 

XU*SD4TAM(PHIU) 

PHI*ATAN2<XUf SD-YU) 

RU*<SD-tU)/COS<FHl) 

R0P«R04EPSf ( 1 . FCOSCPHU ) 

ETAN*RU/ROP 

IF (ETAH.GT.l.) FETAG^^O.O 
IF (ETAH.GT.l.) GO TO 1011 
FCTAG«FCTA(ETAH#ALPG) 

1011 PN1N*DELPS«FETA0 

PhAX*PELPU-UHM*i2f ( yHH«COS(PHlU) )••: 

CALL PNATCH (FNINfPHAX. ALPUG) 

PDAR*(PU* 1 > -P«IN)/(P?1AX-PnlH) 

Figure 69 - Continued 


JET11330 
JETU340 
JET 11330 
JET11360 
JET11370 
JET 11 180 
JET11390 
JET 11 400 
JCrtl410 
JET11420 
JET11430 
JET11440 
JET1I4S0 
JETn460 
JETU470 
JET 11 480 
JETU490 
JETU300 
JETUSIO 
JETUS20 
JETllSSO 
JET11540 
JETU;50 
JETU360 
JET11570 
jETitsao 
JETUS90 
JET11600 
JET1I610 
JET 11 620 
JET11630 
JETU640 
JET11630 
JET1I660 
JIT11670 
J£rU680 
JET 11690 
JET11700 
JET117tO 
JET 11 720 
JET11730 
JET11740 
JET11750 
JETU760 
JET11770 
JET U 780 
JETU790 

jeniaoo 

JETlieiO 
JET11820 
JET11830 
JETl 1840 
JETtieSO 
JET11860 
.JET11870 
JETt I860 
JET118;0 
JETU900 
JET11910 
JETl 1920 
JET11930 
JETl 1940 
JETI1930 
JET11960 
JCTU970 
T1I980 
jtMl990 
JET12000 
JET12010 
JET12020 
JET12030 
JET12040 
JET12030 
JET12060 
JET12070 
JET12080 
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OF 
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quality 


IF {PH1M.QT.PU<1)> GO TO 509 
PBARJ«PU<I)/OELPS 
IF <PeARJ.07.1.) GO TO 309 
CTAJ*a.-F»ARj«f .23>t«(l./ALFQ) 

XPLOTJ< lJ£T>>R0PteTAJ9C0S(PHI) 
YPL0TJ<IJET)«R0PtCTAJ9SIH(PHI) 

XPLOTJ< IJET)-XPL0TJ( 1JET)/XSC 
YPtOTJ( 1JET)»YPL0TJ( IJET)/XSC 
IJET-IJET41 

309 IF <PBAR.OT* 1 . .0R.P9AR.lt. 0. ) GO TO 306 
CTAU«< 1 .'PBAR«t.2S)9«a ./ALPUO) 


JET12090 

JET12iOO 

J£T"2U0 

JET12120 

JET12130 

JETI2140 

JCT12130 

JET12I60 

JET12170 

JETI2180 

JCT12190 


XP*.OTU< IU)*Sl*-YU*ETAU 
rPLOTU< IU>»XU 
XPtnTU< 1U)«XPL0TU( 1U>/X$C 
YPLOTOt IU)»YPL0TU( IU)/XSC 
lU^IUM 

300 PHIU^PN1U4DPHIU 

200 CONTINUE 

201 DPMI''(PI-PHI0)/XPTS 
PHI^PHIO 

DO 600 lK«tfNUPTS 
PDARJ«PU< n/DELPS 
IF (PDARJ.OT.I.) 00 TO 600 
CTAJ«< 1 .•PDARJ99.23)t«(t ./ALPO) 
R0P-R0fEP8t< 1 . ♦C09<PHIM 
XPLOTJ< lJET>>ROP«ETAJtCOS(PHl) 
YPL0TJ(IJET)>R0P9ETAJ$SIN<PH1> 
XPLOrJ(IJET)*XPLOTJ(IJ£T)/XSC 
YPLOTJ< IJET)»YFLOTJ< IJET)/XSC 
IJET-IJET41 
PM1*PMII0PMI 
600 CONTINUE 

IF (IFLA6.E0.0) GO TO 313 
IF (PUin.EQ.O.) GO TO 313 
IF<POF.LT.PHlUO) GO TO 510 
00 TO 311 

310 XPLOTUl tU)«$D 
YPLOTUl IU)«XPUF 

XPLOTUI IU)«XPL0TU< IU)/XSC 
YPLOTUl IU)»YPL0TU< IU)/XSC 
CO TO 600 

311 DPH1U«<PUF-PHIU0)/X« TS 
PHIU>PHIUO 

DO 700 IJ»1»NUPTS 
YU*CUlt$10I«SD/C0S(PHIU) 

XU«SD«TAN(PHIU> 

RU«SD/C0S<PHIU) 

CALL UALLJCRUfVHN«DUH2»DUN3f DUn4«DELPU) 
PHAX*0EI PU-VHN6t24(V«N*C0S(PMIU) )9*2 
PBAR^PUn ) /PMAX 

ETAU«< 1 PBARt«.23)t«< 1 ./ALPUG) 

XPLOTUl 1U)»SI»-YU$£TAU 
YPLOTU< IU)*^XO 


JETI2200 
JET12210 
JET12220 
JETI2230 
JET12240 
JET 1 2230 
JET12260 
JETI2270 
JETI22eO 
JET 12290 
JCT12300 
JET12310 
JET12320 
JET12330 
JET12340 
JET12330 
JET12360 
JET12370 
JET 12380 
JET12390 
JET12400 
JET12410 
JET12420 
JET12430 
JET12440 
JET124S0 
JET12460 
JET12470 
JET12480 
JET12490 
JET12500 
JET12310 
JET 12320 
JET12S30 
JET12540 
JET12330 
JET12S60 
JET12370 
JET 12380 
JET12590 
JET 12600 
JET12610 


XFLOTU< IU)*XPLOTU< IU)/XSC 
YPL0TU( IU)*YPL0TU( lUl/XSC 
IU«IUfl 

PHIUaPHIUlDFHIU 
700 COHTiNur 
313 IU*IU-’ 

800 lJET*lJt -I 

UR1T£ <6ilD lJ£T»j(U 

11 FORMAT (//ISXr'JET IMPINGEMENT RFGIONf 


JET12620 

JET12630 

JET12640 

JET126S0 

JET12660 

JET12670 

JET12680 

JET12690 

I JET » ' r I 4 * I OX p ' UPUASH OEFL JET 12700 


lECTlON REGIONt IU«'*I4//> JET12710 

IF (IJET.OT.O.OR.IU.GT.O) URITE <6pS32) JET12720 

332 FORMAT < //34X» 'XISOJ' f 6 Xp ' YISOJ' p8X p 'XISOU ' tSXi ' YI80U V) JCT12730 

INlN«tU JETI2740 

IMAX'IJET JET12730 

IF (lU.GT.lJET) 1NAX«IU JET12760 

IF (lU.GT.lJET) IMIN-IJET JET12770 

DO 803 IP>1pIMAX jeTt2780 

IF (IM.LE.IMIN) URIYE<6tl2) XPLOT J< IM ) i YPLOT J . .H ) p XPLOTU< IM > • YPLOT JET 1 2790 
1U<1M) JET12800 

12 FORMAT <30Xp4(Fl0.3t2X) ) JET12810 

IF (IM.CT. IMIN./.IO, IMIN.EO. lU) URITE <6rlS> XPLOT J( IM ) p YPLOT J< IH) JET12820 

15 format (30X.2(F10.3p2X) ) JET12B30 

IF (tU.OT. IMIN.AND. IM1N.F0. IJFT) WRITE (Ap 16> XPI OTU ( I H ) p TPLOTU ( IN JE T 1 2840 

1> ' JET12930 

16 FORMAT (S4Xf2(F10.3t2X) ) JET12860 
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( CONTINUE 

CALL LINE ncPLOrUtTPLOrUi lUf 1 * t I 1 I 1 I . U 
CALL LINE (XPLOTJ.YPLOTJ.IJET»lFl»\.l».i) 

) CONTINUE 
DPHIU«.5*PI/XPTS 
PKlU-0* 

DO 900 IJ-l#NUPTS 
YU«CUI«8IGI9S0/CQS<PHIU> 

XU*SD»TAN(PHIU> 

IF <XU.OT.XLHAX> GO TO 901 
XPL0TU(1J)-SD-YU 
YPLOTU(IJ)-XU 
XPLOTU(IJ)»XPLOTU<lJ)/XSC 
YPLOTUU J)-YPLOTU< IJ)/XSC 
PHIU-PHIUiDPHIU 
) CONTINUE 
[ lU-IJ-1 
URITE (6iB07> 

^ fOPHAT (///ISXr 'UPWASH DEFLECTION ZONE LINE» PDAR-0 OUTSIDE 
ICTION KEGIONV/) 

WRITE (AtdOai 

) FORHATi (3SXf 'XUP' »0Xr 'YUPV/) 

DC B13 IH>lfIU 

WRITE (A»01l) XPLQTUMH) fYPLOTU( IN) 

[ CONTINUE 

i FOftNAT (30X.2(F10.5»2X)) 

CALL LINE <XPLOTUfYPLOTU.lUfl*l.i.l. .1) 

CALL ADRAU 

CALL PLOT <-99. f-99. ► 3) 

READ <S»i01) PAUSE 
1 FORMAT <IX»F10.5) 

RETURN 

END 


SUDROUriNE WALL J<Rf UHNtDUHf DU»CU2rPELPU) 

REAL KDELO»NOtM<ELFD»NFDtNf KDEL 

THIS COMPUTES THE WALL JET PROPERTIES OIUEN A WALL RADIUS 

COHNON/WALL/ DELS r NO » KPELO » ALPG * ALPWO * ALPWFD » RO i RGHr UO • DUHO 

riAXlHUH UELOCITY DISTRIDUTION IN DEFLECTION REGION 

AU3« .09 
NFD*7. 

KDELFD»1 ./9. 

BU3HO«AU3«RO 

DELDRH><DU3rv.r-DUHQ)/(R0-RGH) 

ALPWFD=>t.5 

DELAR-(ALPUFD-ALPWO) /<RO-RGH) 

DELNR» (MFD-NO>/(RO-ROH) 

ADELDT* .0175 

ADeLDL*(APELBTt(RO-C. ) -DELS ) / < RO-RGM > 
DKDELR»(RDELFD-KDELO>/(RO-RGH) 

IF (R.LE.RO) F:TA«U.-(R/RO>f»ALPG)**A 
IF (R.OT.RO) FETA^O.O 
IF (R.LE.RGH) VMG*SQRT( I . -FETA) 

IF (R.LE.RGH) UMN»UMG«VG 
IF (R.LE.RGH) 00 TO AOl 
IF (R.LE.RO) 0WH»DWHO^D£LBRH*(R-RCH) 

IF (R.OT.RO) BUH«AW3tr 

IF (R.LE.RO) ALPU«ALPUOiDELAR«(R-RGH) 

IF (R.GT.RO) ALPU>ALPUFD 
IF (R.LE.RO) N=N04PCLNRF(P-RGH) 

IF (R.OT.RO) M^NFD 

IF (R.LE.RO) DELDL-DELStADeLDL*(R-RGH) 

IF (R.OT.RO) DELDL--ADELBT*(R-2. > 

ALAN- t (2. -S0RT(2. ) )/2. ) t|{ 1 ./ALPW) 

KDEL*DELBL/( DELBLt< DUH-DELBL ) /ALAM) 

CALL £INU(KDEL .Nf ALPH.RAT »C J2»CP) 

IF (R.LE . A . ) VISMOH-- 1 .0 

IF (R .GT . « . ) VISNON*! ./(R/4.)*l,24 

FI*VlSMOM*RATf ( I .-FETA>/(R*PWH) 

F2«CP*FETAtVG»i: 


JET12870 
JET12809 
JET12890 
JE712900 
JET12910 
JET 1 2920 
JET12930 
JET12940 
JET 1 2950 
JET12960 
JET12970 
JET12980 
JET12990 
JET13000 
JET13010 
JET13020 
JETI3030 
JET13040 
INTERAJET13050 
JET13060 
JET13070 
JET13080 
JET 13090 
JET13100 
JET13110 
JET13120 
JET13130 
JET13140 
JET13150 
JET13160 
JET13170 
JET1310O 
JCT13190 
JET13200 
JET13210 
JET13220 
JET1323A 
JET13240 
JET13250 
JETI3260 
JCT13270 
JETI3280 
JET13290 
JET13300 
JET13310 
JET13320 
JET13330 
JET13340 
JET133SO 
JET133A0 
JET13370 
JETJ33CC 
JET 1 3390 
JET13400 
JET13410 
JET1342Q 
JET13430 
JCT13440 
JET1345C 
JET134&0 
JET13470 
JET1340O 
JET13490 
JET13500 
JET13S10 
JET13520 
JET13530 
JET13540 
JET 13550 
JET13560 
JE 113570 
JET135B0 
JETI3590 
JCT13A00 
JET13610 
JET13620 
JET13A30 
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yHN*S0RT( (FI-F2)/(?.$CV2) ‘ JET 13640 

VHG»VHH/V6 JET13630 

DELPW«fETA#V0**2 + ‘;r‘S442 JCT 13660 

fWa^UM/RAT JET 13670 

DEXTA-KOELtOU JET13680 

401 CQNTlHUE JCT13690 

IF (R.LT.ROH) PELFW»VO**C JET13700 

IF (R.LT.RGM) PVlM«BiyHO JET13710 

RETURN JET13720 


END 


•U»R0UT1HE JCTFCT < TH • ETAC » PVH # AtPV » TL AN » CH » CT « RV » RT ) 
OELNlN-l.E-5 

CAUL JETfCNtETACrBVHiALRVt »RV1»CV2) 

C WRITE C6iSSS> CW2»AtPVI 

C CALL JETPC <ETAC»BVH»ALPVI»R0I»CV2) 

C WRITE <6»SSB) CV2»ALPVX 

555 FORHAT ( IX t ' C02« ' » F 10 . 6i 2X * ' ALPVI • ' »Fl 0 . 6 ) 

C WRITE (6»123) TflrETACtEVH 

123 FORHAT (IX 
ALPV*l.«ALPv: 

TLAM-l. M 

R2»<80RT<2. )«l » )/S0RT(2.> 

REL«. 1000000 
00 100 X«ir999 

RV-ETAC4U .-£TAC>tR2»t(l./ALPV) 

RVAV»(ETAC-1 . )t( (RV-ETAC)/( I . -ETAC) )4AL0Q< (RV-ETAC )/( 1 .-ETAC) )/ 
lALPW 

CALL SI NTPC ( TH » ETAC » ALPV > TLAM rCHfCT»CNAW> CHAT fCTAMiCTAT) 

F«2.T0VH$»VH6CT-RVtRV 

0-2 * «RVH»BWHtCH-RV$ltV 

FAV-2. TP0HiPWHtCTAM>2. FROfROAW 

FAT-2. *BMHABVHaCTAT 

OAV-2 . $OVHf BVHtCNAW>2 . tROtRUAW 

OAT-2. tfiVH«BVH9CHAT 

0ET*FAg«6AT-GAWlFAT 

PALPV-(0*FAT“FtGAT>/DET 

DLAn-(FtGAO-6tFAO>/OET 

C WRITE <6»600) BOH »CT »RO»CH r ALPWr ALPT r F t 0 
BOO FORHAT <lXfBFl4.6) 

C WRITE (6*700) I • DALPO * PALPT 

700 FORHAT < 1 X » 1 3 * 2X t 2F 15 . 6 ) 

ALPV-ALPVTRELIDALPV 

TLAM-TLAH^RELFPLAH 

IF (AI>S(F ) .LT.OELHIN.ANO.ABS(G) .LT .DELHIN) 00 TO 200 
100 CONTINUE 
STOP 

C 200 WRITE <6rS00) I.ALPU*ALFT 

200 cohtihl: 

500 FORHAT < 1 X » 1 3 » 2X » ' ALFV- ' *F 1 0 . 6 * 2X * ' ALPT* ' *F 10 . 6 ) 

RT«RU 

RETURN 

END 


SUBROUTINE JETKN ( ETAC rR JMi ALPC * R . CU2 ) 

C 

C THIS ROUTINE COHPUTES THE EXPONENT OF THE JET UELOCITY PROFILE 
C IN TME>0TEH’»1AL CORE REGION 
C 

ALP-4.0 
DCLHIN-I .E-5 
REL'^.SO 
no 100 I*. *500 

CALL SIHN(ETAC.i -P*R*CV2#CU2^ *RALP) 

FUNC^CV2~.5t(R/RjH>tl2 
FUNCD*CU2A- (RtRALP) MRJH4 42) 

alpn.alp-func/funco 

C WRITE (6»700> ALPN»FUNC 

700 FORHAT ( IX* ' ALf < ' «F10.5*2X* 'FUNC-' *f lO.S) 

ALP-ALPfRELK ALFN-ALP) 

IF < AIS(FUHC) .LT .DELHIM) 00 TO 101 


JET13730 
JET13740 
JET 13750 
JET13760 
JET13770 
JET 13780 
JET1379Q 
JET13H00 
JCri38tG 
JET13820 
JET13830 
JET13840 
JET 13850 
J€T13e60 
JET13870 
JET13880 
JET138tO 
JET13F00 
J2T13910 
JET13f20 
JET13930 
JET13Y40 
JETi3f3« 
JCT13Y60 
JET13Y70 
JET13Y80 
JCT13YY0 
JET14000 
JCT14010 
JET14010 
JET14030 
JCT14040 
JET14050 
JET14060 
JET 14070 
JET14O0O 
JET140Y0 
JET14100 
JET14110 
JET14120 
JETM130 
JET14140 
JET14150 
JE 114160 
JET14170 
JCT14180 
Jtf ;<190 
JET14200 
JET14710 
JCT14220 
JET14230 
JET14240 
JET14230 
JCT14260 
JET14270 
JET14280 
JET142Y0 
JET14300 
JET14310 
JET11320 
JET14330 
JET14340 
JET14350 
JET14360 
JET14370 
JCTM380 
JET14390 
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too CONTINUE 
STOP 

101 ALPC^ALPN 
RETURN 
END 


SURROUYINC SlHN(ETAC»ALPfR»CU2>CV2ArRALP) 

C 

C INTEGRAL OF ^ET UELOCITY PROFILE FUNCTI0N-CV2 

6 rbratio of rjH to ij 
C 

Al*,3-<» ^2. )F3*/(ALPtl . 1 / <3 . FALPF7 . > ♦ I . / < ^ . SALPF2 > 

A2ri ,-4./<HLPtl. )*6./(2.*ALPFl . )-4./(3.$ALPFl « ) 4 1 . / ( 4 . OALPF 1 . ) 
CU2».S«ETACt«2tA19(l ihET AC ) »$ 24ETAC> ( 1 . -ETAC ) SA2 
AIALP=»4./UALF*2. )442)-3./( (ALP4l)*421 K2./< (3.4ALPI2. )4t2> 

1-4. /( (4.4ALP42. )$I2) 

A2ALP'4./< <ALP+ntf2>-l2./( <?r*ALP41 * ) 4*2 ) + l 2 . /< i 3 . 4ALP41 .)$»?) 
1-4. /( (4.4ALP+1 )*•?> 

ETAP1»< 1 .-ETAC)*42 
ETAP2*ETAC4< 1 .-ETAC) 

CU2A«CTA>19AlALP4ETAli2tA2ALP 

R-ETAC4( l.-ETAC)*( <2.-S0RT<2. ) )/2. )44< 1 ./ALP» 

RALP«( (ETAC-R)/ALP)$ALOG< ( R-ETAC ) / ( 1 . -ET AC) ) 

RETURN 

END 


SUBROUTINE SINTFC ( TN • ETAC r ALPUf TL AH » CHf CT rCHALP » CHTt AH* CTALP r 
ICTTLAH) 

FUP(E'^AP)*< 1 .-ETAPI4ALPV)tt2 
FTF (ETAPT)-( I .-£TAPTMALPT)«2 

FUAVP(ETAF*FV)«2.»S0RT(FV)»(S0RT(FV>-1 . )*ALOO(ErAP) 


JET14400 
JCT14410 
jET: «?420 
JET14430 
JET14440 
JET14450 
JETt44A0 
JET14470 
JET14460 
JET14490 
JET14S00 
JET14S10 
JET14S20 
JET14330 
JET14S40 
JET14S90 
JET14360 
JET14370 
JETMSeo 
JET14390 
JET14400 
JET14414 
JEU4620 
JET 14630 
JET14640 
JET} 4630 
JCT14660 
JET14670 
JET146R0 
JET14696 
JET1470O 
JET14710 
JET 1 4720 
JET14730 
4ET14740 


C 


FTATP(ETAPT»FT)a2. tSORT(FT)V<SORT<FT >- 1 . XALOOt ETAPT > 
FTLAH(ETAPT>ETAP»TLAH)«2.*ALPTiM.-ETAPT«4ALPT>4(ErAPTtt(ALPT-l ♦ 
l))t(ETAf TLAM4TLAH)) 

F0HC(FU.FT)^FV/{FT4< 1 .-FT)/Tri) 

FUN(FT)«FT/(FT4( I .-FT)/TH) 

FUNCT<FV,FT)«FV/( <FT4(1 .-FI )/VH)»*2) 

N*24 

WRITE (6*e00) 7ilfETAC»ALPV»TLAH 
800 rORHAT (SX*4FU.3> 

!»TH*< 1 .-TH)/TN 
Nl-M-1 
ALPT-ALPU 
H2«N-2 

0ETA*<l.-ETAC)/Nl 

ETA-ETACfOETA 

SUH0H*0. 

SUNOHA«0* 

SUNOHT-0. 

SUHOTA«0. 

8UH0TT«0c 

SUHOT«0. 

DO too X«1 *N1 *2 
ETAP*(tTA-ETAC)/< 1 .-ETAC) 

ETAPT-CTAP/TLAH 

FV-FUP(ETAP) 

FT-FT*^<ETAPT) 

FUAV*FVAVP iCT AP *f V> 

FTA^*FTATP(ETAPT*FT) 

FTTLAN«FTLAN(ETAPT*ETAPt UAH) 


JET14730 
JCT14760 
JET14770 
JET147eO 
JET 1 4 >90 
JET14800 
JET1481C 
JETK82C 
JET14830 
JET14840 
JET148S4 
JCT14B60 
JETI4870 
JET14880 
JET14890 
JET14900 
JET14910 
JE 114920 
JET14930 
JET14940 
JEV 14950 
JET14960 
JET14970 
JET14980 
JET14990 
JEY13000 
JET15010 
JET15020 
JFT!5030 


SUH0H«SUH0H4FUtETAtFUNC<rU«F t ) 
SUnOT«SUNOT4FT«CTAiFUKC<FV*FT) 


JEM 5040 
JE1 ISOSO 


SUH0HA*SUH0KAi2. IFUAVtETAtFUMCCFV.FTlFFUtFTATfETAtPTMfFUNCTtFVtFT) JETIS060 
SUrtOHT«SUnOHT4FUtFTTI.AH»ETA«FUNCT(FV*FTi JEY15070 
5UN0TA-SUH0TA4FVAViETAfFUN^FT)iFTATfCTA*FUHC<FT .FU) ♦FTAT*PTHtFUNCTJET13O0O 
1<FV»FT)*ETA JET15090 
3UH0rT-SUM0TT + FTTLAM4ETA$(FUNC(FVfFT)4PTN4FUHCT(FV»FT) ) JET 15100 
FTA«ETA42.4DETA JET 151 10 
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CF PCw:\ 

100 CONTI HUE 

ETA>CTACf2. tOETA 
SUHCH«0. 

SUHET’Q, 

SUNENA«0. 

SUHEfif»0. 

SUNETA«0. 

SUHETT^O. 

DO 200 I^3rN2»2 

ETAF= iETA-ETAC)/< 1 . -ETAC) 

ETAF‘T*£TAP/TLAN 

rv*FVP<ETAP' 

ET«FTP(ETn 
FVW=FVAVP It TAP »FV) 


JCMS120 

JETI3130 

JETISMO 

jcnsiso 

JETlSiAO 
JETt3!70 
JCTlStao 
JET15190 
JET15200 
JETt3210 
JET15220 
JET1S230 
JET13240 
JET 15230 


f lAT^FTATPlETAPTtFT) 
FTTLAH-FTLAH<ETAPTrETAP»TL^H> 
SaHEH*8U«EN4FVfCTAfFUMCtFV»FT) 
SUntTsSUNETFFTAETAtFUHC(FVrFT) 


. CTt37A0 
JCT13270 
JET 13280 
JET 13280 


5UHEHA«SUNEHAf 2 . *F\;A0tETA9FUNC ( FU . FT : iFUtFTAT«ETAtDTH«FUNCT ( F Vf FT ) JET 13300 
SUN£NT«SUHEnTTFV«FTTLAN$ETA»FUNCT(FV»FT) JET 13310 
SUHETA»SUNETATFVA0fETA*fUN(FT>*FTATtCTAfFU\c<FT.FVliFTAT*DTMtF0MCTJETl3320 


KFV*FT)tETA 

SUNETT«SUHC r TTLANtE i At 'FUNC I FU.FT ) ♦DTHf FUNCT(FV»FT > ) 
ETA*ETAi2.t0c(A 
200 CONTINUE 

CN««3tETACtETAC*<DETA/3.)t(ETACTA.tSUttOHT2.tSUHElf) 
CT«.3tETACt£TACT<rETA/3. ) t < ETAC+4 . tr»JH0T*2 . tSUNET) 

CNAtF' I DET A/3 . > t M . tSUHOHA4^2 . tSUNENA ) 

CHTLA«*(DETA/3, )t( t . tSUH0nTF2.fSUHErtT)tDTtl 
CTALf >(tiETA/3. )t(4. tSUNOTA4 2.tSUHETA) 

CTT!.AH^(0ETA/2. ) t(4. tSUH0TTF2.CSUNETT) 

UKITE U*700) Tn>ETAC»AlPM»TLAH>CM.CT.CNALPrCHTLAN#CTALPtCrTtAH 
700 FORNAT <1X»IGF11.3> 

AC TUAN 
END 


JET 13330 
JCT13340 
JET 13330 
JET15340 
JET13370 
JET 13380 
JET15390 
JET1340Q 
JET13410 
JET13420 
JET15430 
JET13440 
JET134S0 
JCT13460 
JET13470 
JETIS480 
JETiS490 

SU8P0UTINE QHAEF ITHD^ET . ALPV»TtAN,'^Q> JET 13300 

FVP(ETAP'«U .-eTA^t*ALPV>tf2 JET1331C 

FTP(ETAP>»(1 .-(ETAP/TtAN)tiALPT>tt2 JET1S320 

FVETAPl£TAF>*-2.t< 1 .«CTAPttALPV)tALPVtETAPtt<ALPM-l . ♦ JET13530 

FTETAP'ETAPM-2. tl 1 .-(£TAP/TLAH)ttALPT)fALPTt( (£TAP/TLAH)tt(ALPT-l JET133^0 
1. 1 :tf 1 ./TLAN) JET 15330 

REL*.3 JET1S3A0 

AtPT«AlfV JET1S370 

C WRITE (4*7v0) THDf£TAc*AtPVf ALPT JET13S60 

700 FORNAT (in. 'T hD* * ,FlO,3t2X. 'ETAC*' »f l0.3»2Xr ' ALPV*' »Fl0.3.?Xt * ALPT JET 15390 
l»'*FlC-3> JET13AC0 

£TAP'>0* JLT13A19 

DETAP^l./24. JET15420 

DO 300 l>l»23 JET1S630 

FU»FVP(ETAP) JET13640 

FT=FTP<ETAP> JET1S6S0 

QOLD«G JET13040 

0*FvtFVtTND/( 1 . *FTt( TH6-1 . i) JET 13670 

IF (0*LT. .3.AMD.00LD.GT. ,5) GO TO 367 JET13680 

ETAraETAPfDnAP JET15690 

300 CONTINUE JET 15700 

STOP JET13710 

C 367 WRITE (6*337) QQLD^Q JET1S720 

337 FORHAI (2X, '00«' *F10.3*2X* *0*' #F10.3> JET15730 

367 DEDQ-DETaP/(Q-OOLD) JCT13740 

ETAP>ETAPTDEDQt( .5-0) JETI373 

C WRITE (6.373) ETAP JET13740 

173 FORHAT (2X. *ETAP«^ .FIO.3) JET15770 

DELnlM-l,t“3 JET15780 

DO 100 I«lf300 JET15790 

Fy*FUP(ETAP) JET1580C 

FT*FTP(ETAP) JET15810 

FUETA»FV€TAP(ETAP) JET13820 

FTETA-FTETAP(E^ AP) JETI3B30 

DEH-J , +FTt(T«D-l . ) JET13840 

FUNC» , 5-FUtFUtTHD/DEK JET138S0 

FUNCD*FVtFVtThDf ( THD-1 . ) tFTET A/ ( DE^tPEN) -2 . tFVtFVETAt TM8 /PEN JET 13860 

ETAH-ETAf -fUHC/FUHCD JET 13870 

ETAP-ETAP+RELt (ETAN-ETAP) JETiSSBO 

C WRITE (6.300) ETAP. nine JETI5P90 
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300 FQRHAT (1X*2F10.3> 

IF (A^SCFUMO -LT. I .E-5) 00 TO 200 
100 CONTIKUE 
STOP 

200 RQ«ETACF(1 .-ETAOtCTAP 
RETURN 
END . 


SUBROUTINE SIHTU < TN» TH> ALPU » TLAH> CM tCT » CHTNt CTTN) 

FVP(CTA)«< 1 .-£TA*0ALPV)«*2 
FTFCETA)«(1 .-(ETA/TLAH)ttALFT)tt2 
FUNC(ETA.FV»FT)=FVfETA/(FTtC I . -FT) /TNI*) 

FUNCP<ETA.FViFr) ^FU*ETA»( ( 1 . -FT ) /< TMBtTH ) > / < (FTF < 1 . -FT ) /THB ) ••2) 
M-24 

ALPT»ALPV 

C URITE (4t777) TM • TM » ALP\> » ALPT . CH . CT 

THB^THtTN 
Hl^N-l 


DETA«1./N1 

ETA=0ETA 

suNOn^o. 

SUN0T*0. 

sunOTN=o. 

SUH0CH»0. 

DO lOO 1*1. HI. 2 
FV*FVP(ETA) 

FT*FTP(ETA) 

SUHOM=SUHOH+FUtFUNC(ETA.FV,FT) 
3UH0T»SUH0nFTtFUNC(ETA.F .FT) 
SUH0Cn*SUM0CnfFV^rU?4CD(ETf .FUrFT) 
SUrtOTM*w*‘ " TfFTtFUHCP<ETAi FV.FT) 
ETA^ETA42. ?DETA 
100 CONTINUE 

ETA^2.tDETA 

SUMEH^O. 


JET1S900 

JET13910 

JET13920 

JET 13930 

JET1S9A0 

JETir950 

JET13960 

JET15970 

JET1S900 

JET15990 

JETIAOOO 

JCTlAOlO 

JET16020 

JET1A030 

JET1A040 

JET1A050 

JET14040 

JET1A070 

JET 16080 

JET 16090 

JET16100 

JET16110 

JET16120 

JET16130 

JCT16140 

JET161S0 

JET16160 

JET16170 

JET16180 

JET16190 

JET16200 

JET16210 

JET16220 

JET1623C 

JET16240 

JET 16230 

JET16260 

JET16270 


SUHE r*o. 

SUNECN=0. 

S*JNETH«0. 

DO 200 1-3»N2»2 
FV*rVP(£TA) 

ft-.ftp<eta> 

Sb.iEn*SUnCN4FUfFUNC<ETA.FV»FT) 
SUHET«SUHET4FT$FUHC<ETA.FV.FT) 
SUHECH*SUNECN4FV«FUNCD(ETA.FV FT) 
SU«ETN«SUHETNfFTtFUHCD<ETA.FV.F. ) 
ETA^ETA42.tDETA 
200 CONTINUE 

CN« DETA/3. )t( 4.«SUH0H42.tSUHEH) 
rT«CDETA/3. )t< 4.4SUH0T42.4SUMET) 
CNTH*< DETA/3. > t ( 4 . tSUNOCNi 2 . tSUNCCN ) 
CTTrt*C DETA/3. ) • ( 4 . «SU«0Tf1t2 . tSUnETH ) 
C URITE C6.777) TN. TN. ALPV» ALPT . CH* CT 

777 FORHAT ,'1X»6F10.5) 

RETURN 

END 


SUBROUTINE GPREST <RO»UG»ThGf ALPC.R.CPO) 

FUNC<AtP)» .5-4./<ALP42. )43. /<ALP4l . ) -4 . / 4 3 . tAtP42 . ) 41 . /< 4 . «ALP 
142. ) 

FUNC iM ALP)*4./< (ALP42. ) 442 ) -3 . / ( < ALP4 1 . ) 442 ) 4 1 2 . / < < 3 . 4ALP42 . ) 442 ) 
1-4. /( ( 4.4ALP42. )442) 

THIS ROUTINE SOLVES GROUND PRESSURE INTEGRAL FOR GROUND PRESSURE 
EXPONENT 

ALPG*1 .5 
REL».3 

PELNIN«1 .E-5 
DO 100 1*1 .49 
CPG-FL' C(ALPO) 

CPGP-4 jNCO(ALPO) 

F*THG-R04R04VG4O04CPG 

FP»' R04R04VG4UC4CPGP 

ALPC*A‘.PG-REL4F/FP 

URITE (6.500) ALPG.CPG*CPGP«. 


JET16280 

JETI6290 

JCT16300 

JET16310 

JET16320 

JET16330 

JETI6340 

JET163S0 

t£T163«0 

JET16370 

JET 16380 

JET16390 

JET16400 

JET16410 

JET16420 

JET16430 

JET16440 

JET16420 

JET16460 

JET16470 

JET16490 

JET16490 

JET16300 

JET16310 

JET 16520 

JET1633C 

JET16S40 

JET 16330 

JET16360 

JET16370 

JET16S80 

JET16S90 

JET16600 

JET16610 

JET 16620 

JET16630 

JET16640 

JET1663G 

JET16660 

JET16670 

JETU680 

-ET16690 
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300 FORMAT (IXr^FlS.S) 

If (AiSif > .L i . uetMiN) GO TO lOt 

100 COMTINUf 

srof 

101 RCrURN 
END 


SUiROUTiNC Slny(KDeL.»N»ALFU»Rf CV2»Cf > 


THIS ROUTINE COMPUTES THE INTEGRAL Of THE VELOCITY SQUARED PROFILE 
AND STATIC PRESSURE PROFILE FOR WALL JET FUNCTION^ 

CU2»VEL0CITY SQUARED INTEGRAL 
CP«STAIC PRESSURE INTEGRAL 
R«RAriO OF DWH TO »W 


JET16700 

JET1A710 

JET14720 

JEY14730 

JIT1A740 

JETIA7S0 

JET1A7A0 

JET14770 

JET1«7B0 

JET16790 

JET1A800 

JET16BI0 

JCT1A020 

JET1&830 

JET1A840 

JET1A8S0 

JETIABAO 

JCTU870 

JETU880 


REAL RDEL*N 
XN<N 

xil«l .- 4 ./< ALPUTl . >TA./( 2 .tA( PWTI • ) -4 . / < 3 . t A< P^f 1 . > f 1 . / ( 4 . t Al PH » I . JFT 1 AR90 
t> JCT1A900 

CU2«f XN/(2.4XH) >fKDEL4U . -KDEL >f XII 
CP«KDEL 4(1. -KDEL ) tX 1 1 
F*( (2.-S0RT( 2. ) T/2. )#t n ./ALP«) 

R-KDEL4(1 .-KDEDtF 
URITE (6*300) CU2tCP 

300 FORMAT I IX* ' CV2» ' *F10 .5. 2X i ' CP* ' *F 1 0 . 5 ) 

RETURN 
END 


JET16910 
JET I 6920 
JETI6930 
JET16940 
JCT16930 
JET 16960 
JET16970 
JET16980 
JETI6990 
JET 1 7000 
JET17010 
JET17020 
JET17030 
JET17040 
JET 17030 

CV2U*1 .-4./(ALPU41 . )46./(2.«ALPU41 . ) -4 . / ( 3 . 4ALPU41 . ) 4 1 . / ( 4 . tALPUEl JET 1 7060 
l.> JET 17070 

F^(SQRT(2. )-l . > /SORT (2. ) JET17080 

R*Ftt( I ./ALPUJ JET17090 

RETURN JET17100 

END JET17110 

JETX7120 

JET17130 

JET17140 

subroutine GMATT (TNfRGrDELPS.PHtI *TNI,VMlrALPO*ADELfN«ErAUG*DELS JET171S0 


SUBROUTINE S1MU6<ALPU*R*CV2U) 

THIS ROUTINE COMPUTES VELOCITY QUARED INTEGRAL FOR UPUASH 


1iALPU*R*BU0*BU0H.TLAH) 

THIS ROUTINE IS USED FUR INITIATING HALL JET REGION 

COMPUTES The INTIAL INVISCID EXPONENT FOR WALL JET PROFILE 

REAL KDEL.H 
»V0*DEL5/KDEL 
TLAM^l.l 
ALPU* .3 

CALL SIMUH (TN.THItTLAMfPMBI >KDEL *N« ALFa* ALPO*ETAUG*DELPS. VMT • 
lBy0»RQ> 

RC0N»<2.-S0RT(2. > )/2. 

R<KD£L4U . 'KDED tRCONtf < 1 . /ALF U) 

BUQH*RtI>WO 


TLAM. ALPU*BU0«BU0H 
. 'TLAN-' .FI0.3.2X. 'ALPW« ' tFlO.S* 2X» 


URITE ( ' - ' 
709 FORMAT ( ■ 
.F10.3) 
RETURN 
END 


SUBROUTINE SIMUM(TN»TMtTLAM*PMB*KDEL*Hf ALPU»ALPO*ETAUO*DELPS*VM i 
1BU.R6) 

REAL XD€L*N 

TM>«TM4fM 

FEL-.3 

DELHIM*! ,E-3 

^ETAG«( I . -CTAUG9fALP6U44 
FMOM-1 , -FETAG 

URITE (6*400) FETAGtFMOM*BU«RG 


JET17160 
JET17170 
JET17180 
JET17I90 
JET17200 
JET17210 
JET17220 
JET17230 
JCT17240 
JET17250 
JET17260 
JET17270 
JET17290 
JET17290 

JCT17300 

BWO-' *ri0.3*2X. 'BU0HJET17310 
JCT17320 
JCT17330 
JCT17340 
JE717350 
JET 17360 
JET17370 
JCT17380 
JET 17390 
JET17400 
JET17410 
JET17420 
JET17430 
JET17440 
JET17430 
JCT17440 
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400 FORMAT CIX» FETAO ‘ *4F10 .5) 

DO 100 1*1*200 

CALL SIMWTl <TMrTH*TLAM*PM|i.KDEL*H*ALRll*Cy2*CR*CT*CWrAiCV2L 
1 iCFA*CTA*CTL*R*RALP) 

WRITE <4*740) CV2.CP*CT*CV2A*CV2t*CPA#CTA*CTL.R 
740 FORMAT UX*9F10.3) 

F*FM0M-RCiiW4<2.«PM>tVM*VM$CV2/TM4DEtP8«FETA0fCP) 

0*FM0MtTH«( TN-1 . )-2. tPH»4VNt(THi-l . )iROt»W»CT 
FALf •'R0t0W9 ( 2 . tPMI«MMt0H«C02A/TMf DELPStFETAOiCPA ) 
FLAM-<RGt|>yt(2.9PMfitUMt0MtC02L/Tlt) 

0ACP**2. tTHOtVMttTMi'l . )SFG9RWtCTA 
6LAM--2 . ff MRtVMt < TMP- t . >«R6tRW9CTL 
WRITE (4*750) FAtP*FLAM*CALF*OLAM 

750 FORMAT MX. FA»‘ *F10.5*1X*'FL«'*F10.3*IX* *GA»' *Fl0.5*lX* *Gt*'*FlO 
13) 

0CT-FALPl6LAM-FLAM$GALF 

DALP«(u»FLAM-F90LAIi)/0CT 

»LAM»(FtGAtF>GtFALF)/!»ET 

ALPWALPU4REL9DALF 

TLAH-TLAHtRELtDLAM 

WRITE <4*500) F*0*AtPW*TLAM 

IF (AiS(F) .LT .DELM1N.AMD.A»S(G) .LT.OELMIN) GO TO 101 
300 FORMAT ( 1X*4F15.S) 

100 CONTINUE 
STOP 

tot CONTINUE 
RETURN 
END 


SURROUTINE SIMWTI (TN*TM*TLAM*PM>*K0Et*N*ALPV*C02tCP*CT»CV2A»CU2L 
I *CPA.CTA*CTt*RfkALF) 

REAL KDEL*N 


JET17470 
J€T174tO 
JET17490 
JET 17300 
JET17310 
JCT17320 
JET17534 
JET17340 
JET17S30 
JCT17340 
JET17370 
JET17300 
JET I 7590 
JCT17400 
JET17410 
JET17420 
JET’743« 
JET17440 
JET17430 
JET17440 
JET17470 
JET1740O 
JET17690 
JET17700 
JCT17710 
JET17720 
JET17730 
JET17740 
JET I 7730 
JET17740 
JET17770 
JET17700 
JET 17790 
JET17800 


FVF <ETAF )«( 1 . -ETAf 9«ACFU>tl2 
FTP(ETAF)»U .•{£TAP/TLAH)ttALPT)tt2 

FVA(ETAP) •-2.9(£TAF««ALPV)$( 1 .-ETAPt9AlPV>»AL0G(£TAF) 
FTA(ETAFT)*-2.t(ETAPT*tALF^)t( 1 . -ETAP199ALPT )«ALOG<ETAPT ) 
FPA(ETAF)*-4. •(ETAP*»ALFV>*( C ! -ETAF«9ALPV>$t3)4AtOC(ETAF) 
FTL<ETAPT)-2.9II . -ETAPT^^ALPT ) »Ai. ' rt (EIAPT#4ALPT)/TLAM 
FUNC(FV.FT*FP»«F0$(FP4< I . FP)/PH^ ' T4( 1 . -FT ) /TMi) 
FUNC0(FV.FT.FP)«FUt(FP4( 1 . FP) t (FT4U . -FT >/TMt>t«2 ) 

FUH(F0*FT>*FT9(FP+U . -F5 )/PMB>/ . ♦( I .-FT)/TMP) 

FONT<rV*FT ) *FV9FT/^FT>< 1 . *FT>/TffR) 

NPTS*24 

alpt=»alpv 

C WRITE <6*777) TN* TM * ALPU* ALPT *CM*CT 

rM9«THtTN 
DTM«<1 .-THR)/T«i 
0FM*(PHB-1 . ) /PMR 
Nl-NPTS-1 
N2«NPTS'2 
XNt«Nl 

0ETA-(1.-KPEL)/XN1 

ETA«XDEtfP£TA 

SUMOM-0. 

SUM0T«0, 

SCV2AO-0. 

*’CV2L0*0. 

SCTAO-0. 

SCTLO'O. 

00 100 l-WNl.2 
ETAP»(ETa-KD£L) /( l.-RDEL) 

CTAPT-ETAP/TLAM 
IF (ETAPT.OT . 1 . ) ETAPT*l. 

FU»FWP(ETAP) 

FT*FTP<ETAP) 

FP«FV*FV 
FUALP»FVA(ETAf ) 

FTALF‘*FTA<ETAFT> 

FPALP^FPACETAPJ 

FTLA=4-FTL<ETAFT) 

SUMOM»SUMO«fFUrFUNCCFO»FT.FP> 

SUMOT *SUMOT + FTtrUNC < FV*FT*FP> 


JET17014 
JET17820 
JET 17830 
JET17840 
JET 17830 
JET17840 
JET1*’870 
JET17880 
JET17890 
JET17900 
JET17910 
JET17920 
JET17930 
JET17940 
JET17930 
JET17940 
JET17970 
JET 17980 
JET17990 
JET 18000 
JET18010 
JCT18020 
JET18034 
JET10O4O 
JET18030 
JET18040 
JET18070 
JET 18080 
JET18090 
XT18100 
JET18110 
JET18120 
JET18130 
JET18140 
JET18130 
JET18140 
JET18170 
JET18180 
JET18190 
JET18200 
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SCV2A0*SCV2A0*FUNC(FViFT.FP)»< 2. »FUALP>FVfFPALP*IiPH> JET 18210 

UFUMCO(FV»FT.FF>tFVfFTAtP*DTH JET18220 

SCV2L0«SCV2L0^FUNCP(FV.FT»FP)*Fy*FTLAH«0TH JET 18230 

SCTAO»SCTAC*FVALF*FUHtFTf FF )frUMC(FVrFT.FP)«FTALP JET 16240 

lFFUNT(FU.FT)iFFALP«0FHFFUMCD<FO»FT.FP)fFTiFTALPF0TH JET I 8230 

5CTLOsSCTLOFFUKC(Fy.FT.FP>»FTLAM*FUNCD(FO.FT.FP) tFTtFTLAM«0TH JET182A0 

ETA-ETaF2.*0ETA JET18270 

100 C3NTIMUE JCT19280 

ETA’=hrELF2.*PETA JET 192^0 

SUMEn* 0 . JET18300 

SUHET«0. JET18310 

5CO2AE*0. JET18320 

SCV2LE^0. JET18330 

SCTAE^O. JET18340 

SCTLE«0. JET183S0 

DO 200 I«3tN2i2 JCT183A0 

ETAP«(ETA-KDEL)/< 1 .-KDEL> JET 18370 

ETAPT-ETAP/TLAM JET 18390 

IF (ETAPT*GT.l. > ETAPT^l. JET183tO 

FV*FVP<ETAP> JET 18400 

FT«FTP(ETAP) JET18410 

FP-FOAFV JFT18420 

FVaLP«FVA<ETAP> JET 18430 

FTAtP«FTA<ETAPT> JET18440 

FPALP*FPA<ETAP> JET I 8430 

FTLAH-FTL(ETAPT) J€Tia4A0 

SUHEH-SUHEMFFViFUNC ( FVt FT »FP ) JET 18470 

8UMET*SUM£T4FTtFUMC(FO»FTtFP) JET18480 

SCV2AE»SCU2AE4FUNC<FV»FT,FP)»<2.«FOALP4FVtFPALP#DPH) JET 18490 

I4rUMCD(fV»FTtFP>*FVtFTALP»PTM JET 18300 

SCV2tE*SCW2LE4FUMCD<FV»FT*FF )4FVtFTLA«4DT« JET 18310 

SCTAE-SCTAE4FUALPfFUH(FT»FP)4FUNC<FU#FT»FP>«FTAtP JET 18520 

14fUWT<FVfFT>#FFALPtDPK4FUHCD<FV,FT»FP)tFT«FTALPtDTH JET 18330 

SCTLE=SCTLE4FUKC < FO t FT rFP ) •FTLA«4FUHCD< FV. FT tFP)tFT«FTLAHfDTM JET 18340 

ETA«ETA42.iPFTA JET 18330 

200 CONTINUE JET 18360 

CV2»M4KDEt/( 2. ♦N)4(PETA/3. )f ( 1 . 44 . TSUM0H42 . tSUKEH ) JET 18570 

CT**H»KDEL/a .4N) 4 (DETa/3 . )«< I . F4 . «SUM0T42 . 4SUNET ) JET 18380 

CV2A*a»£TA/3. )#< 4 . •SCV2A0f2 .fSCy2AE) JET 18590 

CV2L*a»ETA/3. > f < 4 . 4SCV2L04 2 . •SCV2LE ) JET18600 

CTA*<l»ETA/3. )*<4.4SCTA042 .iSCTAE) JET 18610 

CTL <DETA/3. )t<4.«SCTL042.tSCTLE> JET 18620 

CPF*1 . - 4 . / ( ALPUf 1 . ) f6. /(2.4ACFV41 . ) - 4 . / ( 3 . t ALFV+ 1 . >41 . 4 . *ALPU41 . JET 18630 


CPF D«4 ./( ( ALPV41 (J.*ALf V4l .>tf2)4l2./( < 3.4ALPV41 . ) * 1 2 ) JF T 1 0650 
l-4./((4.4ALPV4X. )9t2> JET 18660 

CP-»M»EL4 ( I . ->PEl )»CPF JET 18670 

CFA>( 1 .-KPEDfCFFP JET18660 

RCON»(2.-SOKT<2. >>/2» JET18690 

R-KPEL4( 1 . -KPEL )*RCON«#< ! . / ALF V' JET18700 

RP*(R-MtEL)/M .-KPtL) JET 18710 

RALP«tKPEL-l . )4RP«AL0G(RF )/ALF U JET 18720 

yRITE (6»777> TN.TnrALPV»ALPT,Cv2.CPrCT JET18730 

777 FORMAT <1X»7F10.3) JET18740 

RETURN JET18730 

END JCT18760 

JET18770 

JET18780 

JET18790 

SU8RQUTINE SIMMTV ( 1 N * TM . Tl AH » P nD . K DEL t H t ALPU • CV2 » CP » CT t CV27M . CTTH JET IBPOO 
IrR) JET18810 

REAL KDEL»N JETX8820 

FUP<ETAP)'t(l . -ETAF ttALPV)«*2 JET 18830 

FTP<ETAP>»< 1 . >(ETAP/TLAH)4«ALPT)*t2 JET 18840 

FUNCCFU.FTfFP)«FVt(FP4( I . -FP ) /PH» ) / < FT4 < I »-FT)/TH8) JET 18630 

FUNCDiFV.FT.FP>»FV»( ( 1 . -FT )/( THDtTH >)•( FP4 ( 1 . -FP ) /PHD )/ < < FT4 ( 1 . -FT JET 1 8860 


1 }/TH|)tf2) 

NPT$«24 

alpt»alpv 

MRITE :6.777» TH.T 

THl*TH«TN 

DTH<»< 1 ,-TnD)/THD 

DPH»(PHi* 1 . > /FHP 

N1«NPTS-1 

N2-NPTS-2 

XN1*N1 

DETA*( I . -RPEL> 'XHl 

ETA*KDEL4DETa 

SUnOH-0. 


TH»TH.^LPV»ALP^ *CHiCT 


JET18870 

JET18680 

JETiee90 

JET18900 

JET18914 

JET18920 

JET18930 

JET1B940 

JET189S0 

JET18960 

JET18970 

JET1098O 

JET18990 
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SUHOT-0. 

SCV2T0*0. 

SCTNO^O. 

DO 100 l>lrNI»2 
ETAP-(ETA-K0eL)/« 1 .-KDCL> 

ETAPT«£TAP/TLAM 
IF (ETAPT.OT. 1. ) ETAPT«1. 

FV*FVP<£TAP) 
f T*FTPtETAP> 

FP*FOtFV 

SUnOrt«SUHOHFFV9FUHC(FOtFTfFF) 

SUHOT«SUMOTfFTtFUNC(FVfFT»FP) 

SCV2TO*SCV2TO+FUMCD(FV.FT.FP)9FV 
SCTHO-SCTHOfFUNCD<FV,FT.FP)*FT 
CTA«£TAf2.«D£TA 
100 CONTINUE 

ETA»KDCL42.tnETA 

SUHEH=0. 

SUnCT«0. 

SCV2TE-0. 

SCTHC^O. 

DO 200 I*3»N2»2 
ETAP»(ETA-KDEL )/( 1 .-KPEL ) 

ETAPT*ETAF/Tl AH 
IF (ETAPT*6T.l . ) ETAPT=1. 

FV»FVP<ETAf ) 

FT»FTP<ETAP) 

FP-FCtFV 

SUHEH»SUHEHTFV9FUNC<FV»FT»FP) 

SUHET*SUnETfFT$FUNC(FV.FTtFP) 

SCV2TC»SC02TETFUHCD<FVfFT.FP)tFU 
SCTHE»SCTHE*FUHCP(F0»FT.FP>fF7 
ETA^ETAT2. tDETA 

200 continue 

CV2^H»KDEL/(2.iN)T(0ETA/3. )9U . 44 . tSUHOHf 2t OSUHCH) 

CT'<HtKDEL/( l.iN)4iDETA/3. )«( K44*tSUH0Tf2.tSUHET) 

CV2TH-1DETA/3. )t(4.tSCV2T042.tSCV2TE) 

CTTH* (DCTA/3. > • ( A . tSCTHO^D . tSCTHE ) 

CPF •'I .-A./<ALPV41 . /<2.«AUfV4l . > - 4 . / < 3 . $AtPV4 I . )41 ./<4.tAtPV4l 

1 ) 

CP«KPEL4(1 .-KDEU9CPF 
RCOM«(2»-SOPT<2.)>/2. 

R*KDEL4<1 .-KDEDtRCOHftU ./ALFO) 

WRITE (A»777) TM. TH . ALPU . AtPT »CV2 »CP *CT 
777 FORHAT (IS.7F10.!;) 

RETURN 

END 


SUDROUTINE P6UNAT(CUFD»DELPStDELPW»R0f 3Dr ALPOr ALPUG) 

DINENSION CRRUOl ) »SIGU01 ) 

THIS ROUTINE CONFUTES EXPONENT OF UPWASH GROUND PRESSURE PROFILE 

SIG< 1 )». 10 
S1G(2)»1 .0 
ALPUFD«! .5 
DO 100 I>1«100 
FR»(DELPS/DELFW)t*( .23) 

ETAG>SDt\ I .-SIC(I> )/R0 
ETAU-SIC( n/CUFP 
Fl*PB«< 1 .'ETAG99ALPG) 

F2-1 .“ETAlittALf UFO 

ERR(ll*Fl-F2 

IF (I.EQ.l) GO TO 100 

S»IS1G« J )-SI5< I-l ) )/(ERRn )-ERR<I-l) ) 

IF <ADSIERRM ) ) .LT. I »E-5) CO TO 101 
S;G(141)*-SIG( I>-S*ERRUT 

too CONTINUE 
STOP 

101 SIGUO«StG<I> 

ETAH-SIGUO/CUFD 
IF (ETAH.GT . 1 .0) ETAH^l.O 
PH IN* DEL put ( 1 . >£TAKt«ALPUFD) tt4 
CALL PHATCH t PHI N . DELPU • ALF UG ) 

RETURN 

END 


JET 17000 
JET17010 
JET17020 
JET17030 
JET17040 
JET17050 
JET17000 
JET 17070 
JET170B0 
JETI7090 
JET17100 
JET17U0 
JET17120 
JET 19130 
JET19140 
JET19130 
JET19U0 
JET19170 
JET19190 
JET19190 
JET 19200 
JET19210 
JET19220 
JET 19230 
JET19240 
JET 19250 
JET192A0 
JET19270 
JET192B0 
JET19290 
JET19300 
JET19310 
JET19320 
JET17330 
JET 19340 
JET19350 
JET19340 
JET19370 
. JET19380 
JET19390 
JET19400 
JET19410 
JET19420 
JET19430 
JET19440 
JET 19450 
JET19460 
JET19470 
JET194B0 
JET19490 
JET 19500 
JETt95lO 
JET19320 
JET19530 
JET19540 
JET 19550 
JET 19500 
JET19570 
JET19580 
JET 19590 
JET19000 
JET19010 
JET19020 
JET19630 
JET19040 
JET19650 
JET19660 
JET19470 
JET19080 
JCT19690 
JET19700 
JET19710 
JET19720 
JET19730 
JET19740 
JET19750 
JET19740 
JET19770 


Figure 69 - Continued 


165 



r> n n 


SUI^kQU :N€ PHATCH <PHlNtPHAX.ALPUG) 

DinCNSlON CRR< lot > »AL<101 ) 
c 

C THIS ROUTINE COHPUTCS THE EXPONENT OF THE PRESSURE PROFILE FOR 


JET19780 

JET19790 

JET19800 

JET19810 

JET19B20 

JET19830 

JET19840 


C THE UPHASH DEFLECTION ZONE 
C 

ALPUFD-l.' 

P9-PHIN/PHAX 

ETAH*< 1 .-P8t«.2S)9«U ./ALPUFO) 

CALL SIHUF(ALPUFD»ETaN»CSPUFD) 

AL( 1 )«1.S 

AL(2)*3.0 

DO too I«t»10C 

CALL 91NUF<AL<I>>1.0«CSPU> 

Xl»ETAHt(PDf < 1 .-PD>tCSPU> 

X2-CSPUFD 
ERRU >-Xt-X2 

IF <ADS(ERR( I) ) .LT. l.E-S) GO TO 101 
IF (I.EO.l) GO TO 100 
S«(AL(U-AL(l-l) )/<CRR<n-ERR<I-n > 
AL(Hl>»AL<I>-S>ERRa> 

100 CONTINUE 
STOP 

101 ALPUO»AL(I) 

RETURN 

END 


SUDROUTINE SIHUF( ALPU*ETA»CSPU> 

Al-ETA 

A2«-M./(ALPU41. > )9ETAtt(ALPUil . ) 
A3*<A./12 .tALPUVl >>tETA»t<2.AALPUfl .) 
A4«-M. / (3»tALPU«l * ) ItETAtf (3. tALPUf 1 . ) 
AS«< 1 ./(4.$ALPUit . n«ETAtt(4.«ALPUf 1 . ) 
CSPU-AlfA24A3iA44AS 
RETURN'* 

END 


SU8R0UTINE SIHIETAC • ALP *R » CV2 > 

integral of jet VELOCITY PROFILE FUHCTI0N«CV2 
R«RATIO OF BJH TO IJ 

Al*.5-4./(ALP4 2. )f3./(ALPtl . ) - 4 . / ( 3 . •ALP + 2 . )♦! ,/( 4.tALPf?' 

A2’l .-4./< ALP41 . )fO./(2.4ALr+l . ) -4 . / < 3 . t ALP4 1 . ) ♦ 1 ./C4.4ALP+1 . > 
CV2*.5tETAC*»24Al»U .-ETAC)**24ETAC#( 1 .-ETA0 4A2 
R«ETAC4< 1 t-ETAC>f ( (2.-SQRT(2. ) >/2. )f$( 1 ./ALP) 

RETURN 

END 


JET 19830 
JET198A0 
JET 19870 
JET19880 
JET19890 
JET19900 
JETl99tO 
JET19920 
JET19930 
JET19940 
JET19950 
JET 1 9960 
JET 19970 
JET19980 
JET19990 
JET20000 
JET20010 
JET20020 
JET20D30 
JET20040 
JET20030 
JET20060 
JCT20070 
JET20080 
JET20090 
JET20100 
JET20110 
JET20120 
JET20130 
JET20140 
JET20130 
JET20160 
JET20170 
JET20180 
JCT20190 
JET20200 
JET20210 
JET20220 
JET20230 
JET20240 
JET20250 
JET20260 
JET20270 
JET20280 
JET20290 
JET20300 
JET20310 
JET20320 
JET20330 
JET2034O 
JET20330 


C 


C 


SUBROUTINE SINU< ALPUiRtCV2U) JCT20360 

CV2U3 1 , -4 , /(ALFU+1 . . /( 2. f ALPUf 1 . ) -4 , / ( 3 . »ALPU4 1 . ) ♦ 1 . / ( 4 . tALPUM JET20370 


1.) JCT20380 

F«(S0RT(2. )-l . )/SORT< ) JET20390 

R«Ft4il . /ALPU) JET20400 

RETURN JET20410 

END JCT20420 

JET20430 

JET20440 


A450 


SUBROUTINE Sinun(TN*TH»TLAN.PNB«ALPU»FETA6fDELPy$iVNf BUtROtRtFn »: 

ITNEFF tTAEFF) 0 

NRITE <Ailll) FhBiVn ‘ 

111 FORhAT (IX. PUB*' »F10.3»2X. 1 

TMB*TKtTN * 


PEL*. 3 'O^IJ 
IFLAG«“C r0__0 
DELH1N=»1 E-5 *0330 
TLAH* 1.0 JEr0540 
TLANI»TLAK JET20530 
URITE <6»400) FETAGtFHOMrBUtRQ JET20360 
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400 FORMAT < 1X> *FFTAG' »4F10.S) 
to 100 l<-ltS0 

CALL SIHUTl <TN.TH#TLAh.PHP.ALPU#CV2rCP.CT#CV2LtCTL#R) 

IF (I.eo.l) BU«FM0H/(RG«C02t(2. 

UVN<t240eLFUStF€TAG>> 

IF (I.eo.l) BUI^BU 

C WRJTE (A*740) CV2 » CF . C f t C V2L . CTL 
740 FORHAT <tX»6F10.S) 

C* IF (I.eo.l) URITE (&fl03) BU 

103 FORMAT UXf'BU EST . ^ ' f F 1 0 . 5 ) 

F*FMOM-RGr - 2 .»PhB«VHtVMtCV2/TM4DELPUSiFETAG$CF ) 

QrFMOMtTHK .;V*1 . > -2 . >P«P»VM* ( TMB-1 . )$RG*PUtCT 
FBU=-RGt(2.tfK|i»VM*VM<TV2/TM + PELPUS*FETAGiCP) 

FLAM--RG*I»U«(2. tPH|i*y«*VMtCv2L/TM) 

GBU^-Z . TMB-1 . )tRG»CT 

GLAH«'2.4PMIitvn«(TMB-l . )$RG«BU«CTL 
C WRITE (6»750) FBU»FLAH*GpU»OLAH 

750 FORMAT (IX. FA» '»F10.5.1X,‘FL=^'»M0.5flX» 'QA* ' tF10.5»lX» 'GL* 
13> 

OET-FfiU«GLAH-FLAN*GBU 
IF (ABS(I'ET).LT.l.e-lO) ILFAG^l 
IF (ABStOET) .LT . 1 .e-10) GO TO 101 
0BU*<6*FLAh-FtGLAMi /pfT 
DLAM*(FtC5«U“G»FPU)/t'tT 
BU«BUfREl*D»U 
lLA«»^TLAM4RELfDLAM 
IF (TLAM.LT.O.) IFLAG»l 
IF (TLAM.LT.O. ) GO TO 101 
C URITE (AtSOO) F.OtfiU.TLAN 

IF (ABS(F> .LT.DELMIN.AHD.ABS(G) .LT.DELMIM) GO TO lOl 


JET20370 
JET20S80 
JET20590 
JET2OA0O 
JET20610 
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JET20A30 
JET20640 
J£T20650 
JET20AA0 
JCT20670 
JET20&80 
JET20690 
JET20700 
JET20710 
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JET20730 
fFlO. JET20740 
JET20730 
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JET20770 
JET2079C 
JET20790 
JCT20800 
JET20810 
JET20820 
JET20830 
JET20840 
JET20850 
JET20860 


500 FORMAT (IX»4F15.5) 

ICO CONTINUE 
STOP 

101 IF (IFLAG.EQ.l) BU«BUI 

IF (IFLAG.eO.l) TLAM^TLAMF 

Rf TURN 

END 


SUBROUTINE SIMUTI ( TN» TM » TLAM » F MB » ALPU »CU2 i CP . CT » CU2L »CTL • R ) 
fVP(ETAP)-^( 1 . -ETAP*9ALPU)f92 
FTP(ETAP) ■( I . -(ETAF /TLAM)t*ALPT )tt2 

FTL(ETAPT)*2.t( I .-ETAPTttALPT)4ALPT»(ETAPTtf ALPT)/TLAM 
FUKC^ ^VfrTiFP>=FV»t(FPF( 1 . -FP ) /PMB ) / ( FT t ( I .-FT)/TMB) 
FUMCD(FV.FT»FP) ^FV^(Fr+( 1 ,-FP) /pMB)/( ( FTF ( 1 . -FT > /TMB) «*2 ) 
NPTS»24 
ALPT*ALPV 

C URITE (4#777> TW f TM . ALF V » ALPT »CM . CT 

TMBrTMtTN 
OTM»( t .-TMB)/TMB 
DPM»(PMB-1 ♦ )/PMB 
Hl*MPT5-l 
N2*NPTS-2 
XN1«N1 
OETA-1 ./XHl 
ETA*DETA 
SUM0M«0. 

SUMOT-0. 

SCV2L0»0. 

SCTLO»0. 

BO 100 I«ltNl»2 
ETAP*ETA 
ETAPT«ETAP/TLAM 
IF (E7APT.GT.1 . > ETaPT = 1. 

FV«FVP(ETAP) 

FT-FTP(ETAP) 

FP*FVtFV 

FTLAM*FTL(ETAPT) 

SUMOM:»SUNOMFFVtFUNC(FV.rT.FP) 

SUMOT-SUMOTfFTtFUNC<FV.FT,FP> 

SCV2L0*SCV2L04FUKLP(FV.FT.FF >tFV»FTLAMBDTM 

SCTL0«SCTL0FFUHC(FU,FT,FF >»FTLAh+FUMCO(FV»FT»FP)tFTtFTLAMtPTM 
CTA»ETA+2.tOETA 
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Figure 69 - Continued 



c 


c 


OWG\NAL 
OF POOR QUAUrr 


I 


100 CONTINUE 

ETA-2. tOCTA 
SUNEH-0. 

SUNET-0. 

SCV2LE-0. 

SCTLE-0. >. 

no 200 l«3»N2.2 

CTAP-CTA 

ETAPT-ETAP/TLAH 

IF (ETAPT.OT. 1. ) ETAPT-1. 

FU-FVP(ETAP) 

FT-fTP(ETAP) 

FP-FVffV 

ftlah-ftlcetapt) 

SUHEH-3UHEHfFV*FliNC(FV»FT*FP) 

SU«ET*SUHET4FTtFUNC(FV»FTtFP) 

SCV2LE*SCU2LEfFUNCD(Fy»FT.FF )*FMIFTLAM«PTH 
SCTLE«SCTLEfFUHC(FV.FTiFP)«FTLANfFUNCl*<FV»FTfFP>*FT$FTLAMtDTH 
ETA-ETAf2.». 

200 CONTINUE 

CV2-<PETA/3. )t( 1 .♦4.tSunON*2.«SUHFN) 

CT-(0ETA/3. >«a.i4,«SUnOTf2.$SUHETJ 
CV2L*(DETA/3. )t< A.iSCV2L0^2.*SCV2LE) 

CTt-(DETA/3. >«<4.«SCTL0F2.tSCTLE) 

.-4 ./<ALFVft . )+0. 2.4 AlPU*1 . ) - 4 , / < 3 . 4ALPU4 1 . ) ♦ I . / < 4 . 4ALPV4 1 

. i 

RC0N-(2.-S0PT<2.))/2. 

R»RC0Ht4(l./ALPU> 

WRITE <A»777) TNf TH.ALPV,AtPT.CV2.CPrCT 
777 FOPHAT <lXf7F10.5) 
return 

END 


SUOROUTINE SIHUTU ( TN# “H» TL AM • PN» • ALPU , CV2 » CP » CT i CV2TH t CTTH # R > 
FVP(ETAP)*( t .-ETAPttAtPVUt2 
FTP<ETAP)»< 1 ♦-<ETAP/TLAH)f$ALPT)tt2 
FUNC<FVfFT»FP>*FVt<FPf (l.-FP>/PHlO/(FT4(l,-FT>/TMR> 
FUNCO(FU#FT»FP)»Fy$< < 1 . -FT ) / < TMD4TH ) ) » ' FF4 ( 1 . -FP ) /PMD )/ < (FT4( I . - 


JET21310 
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JET21430 
JET2M40 
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JET21370 
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JET21630 
JET21A40 
JET2U70 
JET2I680 
JET21090 
FTJET21700 


l)/TM8)<f2) 

NPTS-24 

ALPT-ALPV 

WRITE < 6*777) TN*TM* ACPV*ALf T iCM.CT 

TM»-TH«TN 

0TM»<1.-TMR)/TM0 

DPH«(PHt“l . )/Ffip 

Nl-NPTS-1 

N2-NPTS-2 

XNl-Nl 

DETA-l ./XHl 

ETA-2ETA 

SUMOM-O. 

SUMOT-0. 

SCV2T0=g. 

SCTHO-0. 

DO !00 I*l*Nl*2 
ETAP-ETA 
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JET218l<I 

JET2102O 

JET21830 

JET21B40 

JET2t850 

JET21860 

JET21870 

JET21880 


ETAPT-ETAP/Tl.AM 

If (ETAPT.OT, i . > ETAPTrl . 

FV-FVP(ETAP) 

FT*FTP<ETAP> 

SUMOr iCM4FV»FUNC(FMtFTtFP> 

SUM0T=:-UM0T4FT*FUNC<FVfFT*FP) 

8CV2T0-SCV2T04FUNCP(FV*FT»FP)$FU 

SCTM0*5CTM04FUMCD(FV.FTfFP)4FT 

ETA-ETA42.FDETA 
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100 continue 

ETA«2.IDETA 

SUhEMrO. 

SUMET«0. 

SCV2T£»0. 

SCTHE-0. 

00 200 I»3»N2f2 

ETAP-ETA 

ETaPT«ETAP/TLA« 

IF <ETAPT.GT.l.) ETAPT«l. 

FV*FU?(ETAP) 

FT^FTP(ETAP) 

FP»FV$FV 

SU«E«-SU«£HiFV$FUNC<FV»FT»FP> 
SUMET»SUHETtFT*FUNC(Fy iFTfFP) 
SCU2TE«SCV2TEfFUNCP(FV»FT#FF >tFV 
SCTNC*SCTMEFFUNCO<FVrFT»FP)»FT 
ETA*£TA>2.tliETA 
200 CONTINUE 

CV2*<DETA/3, )$U ,44.*SUN0Hf2.*SUHEN) 
CT«(tiETA/3.>4U. ♦4.$SUMOT42. tSUhCT) 
CU2Th^(0ETA/3. U(4.«SCV2T0+2.tSCV2TE) 
CTTH»<tiETA/3. ) ♦ ( 4 . tSCTHO 42 . #SC THE ) 
CPF*l*-4./(ALPVfl. )46./i2.*ALPV+l . )-4. 
I ) 

CP*CPF 

RC0N*<2.'SQRT(2. >)/2. 

R*RCONttn ./ALPV> 

WRITE (6»277> TN » TH > ALPV i ALPT rCU2 » CP » C 
777 FORMAT (1X»7F10.5) 

RETURN 

END 
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Figure 69 - Concluded. 
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